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nucleation  radius  for  boiling 
vapor  core  radius 

resistance  to  liquid  flow,  Appendix  G 

universal  gas  constant 
Reynolds  number  for  liquid  flow 
Reynolds  number  for  vapor  flow 
gas  constant  for  vapor 
parameters,  Eq.  (23)'  (26)' 
thickness.  Chapter  III  and  Appendix  F 
time 

rise  time 

startup  time 

temperature 

room  temperature 

lumped  hot  /one  temperature 

rise  temperature 

startup  temperature 

initial  temperature 

test  chamber  shroud  temperature 


Symbol 

1v 

v 

w 

*g 

xv 

y 

z 
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SPECIAL  NOMENCLATURE  (continued) 

Definition 
vapor  temperature 
specific  volume,  Eq.  (73) 
groove  width 

mole  fraction  of  noncondensible  gas 
mole  fraction  of  vapor 

noncondensible  gas  concentration  variable,  Eq.  (12) 
distance  along  pipe 

dimensionless  z 


Greek  Symbol 

Yv 

« 

AP 

Ar 

AT 


eH 

eP 

eAB 


Definition 

vapor  specific  heat  ratio 
groove  depth 
pressure  drop 
radial  thickness 
temperature  difference 
emlssl vlty 

porosity  of  screen  wick,  Eq.  (49)  and  (83) 
emlssl vlty  of  hot  zone  Internally 
emlssivity  of  pipe  external  surface 
Lennard-Jones  parameter.  Appendix  A 
dimensionless  temperature,  Eq.  (23)'  (26)' 


o 


SPECIAL  NOMENCLATURE  (continued) 


Greek  Symbol 
e 

K 

X 

v 

v 

P 

a 

°AB 

al 

«0,AB 


Definition 

wetting  angle,  Eq.  (44) 

Boltzmann  constant,  Appendix  A 

latent  heat  of  vaporization.  Chapter  III 

parameter,  Eq.  (II) 

dynamic  viscosity 

kinematic  viscosity 

density 

Stefan -Boltzmann  constant 
collision  diameter 

surface  tension  coefficient,  Eq.  (40) 

heat  pipe  tilt  to  the  horizontal,  Fqs.  (41)  and  (4?) 

collision  Integral,  Appendix  A 


Subscript 

ax 

A 

Ar 

c 

C 

E 

eff 


Definition 

axial 

adiabatic 

argon 

cold  zone 

condenser 

evaporator 

effective 
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Subscript 

h 

1 

In 

I 

t 

m 

Na 

NCG 

o 

P 

r,R 

s 

ss 

v 

w 


Definition 

hot  zone 

inner 

input 

insert 

liquid 

mean 

sodium 

noncondensible  gas 
outer,  output 
pipe 

radial,  reservoir,  radiative  (as  noted  in  text) 
sonic,  screen 
stainless  steel 
vapor 

wick,  wire,  wall 


Abbreviation 

DWAHP 

dia . 

eq. 

ID 


Definition 

double  wall  artery  heat  pipe 

diameter 

equation 

inner  diameter 
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SPECIAL  NOMENCLATURE  (concluded) 


Abbreviation 

Definition 

LHS 

left  hand  side 

NCG 

noncondensible  gas 

00 

outer  diameter 

RHS 

right  hand  side 

TC# 

thermocouple  location 

number 

TIG 

tungsten-inert -gas 

VCHP 

variable  conductance 

heat  pipe 
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INIRODUCI ION 

1.1  About  the  Heat  Pipe  In  General 

The  heat  pipe  is  a  wonderful  and  conceptually  simple  two-phase 
(liquid  and  vapor)  thermodynamic/heat  transport  device  invented  to 
circumvent  the  barrier  of  solid  metallic  conductance.  The  heat  pipe 
working  principle  was  first  proposed  in  1944  but  the  fruit  of  the 
research  in  this  area  was  realized  only  during  the  past  couple  of 
decades  starting  from  the  original  work  on  heat  pipes  at  Los  Alamos 
National  laboratory,  which  coined  the  word  "heat  pipe"  [1-3]. 
Though  the  original  work  was  in  conjunction  with  the  thermionic 
converters  (to  work  out -of -core  in  nuclear  reactors),  heat  pipes 
were  used  extensively  In  spacecraft  for  thermal  control  applications 
in  the  pre  Apollo  and  post -Apollo  missions.  The  rapid  development 
of  this  device  caught  up  with  many  other  terrestrial  (solar,  energy 
regeneration  and  electronics  cooling)  applications  as  well.  Today 
this  device  is  a  commercial  off-the-shelf  item  with  dozens  of 
manufacturers  in  the  United  States,  Europe  and  Japan. 

The  name  heat  pipe  should  not  be  confused  with  the  "heat  pump* 
(which  has  a  separate  role  In  thermodynamics)  where  there  is  net 
work  done  whereas  in  a  heat  pipe  there  Is  no  net  work  and  no  moving 
parts.  Ihe  heat  pipe  is  a  closed  tube  or  chamber  with  one  of  many 
different  shapes  whose  inner  surfaces  are  lined  with  a  porous 
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capillary  wick.  The  wick  is  saturated  with  the  liquid  phase  of  the 
working  fluid  and  the  remaining  volume  of  the  chamber  contains  the 
vapor  phase.  Heat  applied  at  the  evaporator  vaporizes  the  liquid  in 
that  section  causing  the  pressure  difference  to  drive  the  vapor  down 
to  the  condenser  section  where  It  condenses  releasing  the  latent 
heat  of  vaporization  to  a  heat  sink.  Ihe  condensate  is  returned 
back  to  the  evaporator  by  the  capillary  pumping  action  of  the  wick 
and  this  process  continues  as  long  as  the  wick  does  not  become  dry. 

This  is  the  basic  working  principle  of  a  heat  pipe.  The  heat  pipe 

is  functionally  a  "heat  transport  pipe"  as  well  as: 

•  A  high  conductance  device.  It  surpasses  the  best  metallic 

conductors  such  as  gold,  silver  and  copper  in  thermal 

conductance.  Its  conductance  could  be  anywhere  from  30  4800 
times  that  of  an  equivalent  size  copper  bar  in  various 
temperature  ranges. 

•  Near-isothermal .  From  end-to-end,  the  temperature  drop  (AT) 

is  very  small;  much  smaller  than  that  of  a  metallic  conductor. 

•  Has  no  moving  parts.  It  is  a  hermetically  sealed  stationary 

device;  works  as  it  rests. 

•  Transports  heat.  From  one  end  to  the  other  end  without  the 

expense  of  any  external  work  done  on  it,  the  transport  mechanism 

being  the  closed  loop  flow  pattern  of  the  evaporating  liquid  and 
the  condensing  vapor  set  up  by  the  capillary  wick  and  the 

working  fluid  contained  inside  it. 
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•  Not  an  engine.  Since  it  does  not  have  any  net  work,  there  is  no 
efficiency  defined  for  this.  Its  performance  is  judged  by  the 
maximum  heat  transport  ability  and  the  distance  of  transport 
under  minimum  temperature  drop. 

•  Analogous  to  a  transformer  (electrical).  It  can  accept  heat 
input  at  certain  flux  and  deliver  at  a  different  flux. 

•  A  maintenance  free,  no-noise  device.  Once  installed,  it  does 
not  require  maintenance  and  does  not  make  any  noise  during 
operation. 

•  A  vapor  chamber.  It  may  build  up  high  pressure  Inside  if  heated 
inappropriately  leading  to  hazardous  explosion. 

Heat  pipes  can  be  categorized  into  several  types  based  on  their 

working  principle,  configuration  of  the  wick  structure  and  working 

temperature  range  [4]. 

(1)  Working  Principle:  •  Conventional  heat  pipe  mode 

•  Thermo  syphon 

•  Variable  conductance  heat  pipe 

•  Osmotic  heat  pipe 

•  Wickless  rotating  heat  pipe 

•  Capillary  pumped  loop  heat  pipe 

•  Electrodynamic  heat  pipe 

(2)  Conf iguration:  •  Shape  -  cylindrical,  flat,  annular, 

bent  in  forms  -  U,l,C,  etc. 

•  Structure  -  rigid  or  flexible 

•  Material  -  metallic  or  ceramic 


•  Wicks  conventional,  composite  and 
arteries 

(3)  Temperature:  •  Cryogenic  (<1?2  K) 

•  Moderate  temperature  (1?;’  6?8  K) 

•  High  temperature  (^6?8  K) 

Figure  1  illustrates  the  heat  pipe  temperature  regimes  ana  working 
fluids. 


Figure  1.  Heat  Pipe  Temperature  Regimes  and  Working  Fluids. 
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1 .2  Outstanding  Problem 

In  recent  years,  the  demand  for  high  capacity,  high 
performance,  and  reliable  heat  pipe  systems  has  Increased  with  the 
advent  of  newer  space  missions  such  as  space  stations  and 

terrestrial  applications.  Including  power  plants,  electronic 
cooling,  highway  deicing,  etc.  Fairly  well  established  design, 
fabrication,  material  compatibility,  and  steady  state  performance 
testing  procedures  for  all  types  of  heat  pipes  exist  In  the  industry 
[4  8].  However,  a  thorough  understanding  of  the  transient 
characteristics  Is  lacking.  An  area  of  major  concern  for  heat  pipe 
applications,  particularly  liquid  metal,  high  temperature  heal 
pipes.  Is  the  behavior  of  these  devices  under  transient  load 
conditions  during  startup,  shutdown  and  operational  variation  of  the 
thermal  loads.  The  startup  process  may  not  be  successful  when  the 
working  fluid  Is  In  the  frozen  state  since  there  will  be  no  liquid 
return  from  the  condenser  to  the  evaporator.  A  sudden  removal  of 

the  heat  input  to  the  evaporator  (shutdown)  without  removing  the 

condenser  load  may  freeze  the  working  fluid  at  the  condenser  and 

create  evaporator  dryout  [9],  These  transient  operating  problems 
are  further  complicated  by  special  artery  wick  designs  which  provide 
superior  performance  during  normal  operation.  An  example  of  a 
special  wick  design,  shown  In  Figure  2,  is  the  Double  Wall  Artery 
Heat  Pipe  (DWAHP)  developed  by  Ponnappan  and  Mahefkey  [10]. 

The  OWAHP  was  introduced  in  198?  and  a  modified  design  was 
developed  subsequently  [II].  In  the  modified  DWAHP,  a  screenless 
adiabatic  length  is  used  In  order  to  make  very  long  heat  pipes 
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Figure  2.  Double  Wall  Artery  Wick  Concept 


simpler.  Other  unique  advantages  of  this  design  are  entrainment 
free  transport  and  the  availability  of  lean  wick  evaporator 
configurations.  Though  this  design  concept  works  well  In  the  low 
temperature  regime  (25  200°C)  copper  water  system,  It  was  not  known 
whether  the  same  design  could  be  scaled  up  to  the  high  temperature 
regime  liquid  metal  system.  The  Issues  concerning  the  extension  of 
the  DWAHP  development  which  form  the  basis  of  the  present 
Investigation  are  as  follows: 

(1)  The  choice  of  the  wick  structure. 

(2)  The  methods  of  coupling  the  evaporator  Input  and  the  condenser 
output. 

(3)  The  thermal  load  requirements. 

(4)  The  Initial  conditions  of  thawing  or  priming  status  of  arteries. 

(5)  Whether  a  noncondensible  gas  filling  will  help  startup  from 
frozen  state. 

(6)  What  effects  the  gas  will  have  on  artery  priming. 

(7)  The  modeling  of  the  startup  from  frozen  state  using  the  coupled 
diffusion  and  energy  balance  equations. 

(8)  The  design,  fabrication  and  testing  of  an  appropriate  hardware 
to  verify  the  theoretical  modeling. 

1 . 3  Review  of  Literature 

A  literature  search  was  conducted  and  a  brief  review  is 
presented  here.  Though  the  heat  pipe  technology  is  only  25  years 
old,  an  abundant  number  of  technical  publications  covering  the  broad 
spectrum  of  this  technology  exists  In  the  open  literature.  A  closer 
scrutiny  of  works  related  to  transient  behavior,  unsteady  power, 
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pulsed  operation,  startup,  shutdown,  and  analysis/modeling  narrowed 
down  the  publications  to  a  sizeable  number  covering  the  period 
between  1967  and  1988.  Online  Information  Retrieval  (computer 
search)  of  data  bases  covering  Physics  Abstracts,  Aerospace 
Abstracts,  NTIS,  OTIC,  etc.  was  used  for  the  literature  search. 

Cotter  [12]  was  the  first  to  quantify  the  dynamics  of  a  heat 
pipe  startup  using  a  lumped-parameter  one -dimensional  model.  An 
effective  thermal  conductivity  In  the  axial  direction  of  the  heat 

pipe  dependent  upon  the  state  of  the  vapor  for  the  free  molecular 
flow  regime  was  introduced.  Characteristic  response  times 
(condenser  rate  an:l  axial  rate  of  responses)  were  defined.  He 
predicted  the  required  length  of  the  heat  pipe,  once  the  Input  and 
output  conditions  were  determined. 

Sockol  and  Forman  [13]  visually  observed  the  moving 

temperature  front  during  startup  in  a  lithium  heat  pipe 
(1273-1673  K).  Cotter's  model  was  reassessed  and  modified.  They 
predicted  and  confirmed  that  the  temperature  of  the  hot  zone  was 
fairly  independent  of  the  power  Input  to  the  evaporator  until  the 
hot  zone  reached  the  end  of  the  pipe. 

Marcus  [14,15]  conducted  an  extensive  investigation  on  the 
theory  and  design  of  room  temperature  gas  loaded  variable 

conductance  heat  pipe  (VCHP).  His  study  concentrated  mainly  on  the 
steady  state  theory.  Marcus  noted  that  "seeking  an  accurate 

predictive  capability  for  heat  pipe  transient  behavior  is 
substantially  a  more  ambitious  task  than  the  development  of  the 
steady  state  theory." 
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Tower  [16]  performed  an  analysis  of  oxygen  diffusion  Into 
metal  alloys  pertaining  to  the  startup  of  T-lll  (tantalum  alloy) 
heat  pipe  filled  with  lithium. 

Sockol  [17]  presented  a  model  for  the  rapid  startup  of  a  high 
temperature  gas  loaded  heat  pipe.  A  two-dimensional  analysis  was 
used.  Startup  was  not  Initiated  until  the  vapor  pressure  Py  In 
the  hot  zone  reached  a  value  proportional  to  the  Initial  gas 
pressure  .  Through  proper  choice  of  ,  startup  could  be 
delayed  until  P^  was  large  enough  to  match  with  the  thermal  load 
on  the  heat  pipe.  However,  the  analytical  results  were  not  verified 
with  any  experimental  data. 

The  use  of  the  liquid-trap  as  a  secondary  heat  pipe  for 
forward  mode  operation  during  a  heat  pipe  diode  shutdown  was  studied 
by  B&K  Engineering  [18].  Thermal  switching  was  defined  as  the 
shutting  down  of  the  thermal  diode  In  case  of  a  reversal  of  heat 
flow  direction. 

Camarda  [19]  used  a  simple  analytical  technique  to  determine 
startup,  transient  and  steady  state  performance  of  the  heat  pipe 
during  testing.  The  temperature  ranged  from  883  K  to  9 22  K  and  heat 
flux  23.9  to  39.5  W/cm2. 

Colwell  [20]  developed  a  method  for  predicting  the  behavior  of 
cryogenic  heat  pipes  during  startup  or  changes  In  thermal  transport. 

Williams  [21,25)  evaluated  the  forward  (heat  flow  from 
evaporator  to  condenser)  and  reverse-mode  (heat  flow  from  condenser 
to  evaporator)  operation  for  an  ethane  heat  pipe  working  in  170  to 
220  K.  Results  Indicated  that  the  heat  pipe  might  not  reliably 
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startup  In  the  forward  mode.  However,  startup  could  be  initiated 
when  preceded  by  a  diode  reversal. 

Groll  and  Munzel  [22,23,2b]  presented  an  analysis  and  design 
of  all  aluminum  and  all  stainless  steel  ammonia  heat  pipe  diode. 
The  liquid  trap  principle  was  used.  The  shutdown  ratio,  defined  as 
the  ratio  of  the  forward  mode  conductance  to  the  reverse  mode 
conductance  (W/K),  was  384.  Test  program  comprised  transient 
shutdown  behavior  and  restart  from  the  shutdown  condition. 

Alario  [24]  developed  a  stainless  steel  ethane  heat  pipe  diode 
and  used  it  for  investigating  the  startup  dynamics  of  arterial 
cryogenic  thermal  diodes. 

Shukla  [27]  conducted  a  combined  analytical  and  experimental 
investigation  of  the  transient  response  of  a  gas  controlled  heat 
pipe  with  an  absorption  gas  reservoir. 

Groll  and  Supper  [28]  showed  a  shutdown  technique  which 
consisted  of  heating  up  the  condenser  with  a  constant  heat  flux 
while  having  the  evaporator  and  trap  thermally  free-floating. 

Bystrov  et  al.  [29]  presented  the  results  of  experimental 
studies  of  the  startup  of  vacuum  and  of  gas  loaded  sodium  heat  pipes 
by  sudden  and  smooth  application  of  heat  to  the  evaporator.  A 
qualitative  analysis  of  the  dynamics  of  startup  of  such  pipes  was 
also  presented. 

Bystrov  and  Goncharov  [30]  examined  an  approximate  method  of 
calculating  the  unsteady  temperature  fields  and  the  startup 
characteristics  for  a  quick  startup  and  switch  to  the  heat  pipe 
regime,  using  a  lumped  parameter  model  for  the  vapor  predicted  on  a 
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planar  vapor -gas  boundary  and  the  possibility  of  neglecting  the 
thermal  inertia  of  the  vapor  flux.  The  analytical  results  were  in 
good  agreement  with  the  experimental  results  of  an  argon  filled 
annular  artery  sodium  heat  pipe. 

Colwell  and  Chang  [31]  examined  the  transient  behavior  of  a 
slab  type  wick  In  a  low  temperature  heat  pipe.  It  was  found  that 
under  some  conditions,  the  heat  pipe  might  approach  steady  state 
operation  with  a  portion  of  the  capillary  structure  dried  and  that 
rewetting  could  be  quickly  accomplished.  However,  when  the  entire 
wick  structure  at  the  evaporator  dried,  rewetting  was  much  more 
difficult. 

Chang  and  Colwell  [32]  developed  a  computational  model  based 
on  finite-difference  approximations  for  predicting  the  transient 
operating  characteristics  of  low  temperature  heat  pipes.  It  was 
shown  that  the  width  of  fluid  gap  between  screen  layers  had  a 

significant  effect  on  the  transient  operation. 

Colwell  and  Hartley  [33]  emphasized  the  need  for  developing 
the  capability  to  predict  operational  behavior  of  liquid  metal  heat 
pipes  under  adverse  operating  conditions  such  as  when  drying, 
rewetting,  choking,  noncontinuum  flow,  freezing,  thawing,  etc., 
occurred  within  the  heat  pipe. 

Beam  [34]  developed  a  lumped -parameter  transient  model  of  a 
heat  pipe  and  verified  his  model  with  data  obtained  from  a  screen 

wick  copper -water  heat  pipe.  Two  response  times,  one  corresponding 

to  the  thermal  response  of  the  heat  pipe  and  the  other  to  the 

hydrodynamics  of  the  wick  with  an  applied  heat  flux,  were  defined. 
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An  upper  limit  for  the  pulse  size  (so  as  not  to  create  a  dryout)  was 
also  determined  In  this  analysis. 

Culllmore  [35]  described  the  Advanced  Simulation  Analysis  of 
Heat  Pipes  (ASAHP)  program,  that  simulated  the  transient  operation 
of  heat  pipes.  The  transient  capabilities  Included  wicking  and 
entrainment-limited  dryout  and  recovery  simulation  with  the 
corresponding  condenser  flooding.  The  program  provided  a  means  to 
evaluate  the  effects  of  adverse  transients  (rapid  changes  in  orbital 
environments,  or  high -power/short -duty  cycles)  on  a  spacecraft  heat 
pipe  radiator. 

Ambrose  et  al.  [36]  have  presented  a  simplified  one 
dimensional  model  to  predict  the  drying  and  rewetting  of  the  wick 
during  a  pulsed  startup.  Their  model  Included  a  thermal  capacitance 
term.  They  compared  model  predictions  with  the  experimental  results 
obtained  with  a  copper  water  screen  wick  heat  pipe.  It  was  verified 
that  the  ability  of  the  heat  pipe  to  transport  a  pulsed  heat  load 
was  equal  to  the  maximum  design  capillary  heat  transport  rate 
without  dryout. 

Tilton  et  al.  [37]  analyzed  the  transient  behavior  of  a  steady 
operating  Inconel  617  sodium  heat  pipe  with  superimposed  external 
thermal  loading  at  the  condenser.  The  analytical  model  predicted 
the  direction  of  the  mass  flow  of  vapor  and  the  location  of  the 
condensation  and  evaporation  zones  along  the  heat  pipe  as  a  function 
of  time  for  different  condenser  heat  loads.  Heat  pipe  flow  reversal 
and  Incipient  dryout  could  be  predicted.  Experimental  results  of  a 
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45-cm  long  1.27  cm  diameter  screen  wick  heat  pipe  using  radiant 
heaters  were  used  to  verify  the  analytical  prediction. 

Merrlgan  et  al.  [38]  reported  transient  performance  character 
Istics  of  a  4  m  moly-llthlum  heat  pipe  with  annular  screen  wick. 
Shutdown  tests  were  conducted  to  establish  the  requirements  of 
restarting  the  pipe  from  frozen  condition.  A  freeze  plug  could  form 
In  the  condenser  during  the  cooldown  and  cause  liquid  depletion  of 
the  evaporator  thereby  causing  restart  difficulty. 

Bobco  [39]  developed  a  closed  form  expression  for  estimating 
the  location  of  the  minimum  temperature  in  the  Inactive  condenser  of 
a  VCHP  under  steady  state.  Equations  were  presented  to  show  how  the 
gas  reservoir  design  may  be  based  on  the  operating  conditions  for 
maximum  heat  load  and  Incipient  freezeout.  The  experimental  data  of 
Marcus  [14]  were  used  to  verify  the  theory. 

Costello  et  al.  [40]  developed  a  comprehensive  model  but  did 

not  report  calculation  results.  Hall  and  Ooster  [41]  Introduced  the 

transient  computer  code  called  THROPUT  (Thermal  Hydraulic  Response 

of  Heat  Pipe  Under  Transients)  and  applied  It  to  verify  the 

performance  of  a  4  m  lithium  heat  pipe  charged  with  a  noncondensible 

-4 

gas  at  10  Pa  and  300  K.  This  code  required  large  CPU  time  (24 

hours  on  a  VAX  8000  system  for  2  hours  of  transients).  In  addition. 
It  Is  reported  that  the  liquid  does  not  return  fast  enough  to  keep 
up  with  evaporation,  resulting  in  an  ever  decreasing  liquid 

Inventory  In  the  evaporation  section. 

Jang  et  al.  [42]  used  finite  element  techniques  with  explicit 
and  Implicit  solution  procedures  to  solve  the  mathematical  model  of 
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a  liquid  metal  heat  pipe  used  In  leading  edge  cooling.  They 

compared  model  predictions  with  experimental  data  of  Camarda  [19] 
and  found  good  agreement. 

«  / 

Seo  and  El-Genk  [43]  developed  a  two  dimensional  transient 

model  for  liquid  metal  heat  pipes  that  Is  capable  of  predicting  the 

operating  limits  of  the  heat  pipe  during  steady  state  and  transient 

operations.  Their  model  assumes  a  completely  thawed  heat  pipe  with 

no  noncondensible  gas  In  the  vapor  core  of  the  annular  wick  heat 

pipe.  They  also  show  how  the  slow  cooling  rate  (0.1  kw/s)  or  fast 

heat  rate  (1.0  kW/s)  from  a  stable  operating  point  can  encounter  the 

sonic  or  entrainment  limit  respectively  In  a  4  m  long  1.9  cm  outside 

diameter  lithium  heat  pipe  built  at  Los  Alamos  National  Laboratory 

..Vf 

[38]. 

Peterson  and  Tien  [44,45]  found  numerical  and  analytical 
solutions  for  two -dimensional. '  gas  distribution  In  gas  loaded  heat 
pipe  in  the  room  temperature  regime.  No  frozen  startup  or  transient 
behavior  is  stud.ied  here.  The  authors  claim  that  the  errors  In 
"flat -front"  diffusion  models  are  eliminated  and  that  the  condenser 
shutoff  length  Is  computed  accurately. 

It  is  concluded  from  the  literature  review  that  a  combination 
of  analytical  and  experimental  works  on  the  transient  behavior  of 
various  types  of  heat  pipes  exists.  Most  of  the  referenced  papers 
deal  with  thermal  diodes  and  cryogenic  or  room  temperature  heat 
pipes.  There  are  only  three  publications  on  the  gas  filled  liquid 
metal  heat  pipe  transient  behavior;  those  of  Sockol  [17]  and  Bystrov 
et  al.  [29,30] . 
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It  Is  difficult  to  apply  the  analysis  pertaining  to  a 
particular  heat  pipe  design  to  another  constructionally  different 
heat  pipe.  Conventional  or  the  homogeneous  wick  heat  pipes  are  no 
longer  useful  for  high  capacity  application.  The  state-of-the-art 
heat  pipe  designs  have  developmental  and  material  problems.  Making 
the  heat  pipes  suitable  for  adverse  operating  conditions  is  a 
technological  challenge.  In  view  of  this,  It  is  of  Interest  to 
develop  the  Oouble  Wall  Artery  Heat  Pipe  into  a  liquid  metal  high 
capacity  design  and  evaluate  its  performance  theoretically  and 
experimentally  for  steady  and  transient  thermal  loads. 

1 . 4  Need  and  the  Uniqueness  of  the  Research 

Future  space  missions  will  require  large  area  space  radiators 

constructed  out  of  liquid  metal  heat  pipes  to  reject  large  amounts 

of  waste  heat  at  rejection  temperatures  of  1000  K  or  more.  Peak 

electric  power.  In  excess  of  1  MWe,  is  required  which  generates 

waste  heat  in  the  form  of  pulses  with  a  peak -to  average  ratio  of 
4  S 

10  to  10  based  on  the  mission  and  orbital  duty  cycles 
envisioned  [46.47].  Fssentially  these  heat  pipes  have  to  withstand 
directed  energy  Impacts  such  as  micrometeoroids,  or  pulsed  input 
loads.  A  space  nuclear  power  system  Is  a  typical  example  of  a  power 
supply  requiring  such  high  capacity  liquid  metal  heat  pipes  [48]. 

Another  example  is  NASA's  space  station  which  Is  being 
designed  to  be  powered  by  a  solar  dynamic  cycle.  A  solar  dynamic 
system  will  require  thermal  energy  storage  and  some  form  of  energy 
transport  device  (liquid  metal  heat  pipe)  for  transporting  energy  to 
the  power  conversion  unit.  The  transition  from  illuminated  to 
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eclipse  portion  of  the  orbit  will  introduce  transient  thermal  loads 
on  the  heat  pipe  device  [49].  In  addition,  cold  start  or  thawing  of 
a  liquid  metal  heat  pipe  is  a  unique  problem  for  space  or 
terrestrial  applications.  Uniform  distribution  of  the  working  fluid 
within  the  heat  pipe  wick  structure  Is  Important  during  shutdown  of 
any  heat  pipe  device.  An  abrupt  shutdown  during  a  dryout  may  cause 
irrecoverable  damage  [9].  Long  transport  lengths  cause  priming  and 
flow  resistance  problems.  In  essence,  the  need  of  a  reliable  liquid 
metal  heat  pipe  Is  apparent  from  the  above  applications,  and  the 
uniqueness  of  Investigating  a  liquid  metal  DWAHP  for  Its  startup  and 
shutdown  characteristics  is  exemplified  by  the  operational 
requirements  shown  in  Figure  3. 

Though  the  general  problem  of  transient  operation  is  common  to 
all  the  liquid  metal  heat  pipes  as  well  as  low  temperature  heat 
pipes,  the  same  problem  Is  unique  to  the  liquid  metal  DWAHP  proposed 
to  be  investigated  presently  for  the  following  reasons: 

(1)  Long  transport  section  length  without  screen  wick,  no  direct 
contact  between  liquid  and  vapor  flows  in  the  transport 
sections ; 

(2)  Composite  wlcklng  arrangement  with  arterial  grooves 

nonconstant  groove  widths;  and 

(3)  Noncondensible  gas  loading  and  reservoir  wick  in  evaporator  to 
aid  startup  from  frozen  state. 

1 . 5  Statement  of  the  Problem 

The  startup  of  a  liquid  metal  heat  pipe  (sodium,  potasium  or 
lithium)  from  the  frozen  state  Is  an  area  of  concern  for  potential 
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applications  in  space  radiators  since  the  heat  pipe  may  dry  out  due 
to  the  working  fluid  depletion  in  the  evaporator  during  startup. 
Though  active  methods  such  as  electrical  trace  heating  or  prelaunch 
thawing  can  be  adapted  for  starting,  they  are  neither  convenient  nor 
desirable.  A  passive  and  reliable  startup  technique  which  will  not 
require  onboard  electrical  power  is  preferable.  Researchers  have 
recommended  filling  the  liquid  metal  heat  pipe  with  noncondensible 
gas  in  order  to  make  It  start  easily  [17,30].  Ihe  heat  pipe 
geometry  they  used  was  quite  different  from  the  one  tested  in  the 
present  work.  A  composite  wick  high  capacity  liquid  metal  heat  pipe 
of  the  DWAHP  type  with  long  adiabatic  section  has  not  been  built  and 
tested  so  far. 

The  DWAHP  gains  further  attention  for  Investigation  because  of 
its  entrainment  free  and  capillary  wickless  transport  section,  high 

radial  heat  flux  evaporator  d  ign,  and  ease  of  manufacture.  The 

present  problem  Is  to  Investigate  and  find  answers  to  the  following 
basic  questions  on  this  subject  area. 

(1)  Will  the  noncondensible  gas  charging  provide  reliable  startup 
of  a  liquid  metal  arterial  heat  pipe  of  the  DWAHP  type? 

(2)  What  analytical  methods  could  be  c:scd  or  need  to  be  developed 

to  solve  the  transient  problem  of  the  above  mentioned  heat 

pipe,  e.g.,  the  mass  diffusion  and  thermal  problems? 

(3)  How  do  the  vapor  front  propagation  and  the  temperature  vary 

with  time  and  input  power  variation? 

(4)  What  experimental  methods  can  be  devised  to  explore  the  above? 
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1 . 6  Objectives  of  the  Present  Study 

The  scope  and  objectives  of  the  present  study  are  grouped  and 
defined  under  two  categories,  namely,  analytical  and  experimental: 

•  Analytical .  Develop  an  analytical  model  and  solve  to  predict 

the  frozen  state  startup  and  the  transient  thermal 

characteristics  of  a  high  temperature  heat  pipe  with  long 
adiabatic  section  and  filled  with  noncondensible  gas  charge. 
The  study  of  the  effects  of  different  initial  charge  pressures 
and  power  input  is  also  of  interest. 

*  Experimental .  Design,  fabricate  and  test  suitable  experimental 

hardware  to  verify  the  analytical  predictions  on  the  startup 

transients.  Investigate  the  dynamics  of  the  diffusion  zone 
during  the  transient  condition  and  determine  the  condenser 
shutoff  length  during  the  steady  state  condition. 
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CHAPTER  II 


THEORETICAL  ANALYSIS  OE  STARTUP 
2 .  T  Startup  Characterise c s 

In  general,  heat  pipes  are  intended  for  isothermal  operation 
at  a  fixed  temperature  and  energy  transport.  Starting  up  a  heat 
pipe  does  not  present  any  difficulty  if  it  is  heated  isothermal  ly  to 
the  operating  temperature  with  the  power  then  gradually  increased  to 
the  allowable  limit.  Often  in  practical  situations,  this  method  of 
startup  is  inconvenient  or  even  Impossible  and  also  the  condenser 
zone  heat  removal  is  effected  simultaneously  with  the  evaporator 
heating.  Other  aspects  which  concern  startup  are: 

(T)  The  allowable  rate  and  duration  of  heat  supply; 

(2)  The  method  of  heat  transfer  coupling  between  the  source  and  the 
evaporator,  and  the  sink  and  the  condenser; 

(3)  The  initial  state  (solid,  liquid  or  vapor)  of  the  working  fluid 
within  the  wick  structure; 

(4)  The  heat  transport  capacity  versus  temperature  charac¬ 
teristics  of  the  heat  pipe;  and 

(5)  The  wick  design  in  the  transport  zone. 

Three  characteristic  startup  types  have  been  categorized  for 
heat  pipes: 

(1)  Uniform  startup  (e.g.,  water  and  ammonia  heat  pipes); 
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(2)  Frontal  startup  without  noncondensible  gas  (NC6)  (e.g.,  sodium 
heat  pipe) ;  and 

(3)  Frontal  startup  with  NCG  (e.g.,  sodium  heat  pipe  with  argon). 
These  are  very  well  described  in  the  heat  pipe  literature  [4,8]. 
Figure  4  illustrates  these  startup  characteristics  schematically. 

Heat  pipes  with  working  fluids  having  high  vapor  pressures  at 
startup  temperatures  typically  exhibit  uniform  startup,  while  those 
with  low  vapor  pressures  exhibit  frontal  startup.  Use  of  NCG  In  the 
vapor  core  has  been  found  to  ease  the  startup  difficulties  In  liquid 
metal  heat  pipes  [8,17,30].  There  are  however  advantages  and  dis¬ 
advantages  In  charging  a  heat  pipe  with  NCG: 

Advantages 

(1)  The  heat  pipe  is  enabled  to  operate  as  a  variable  conductance 
heat  pipe  (VCHP). 

(2)  The  NCG  acts  as  a  buffer  to  the  working  fluid  vapor  while 

expanding  Into  the  vapor  core  and  helps  the  starting  or 

shutting  down  of  the  heat  pipe. 

(3)  The  size  of  the  gas  slug  In  the  heat  pipe  at  a  given 

temperature  can  be  predetermined  and  this  flexibility  allows 
the  designer  to  control  the  length  of  Inactive  condenser. 

(4)  The  need  for  special  starting  up/shutting  down  arrangements 

such  as  trace  heaters  and  condenser  decoupling  from  a  heat  sink 
is  eliminated. 
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Figure  4.  Heat  Pipe  Startup  Characteristics 
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Disadvantages 

(1)  lhere  Is  a  potential  danger  of  NC6  bubble  trapment  In  artery 
wicks  which  may  lead  to  depriming  of  the  artery. 

(2)  The  choices  of  compatible  fluid  and  NCG  pair  may  be  limited  for 
high  temperature  applications.  The  pair  should  be  chemically 
non-reacting  and  non  absorbing  (non -soluble)  types  even  under 
high  pressure  and  temperature  conditions. 

(3)  Even  with  careful  design,  the  non-use  of  a  portion  of  the 
condenser  (inactive  length)  under  steady  state  operation  may 
have  to  be  tolerated. 

Stages  of  Frozen  State  Startup 

A  gas  filled  heat  pipe  in  general  (or  the  present  OWAHP  in 
particular)  undergoes  several  stages  of  energy  transport  during  the 
startup  from  the  Initial  frozen  state  to  the  operating  regime  [4,30]. 
Stage  1 :  Melting  and  evaporation  of  heat  transfer  fluid. 

With  the  supply  of  uniform  heat  Input  to  the  evaporator, 
the  frozen  working  fluid  thaws  and  eventually  evaporates  at  the 
liquid  meniscus  in  the  capillary  pores  of  the  wick.  Gradual 
withdrawal  of  NCG  from  the  evaporator  to  the  adiabatic  section 
begins. 

Stage  2 :  F ormation  of  axial  vapor  f lux . 

As  the  heat  power  supplied  Increases,  the  temperature  of 
the  evaporator  increases  and  the  vapor  pressure  compresses  the 
NCG  away  from  the  evaporator.  The  speed  of  vapor  flow  at  the 
evaporator  exit  Increases  up  to  sonic  speed. 
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Stage  3 :  Sonic  regime . 

Continuation  of  power  Input  to  the  evaporator  takes  the 
heat  pipe  along  the  boundary  of  the  sonic  limit.  In  this 

regime,  considerable  non -isothermal  conditions  prevail  along 
the  heat  pipe.  Evaporator  fluid  inventory  continues  to  deplete 
as  the  adiabatic  section.  If  any,  thaws  and  no  fluid 

recirculation  from  condenser  commences. 

Stage  4:  Sonic  to  supersonic  vapor  flow. 

With  Increase  In  vapor  pressure,  the  vapor  gas  front 

moves  toward  the  condenser  end.  Vapor  velocity  variations 
result  from  a  variable  mass  flow  through  a  constant  vapor  core 
area.  The  velocity  can  Increase  further  and  become  supersonic, 
or  a  recompression  of  the  vapor  may  occur,  resulting  in  a 

pressure  recovery.  The  amount  of  pressure  recovery  depends  on 
the  heat  removal  rate  In  the  condenser  section. 

Stage  5:  Isothermal  operating  regime. 

Once  the  vapor  flow  returns  to  the  subsonic  levels,  the 
liquid  return  flow  establishes  recirculation,  and  the  condenser 
heat  removal  rate  equals  the  evaporator  input  rate,  the  pipe 
switches  to  steady  Isothermal  operating  condition. 

2.2  Parametric  Analysis 

Homogeneous  and  annular  wick  liquid  metal  heat  pipes  are 
reported  to  have  been  analyzed  or  tested  In  NCG  filled  mode  without 
problems  [17,29,30].  However,  it  Is  not  yet  established  in  the 
liquid  metal  heat  pipe  research  literature  that  if  a  grooved  artery 
heat  pipe  (of  the  DWAHP  type  with  unvented  artery  and  screenless 
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adiabatic  section)  would  start  from  frozen  state  and  work  with  NCG 
filled  mode.  In  order  to  answer  this  question,  a  preliminary 
analysis  addressing  the  design  and  testing  parameters  Is  done  In  the 
following  question  and  answer  format. 

(1)  Whether  freezing  of  sodium  In  a  closed  artery  channel  will 
create  any  problem? 

Sodium  shrinks  on  freezing  since  the  density  of  the  solid  Is  larger 
than  that  of  the  liquid  at  the  freezing  point  (97.83°C),  while  the 
volume  change  Is  only  2.63%.  The  void  created  within  the  artery 
channel  (which  is  fully  primed  during  normal  operation)  may  be 
vacuum  or  NCG  depending  upon  the  propagation  of  the  solidification 
front  of  the  sodium  during  shutdown  of  the  heat  pipe.  On  thawing 
(Initial  startup),  the  void  will  be  filled  back  If  vacuum  or  the  gas 
will  be  expelled  out  of  the  artery  if  NCG  was  present  to  start  with. 

(2)  Whether  the  problem  of  vapor  back-flow  Into  the  artery  exists 
In  the  liquid  metal  OWAHP? 

This  Is  not  likely  under  low  radial  evaporator  heat  fluxes 
2 

( <50  W/cm  )  since  the  critical  superheat  (aTcr^)  required 

to  create  bubble  nucleatlon/boi 1 Ing  in  the  wick  Is  approximately 

40  C  at  700°C  (see  Table  11).  On  the  other  hand,  the  evaporator 

at  across  the  sodium  saturated  stainless  steel  wick  (with  typical 

equivalent  thermal  conductivity,  kp  =  40  W/m°C)  Is  about  3  C  for  a 

heat  Input  of  1797  W  at  800°C  (see  Section  3.5.5).  At  higher  heat 

2 

fluxes  (>50  W/cm  ),  capillary  flow  Inserts  at  the 
evaporator/adiabatic  transition  may  be  needed  to  stop  vapor  back 
flow  (see  Reference  55  for  details). 
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(3)  What  are  the  conditions  under  which  the  NCG  bubbles  can  get 
Into  the  artery?  Under  what  conditions  are  they  deleterious 
to  the  heat  pipe  operation? 

This  depends  on  the  solubility  of  the  NCG  Into  the  working  fluid  and 
the  construction  of  the  artery  at  the  condenser.  If  solubility  of 
NCG  in  the  working  fluid  is  high,  the  gas  will  tend  to  disperse  in 
the  fluid  and  stay  under  subcooled  conditions  within  the  artery.  If 
the  fluid  is  frozen  In  this  condition.  It  may  lead  to  bubble  growth 
problems  later  while  thawing.  In  the  present  case,  argon  is  not 

soluble  in  sodium  since  the  Ostwald  coefficient  is  low 

(a  =  2.27x10  at  500°C)  [50]. 

If  the  construction  of  the  artery  at  the  condenser  section  is 
such  that  the  artery  wall  is  a  fine  porous  structure  (such  as  a 
screen  or  sinter),  the  vapor -gas  mixture  may  pass  through  the  porous 
wall  and  get  trapped  due  to  the  capillary  pressure  difference  across 
the  artery  wall.  In  the  case  of  vented  artery  channels  (such  as  in 
the  OWAHP  condenser),  the  liquld/vapor  meniscus  exists  Inside  the 
groove  at  each  of  the  vent  slots  which  are  large  enough  to  let  the 
gas  vent  back  into  the  vapor  core  while  the  condensation  of  the 

vapor  takes  place.  Hence,  conceptually,  gas  bubbles  cannot  be 
trapped  in  the  artery  channels  of  the  present  design. 

(4)  In  case  the  NCG  bubble  Is  trapped  in  the  closed  artery 

channel,  will  It  affect  the  priming?  Will  the  bubble  grow  and 
prevent  startup? 

Upon  thawing  of  the  sodium,  the  volume  of  the  bubble  will  expand.  A 
further  Increase  in  temperature  could  expand  the  bubble  more  and 
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might  cause  artery  depriming.  Hence,  under  this  hypothesis,  the 
heat  pipe  may  not  start.  The  other  possibility  is  that  the  bubble 
may  be  Initially  stagnant  at  an  axial  position  during  the  thawing 
process,  but  start  moving  toward  the  evaporator  (the  direction  of 
higher  capillary  potential)  soon  after  the  moving  hot  front  croiscs 


the  adiabatic  zone 

and 

the 

heat  pipe 

action  starts. 

Both 

possibilities  exist 

in 

the 

practical 

situation  and 

hence 

experimental  study  Is  necessary.  In  order  to  avoid  the  multiple 
artery  averaging  effects  from  confusing  the  experimental  results,  it 
was  planned  to  construct  a  single  channel  artery  heat  pipe  as  shown 
in  Figure  5.  The  details  of  this  design  are  given  in  Chapter  III. 

(5)  During  startup  mode,  will  there  be  condensate  formation  at  the 
upstream  side  of  the  moving  temperature  front?  If  so,  what 
are  the  effects  on  the  startup? 

This  Is  a  unique  situation  with  respect  to  the  ventless  adiabatic 
artery  channels  as  opposed  to  a  porous  (sinter  or  screen  wall) 
artery.  The  vapor  traveling  down  the  vapor  core  will  lose  its 
latent  heat  to  the  incremental  length  of  the  adiabatic  section  as 
the  vapor/gas  front  moves.  This  movement  depends  on  the  axial 
conduction  and  the  diffusion  rates  of  the  vapor  into  the  NCG. 
Initially,  there  will  be  some  amount  of  condensate  collection  at  the 
vapor/gas  front  but  as  more  and  more  vapor  is  generated  and  the 
vapor  pressure  builds  up,  the  already  formed  droplets  will  be 
vaporized  back  and  moved  toward  the  advancing  front.  The  thawing 
proceeds  progressively  and  the  front  moves  down  the  vapor  core.  As 
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OUTER  TUBE 


Figure  5.  Single  Channel  Artery  Heat  Pipe  Verflcatlon  Model 


the  front  reaches  the  condenser  vent  slot  zone,  the  heat  pipe  action 
starts . 

The  key  analytical  design  issue  here  is  to  satisfy  the 
requirement  that  the  evaporator  should  not  be  depleted  of  working 
fluid  before  the  thawing  of  the  transport  section  is  completed.  The 
diffusion  gas  front  analysis  for  this  situation  will  determine  the 
maximum  possible  adiabatic  length  for  the  case  of  frozen  startup.  In 
addition  to  the  axial  conduction  (along  the  heat  pipe  wall)  and  the 
diffusion  mechanism  available  for  the  energy  transfer  during  start¬ 
up,  there  also  Is  the  radiation  mode  to  help  the  thawing  process. 

2 • 3  Heat  Trans fer  Model  of  the  Startup  Process 

A  gas  loaded  heat  pipe  contains  a  fixed  quantity  of  NCG  in 
addition  to  the  working  fluid.  The  purpose  of  the  gas  Is  to 
facilitate  the  startup  from  frozen  state.  All  the  three  modes  of 
heat  transfer,  namely,  conduction,  convection  and  radiation,  occur 
during  the  startup.  Identifying  and  estimating  the  dominant  mode  is 
the  key  to  the  solution  of  the  startup  analysis.  Figure  6  shows  the 
heat  transfer  paths  of  a  gas-filled  heat  pipe  during  startup.  The 
heat  input  (Q^)  to  the  evaporator  Is  transferred  to  the  adiabatic 
and  condenser  sections  by  1)  axial  conduction  (0  )  through  the 

dX 

wall  and  wick,  2)  convective  transfer  by  hot  vapor  diffusing  into 
the  NCG  and  condensing  at  the  downstream  side  of  the  diffusion  front 
(Qv),  and  3)  radiation  exchange  between  the  hot  and  cold  zones 

(Qr). 

Generally,  In  the  moderate  temperature  (122  628  K)  VCHP 
research,  has  been  reported  to  be  much  smaller  than  the  axial 
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TEMPERATURE 


Figure  6.  Heat  Transfer  Model  for  a  Gas  Loaded 
Heat  Pipe  During  Startup 
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conduction  (Q  ),  with  the  diffusion  rate  assumed  constant 
[4,14,39].  For  high  temperature  heat  pipes,  Sockol  [17]  found  that 
Qv  Is  much  larger  than  and  that  approaches  at  the 

beginning  of  the  startup  and  drops  to  almost  zero  when  steady  state 
Is  approached.  This  theory  justifies  the  stationary  vapor-gas  front 
at  steady  state.  On  the  contrary,  if  axial  conduction  term  were 
dominant,  the  hot  front  may  move  even  after  the  diffusion  process 
stopped.  High  wall  conductance  tends  to  spread  the  front  over  the 
condenser  masking  the  diffusion  zone  temperature  profile.  High 
temperature  wall  materials  such  as  stainless  steel  and  refractory 
metals  have  poor  thermal  conductivity;  hence,  the  temperature 
gradient  at  the  diffusion  front  will  be  sharper  than  that  for  high 
conductivity  materials  such  as  copper  and  aluminum  commonly  used  for 
low  temperature  heat  pipes. 

A  qualitative  picture  of  the  temperature  profiles  at  the 
diffusion  zone  Is  shown  In  conjunction  with  the  heat  transfer  model 
sketched  in  Figure  6.  At  the  beginning  of  the  startup,  th'  gradient 
Is  very  steep  indicative  of  strong  vapor  diffusion  activity  and  at 
steady  state,  it  Is  spread  out  due  to  axial  conduction  effects. 
Hence,  it  appears  that  the  transient  state  is  governed  by  the 
diffusion  process  and  the  steady  state  by  axial  conduction.  This 
hypothesis  is  verified  through  numerical  and  experimental  results 
(see  results  In  Section  5.2.2). 

Ihe  radiation  exchange  (Q^)  between  the  hot  and  cold  zones 

is  less  significant  than  Q  or  Q  since  the  length  to  diameter 

v  ax 

ratio  of  the  pipe  Is  very  large. 
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In  the  present  work,  a  two-dimensional  analysis  of  the  vapor 
diffusion  through  the  NCG  is  performed.  The  diffusive  energy 
transport  rate  thus  obtained  is  coupled  with  a  thermal  model  of  the 
startup  of  the  radiatively  cooled  sodium  argon  double  wall  artery 
heat  pipe. 

2.4  Diffusion  Model  and  Analysis 

The  objective  of  the  diffusion  analysis  Is  to  determine  the 
vapor  diffusion  rate  past  the  hot  front  as  a  function  of  the  length 
and  temperature  of  the  hot  zone.  The  model  is  based  on  tick's  law 
for  binary  gas  system  explained  in  detail  by  Bird  et  al.  [51]  and 
developed  by  Sockol  [17].  Assumptions  used  in  this  analysis  are  as 
follows : 

(1)  The  local  condensation  rate  upstream  of  the  front  Is  assumed  to 
be  negligible. 

(2)  The  cold  wall  acts  as  a  perfect  sink  for  the  vapor. 

(3)  The  transition  region  between  the  cold  and  hot  zones  is 
vanishingly  thin. 

(4)  The  upstr^m  vapor  velocity  Is  uniform  over  the  cross  section. 

(5)  The  gas  Is  stagnant  throughout  the  pipe. 

(6)  The  axial  vapor  pressure  drop  Is  negligible. 

(7)  The  vapor  velocity  during  startup  is  limited  to  much  less  than 
sonic  speed. 

(8)  In  a  long  cylinderlcal  heat  pipe,  the  vapor  concentration 

should  relax  to  a  quasi  steady  profile  in  time  of  order 
2 

(rv/2. 4)  /0.  If  the  distance  moved  by  thp  front  during 
this  time  is  small,  steady-flow  equations  can  be  used  for  the 
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diffusion  analysis.  For  the  sodium  argon  system  at  1000  K, 

2 

with  rv  -  0.635  cm  and  D  -  4.591  cm  /sec,  the  relaxation 
time  Is  15.25  milliseconds.  As  it  will  be  seen  in  the 
numerical  analysis  results  section  later,  the  rate  of  front 
movement  is  0.04  to  0.37  cm/sec  for  a  power  input  of  1000  W 
with  the  gas  charge  pressure  of  0.2  torr.  Hence,  the  front 
movement  in  15.25  milliseconds  Is  very  small  and  the  steady 
state  diffusion  assumption  is  justified. 

Pick's  law  for  binary  diffusion  In  the  form  of  molar  flux  is  given 
by  [51] 


N  =  x  (N  f  N  )  -  cD  Vx 
v  v'  v  g'  vg  v 


(1) 


lor  a  stagnant  gas  Ng  *  0  and  Eq.  (1)  becomes, 


N  -  x  N  -  cO  vx 
v  v  v  vg  v 


or 


<'  ’  \>"v  '  -c0,g'\ 


Total  mass  fraction  of  the  system,  x  +  x  *  1 

g  v 


Therefore,  the  vapor  flux  i$v  can  be  written  In  terms  of  the  mole 
fraction  of  the  gas  x^  as: 
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cDya'(1 

xg 


x  ) 
.  g 


or 


vg 


which  also  can  be  written  as 


Nv  •  cDvg7  an  xg 


(2) 


From  the  principle  of  conservation  of  mass  applicable  to  the  vapor 
flow. 


dlv  Ny  -  0 

or  V  •  Nv  =•  0  (3) 

Substituting  for  Ny  from  Fq.  (2), 

cIV2  *n  xg  -  0 

Since  cDvg  /  0  ,  7^  an  xg  -  0  (4) 

In  order  to  solve  this  equation  and  to  set  up  the  necessary  boundary 
conditions,  a  cylindrical  coordinate  system  is  Introduced  as  shown 
in  Figure  7.  According  to  Assumption  (3),  the  axial  distance  over 
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Figure  7.  Coordinate  System  for  Diffusion  Analysis. 


which  departs  from  0  or  1  is  small  compared  to  the  hot  zone  or 
cold  zone  lengths,  from  the  location  of  the  front,  the  boundaries 
of  z  are  at  ■»  <®.  The  boundary  conditions  are  as  follows: 
z  <  0: 


(1)  As  z  »  <»,  the  vapor  flux  is  given  by  an  unknown 

constant  flux  N  . 

o 


-» 

N  -> 
v 


-» 


That  is  f rom  \ q .  (?) , 


cD  *  (In  x0)  »  N0  as  z  ->  ®  (5) 

dz  9 

(2)  At  the  wall,  the  vapor  flux  (Ny)  in  the  radial  direction  is  zero 
as  there  is  no  condensation  in  the  hot  zone. 

3  (in  xq)  --  0  at  r  =  rv  (6) 

3r  9 
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(7) 


z  _>  0: 


(1)  At  the  condenser  end,  there  is  no  vapor;  only  gas  exists. 


x  -*  1  as  z  -»  t  ® 
9 


(2)  At  the  wall  (perfect  sink),  there  is  no  sodium  vapor  due  to 

complete  condensation,  x  »  0  and  x  -  1 

v  g 


x  =  1  at  r  =  r  (8) 

g  v 


A  mass  balance  on  the  gas  is  done  to  obtain  the  unknown  vapor 

flux  Nq  in  terms  of  the  initial  molar  concentration  c^  of  the  NCG  at 

T  . 
c 


(mass  of 
gas 

in  pipe 


(mass  of 
gas 

in  hot  zone 


+ 


(mass  of 
gas 

in  cold  zone 


mass  of 
gas 
in  hot  zone 


(mass  of 
vapor  +  gas 
in  cold  zone 


(mass  of 
vapor 
in  cold  zone 

(9) 


If  the  total  length  of  the  heat  pipe  is  Lp  and  the  hot  zone  length 
is  L,  then  Eq,  (9)  can  be  expressed  in  the  Integral  form  as. 


-l 


dz 


* 


a* 


i 

—a> 


cxg  dz  + 


0 


dz 
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where  c  Is  the  total  molar  density  of  the  mixture  and  the  bar 
indicates  an  average  over  the  cross  section.  Simplifying  the 
integrals  and  bringing  the  constant  c  (as  cc  and  c^  evaluated  at 
and  T,  respectively)  outside  the  Integral,  one  obtains 


0  00 

c1Lp  *  ch  j  *g  1)2  +  cc(Lp  -  L)  -  J  <’  -  V  d2  (10) 

-<D  0 

From  ideal  gas  law,  the  molar  densities  at  various  temperatures  are 
expressed  as. 


P  (T) 
v 

RT 


Ch  ' 


P  (1) 
v 

RT 


where  -  Initial  charge  pressure  of  NCG. 


Introduce  dimensionless  variables  as  follows: 


N  r 

e  =  I_  ;  g  -  -JL-V  (Diffusion  Parameter) 
T  cD 

c 
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ni) 


c  (L  -  L)  -  c.L  l  P. 

\(6,t)  -  _?_JL  .  --  I-B  -  lD  -  l  .P .JL 

c  r  H  P  (1) 

C  V  V 


A  new  dependent  variable  y  Is  Introduced  In  place  of  in  x^  as. 


y(r\z) 


m  x. 


(12) 


Eq.  (4)  in  cylindrical  coordinates  transforms  to 


1  a_  p  3y  +  i_y  *  0 
r  ar  3r  a£2 


(13) 


Boundary  conditions,  Eq.  (5)  to  (8)  become, 

z<0:  3Y  •>  -  l  as  z  -  ® 

3z 

3y  =  o  at  r  =  1 

ar 

z>0:  y  -»  0  as  z  ® 


(14) 


y  *=  0 


at  r  =  1 


Rearranging  and  normalizing  Eq.  (10), 


x(e,l) 


1 


(1  -  Xg)  dz 


--  I 

cc  3 


Xg  dz 
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Averaging  the  gas  concentration  over  the  cross  section. 


x(e,t) 


2 


1 


0 


,-gy 


)di 


Letting 


fc(g>  -  2  J 

0 

1 

F„(g)  *  2  J 

0 


r  dr 


r  dr 


j  (1-e‘9y)dz 
0 


0 


Eq.  (15)  becomes, 


X(0,t)  -  Fc(g)  -  1  Fh(g) 

e 


05) 


06) 


07) 


08) 


So lotion  Method 

The  partial  differential  F.q.  (13)  and  the  boundary  conditions 
fq.  (14)  are  solved  for  y(r,z)  In  the  closed  form  using  Bessel 
functions  as  described  by  Sockol  [17].  Once  y  is  known,  g  can  be 
obtained  by  solving  Fq.  (18).  The  mass  flux  In  terms  of  0  and  t 
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are  obtained  by  performing  the  radial  integral  numerically.  A 

simplifying  approximation  is  used  for  this  purpose.  For  large 

values  of  g,  F.(g)  is  negligible  compared  to  F  (g)  and  the 
h  c 

radial  Integral  yields. 


Fc(9)  ^  1  Un  g  -  1.339)  (19) 

Pi 

where  p^  -  2.4048  is  computed  from  the  zeros  of  Bessel  function. 
Now,  the  approximate  equation  for  Fc ( g )  is 


Or 


Fc(g)  •*  \(e,i)  *  I  (tn  g  -  1.339) 
0 


g  =  exp(l  .339  +  p^) 


(20) 


(21) 


The  energy  transport  by  vapor  diffusion  from  the  hot  to  cold 
zone  Qv  is  given  by 


0v  *  AyMyNghfg 


7 

[cm  x 


g _ x  gmol  e 

gmole  2  , 

”  cm  sec 


-J] 

g 


=  *rv"Ahfg  W 

Or  in  terms  of  cO  and  g. 


Q 


v 


=  wr  M  cOh 
v  v 


(?2) 
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The  values  of  cD  for  sodium  argon  mixture  as  a  function  of 
temperature  are  computed  from  the  kinetic  theory  of  gases  using  the 
Chapman  Fnskog  relation  as  outlined  In  Appendix  A.  The  diffusion 
parameter  g  has  to  be  very  large  at  the  beginning  of  the  startup  In 
order  to  have  large  and  must  drop  to  zero  toward  the  end  of 

startup. 

2.5  T ransient  Thermal  Analysis 

The  transient  state  of  the  gas  loaded  liquid  metal  heat  pipe 

with  a  long  adiabatic  section  at  time  t  after  the  application  of 
evaporator  power  input  Is  shown  in  Figure  8.  The  top  sketch 

Figure  8(a)  shows  the  heat  pipe  during  the  thawing  process  when 
there  is  no  heat  output  at  the  condenser.  Figure  8(b)  shows  the 

condition  when  the  pipe  Is  transporting  and  radiating  out  Qq  at 

the  condenser. 

2.5.1  Assumptions 

(1)  Axial  conduction  Is  negligible  compared  to  Qy  during  startup. 

(2)  The  evaporator  and  adiabatic  sections  are  perfectly  insulated 

(shielded)  such  that  no  heat  loss  occurs  in  the  radial 
direction. 

(3)  Axial  and  radial  temperature  variations  within  the  same  zone 
(hot  or  cold)  are  not  considered.  The  hot  zone  is  at  a 
uniform  temperature  T  and  the  cold  zone  (gas  slug  region)  is 
uniformly  at  1  .  Only  the  time  dependence  of  temperature  and 
the  location  of  the  hot  front  are  analyzed,  i.e.,  T  -  T(t)  ; 
L  -  L(t). 
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(4)  There  Is  no  liquid  return  flow  to  the  evaporator  during  the 
thawing  stage  and  there  Is  enough  fluid  Inventory  to  last  till 
the  pipe  starts  up. 

(5)  Initially,  the  wick  structure  is  saturated  with  the  working 
fluid  which  Is  in  a  frozen  state. 

(6)  The  heat  pipe  vapor  core  is  filled  with  a  NCG  (argon)  charge 
at  a  pressure  and  temperature  1  .  There  is  no  separate 
gas  reservoir  In  the  pipe. 

(7)  No  chemical  reaction  occurs  between  the  fluid  and  gas. 

(8)  Ihe  mechanism  of  vapor  diffusion  Into  the  gas  within  the  vapor 
core  is  as  outlined  In  the  previous  section. 

(9)  Horizontal  orientation  of  the  pipe  is  assumed  during  startup. 

(10)  Heat  input  to  the  evaporator  is  uniform  along  its  length  and 

it  can  be  conductively  or  radlatively  coupled  whereas  the  heat 

removal  at  the  condenser  Is  radlatively  coupled  to  the  cold 

sink  at  T  . 

s 

(11)  lhe  reservoir  wick  at  the  end  of  the  evaporator  is  in  good 
capillary  contact  with  the  evaporator  wick. 

2.5.2  Definitions 

( 1 )  In U la.  1  hot  zone  temperature,  p r  rise  tempe ratu re  ( 1  . )  is 

HI 

the  temperature  at  which  the  sodium  vapor  pressure  Is  equal  to  the 

adjusted  Initial  NCG  charge  pressure  such  that 

P v ( 1 H 1 )  ^  Pi  <j- ---»  -) 

p  '  4 
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(2)  Startup  temperature  (THs)  is  the  temperature  of  the  hot  zone 
when  the  hot  front  has  reached  the  beginning  of  the  condenser. 

(3)  Rise  time  (t,.)  Is  the  time  requ'red  to  raise  the  temperature 

of  the  evaporator  to  1  ,  from  1  (room  temperature)  after  the 

HI  C 

Qi  is  applied. 

(4)  Startup  time  (t§)  is  the  time  taken  after  t  for  the  hot 

front  to  move  to  the  beginning  of  the  condenser,  i.e.,  to  move 

through  the  adiabatic  length,  L. . 

A 

2.5.3  Governing  Differential  Equations 

Energy  balances  are  conducted  for  the  control  volumes 
enclosing  the  hot  and  cold  zones  of  the  heat  pipe  during  the 

transient  states  of  hot  zone  diffusing  Into  the  cold  zone  as 
depicted  in  figure  8.  Two  sets  of  equations  corresponding  to  the 

thawing  stage  and  transporting  stage  for  each  of  the  control  volumes 
are  written  with  the  initial  and  boundary  conditions.  The  symbols 
are  explained  in  the  nomenclature. 

Hot  Zone: 

Thawing  stage:  0  <  t  <  t$  ;  i  L  s  (Lf  f  L^) 

[CElE  +  CA(l  -  lE) ]  41  ,  Qi  -  Qv  -  Or  (?3) 

dt 

with  boundary  condition  :  T  =•  Tc  at  t  -  0 
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Transporting  Stage:  t$  <  t  <  °®  ;  (L£  ^  Lft)  S  L  i  Lp 


[CELE  f  C^La  *-  Cc( L  '  tE  -  La>}  ~  *r  °1  "  °v  Or  '  Qo  <24> 


with  boundary  condition  :  L  -  (L^.  «■  l-A>  at  t  --  t 


Cold  Zone: 


Thawing  Stage:  0  <  t  <  t  ;  L.  z  L  £  f  LA^ 


Ca(1  -  Tc)  jjk  .  Qv 

with  boundary  condition  :  L  :  I t 


at  t  =  0 


Transporting  Stage:  te  <  t  <  ®  ;  (Lt  f  1  6  $  L 


CC(T  -  1C)  £  '  Ov 
at 

with  boundary  condition  ;  l  =  (L^  4-  L^)  at  t 
Other  known  relations  are, 


Q  --  «r  M  cDh  g 
v  v  v  fg 


(Diffusion 

Equation) 


Q  ■  irD  (I.  -  Lr  -  U)  cnu(T4  -  T4)  (External 
00  E  A  p  Radiation) 


0r  “  2irrvL^HCcHc,^T  ”  V 


(Internal 

Radiation) 


Eqs.  (23)  to  (29)  are  normalized  using  the  following  nondimensional 


variables . 


6  -  J  ;  e$  .  s 


45 


L  L  I--  l 

1  -  L  ,  lE  =  .1.  .  lA  =  A  ,  4C  =■  A  .  ftp  '  -& 

rv  rv  rv  rv  rv 


Qi  = 


CE  rv  7c 


<lv  " 


Q  Qr 

_v —  ,  qr  *• 

Cr  r  T  C.  r  T 

Eve  t  v  c 


,0  -  ^  ■  - 


s2 


The  normalized  equations  and  the  boundary  conditions  are. 


[If  f  st(4  -  i£)]  «  <11  -  dv  '  Qr 

dt 


(23) 


with  boundary  condition  :  0-1  at  t  =  0 


or 


0  <  t  <  t,  ;  lt  U  <  (lE  ^  lft) 


[lE  f  s-j 8-a  +  s2(l  "  lE  ~  1a)  1  ^  K  Ql  qv  *  ^r  % 


(24) 1 


with  boundary  condition  :  l  -  (4^.  +  4^)  t  -  t 


for  t  <  t  <  ®  ;  (4f  Ma)  U  Up 


si (e  -  1 )  ^  -  dv 
1  dt 


(25) 1 


boundary  condition  :  4  -  4^ 


at  t  «  0 


s2(e  -  1)  -  qv 

£  dt 


(26)' 


boundary  condition  :  4  =■  (4^  f  4ft)  at  t  =■  0 
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(27)' 


q  *-  A  g 
v  v 


where 


*H  cDh, 
Av  = ¥ — -f-fl 

Vc 


%  '  Ao(l  '  lA  -  *E>  (e4  -  6s> 


(28)* 


where 


tOcaT. 

A0  =  --SJL-C 


qr  =  Arl(e  -  1) 


(29) 1 


where 


Ar  =  2*FHCcHqTc 
CE 


from  the  diffusion  analysis  we  have, 


g  =  exp(l  .339  «■  2.4048  \(6,l)] 


(21) 


I  P. 

X(0,l)  -  In  “  l  -  -£-± 


PV(T) 


(ID 


With  this,  the  formulation  of  the  transient  problem  Is  completed  and 
the  numerical  solution  technique  is  described  in  Section  2.6. 

2.6.  N ume r leal  Computation s 

The  governing  differential  equations  and  the  Initial 
conditions,  F.qs.  (23)'  to  (29)',  are  solved  by  numerical  Integration 
using  the  fourth  order  Runge  -Kulta  method  for  the  2  m  sodium  argon 
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heat  pipe  design  described  in  the  next  chapter  [52].  A  FORTRAN 
program  called  SOOART  was  developed  for  this  purpose.  SOOAR1  works 
In  conjunction  with  several  subroutines  written  for  specific 
purposes.  Table  1  lists  all  those  subroutines  and  their  functions. 
A  flow  chart  illustrating  the  sequence  and  logic  is  given  In 
Figure  9.  A  listing  of  the  SOOART  program  and  subroutines  are  given 
in  Appendix  B. 

2.6.1  Input  Data 

The  input  data  required  for  the  program  is  given  in  Table  2. 

The  temperature  dependent  properties  such  as  the  diffusion 

coefficient  (cD)  and  the  enthalpy  of  vaporization  (h.  )  are  input 

fg 

as  coarse  table  and  Interpolated  within  the  program  for  fine 

intervals.  Since  P  (T)  is  an  exponentially  varying  property,  the 
empirical  relation  [Eq.  (30)]  by  Dltchburn  and  Gilmour  [53]  is 
used . 

Py  --  2.29  x  1011  (T~05)  (10-5%//,)  (30) 

2 

where  T  is  in  degrees  K  and  Py  in  N/m  .  A  sample  input  data  file 

is  also  given  along  with  the  program  In  Appendix  B. 

2.6.2  Computation  of  Startup  Parameters 

The  startup  parameters  to  be  computed  are  the  rise  time  and 

temperature  (t^,  T^),  startup  time  and  temperature  (t^,  T  ) 

and  the  sodium  mass  depletion  rate  from  the  evaporator  (m).  The  NCG 

Initial  charge  pressure  [P.(l  )]  and  the  evaporator  power  input 

1c 

(0.)  are  the  main  Input  parameters  to  start  the  computations.  For 
the  given  ,  the  startup  vapor  pressure  is  found  from: 
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TABLE  1.  LIST  OF  SUBROUTINES  ANO  THEIR  PURPOSES 


NO.  SUBROUTINE  FUNCTION 


1.  INPSUB  Reads  the  input  data  and  transfers  to  the  main 

program.  Calls  WRITEUP  to  printout  the  input 
data . 

2.  RUNGE  Performs  the  numerical  integration  of  the 

differential  equations  (first  order  time 
derivatives  of  length  and  temperature  dL/dt  and 
dT/dt)  with  the  precision  of  time  step  (TSTEP) 
specified.  Runge-Kutta  method  is  used. 

3.  CALC  Computes  g,  Qv,  Qr.  and  Q0  from  the 

governing  differential  equation.  Uses  INTPOLT 
and  PVCALC  to  compute  the  temperature  dependent 
properties,  hfg,  cD  and  Pv. 

4.  OERV  Converts  the  differential  equations  to 

differential  quantities  for  numerical 
Integration. 

5.  INTPOLT  Interpolates  property  data  for  any  temperature 

during  the  integration. 

6.  WRITEUP  Prints  the  input  and  output  data. 

7.  PVCALC  Uses  an  empirical  relation  to  calculate  Pv  of 

sodium  as  a  function  of  temperature. 

8.  STARTUP  Calculates  the  starting  temperature  for  the 

given  gas  charge  pressure  (P^). 


49 


t-TIME  ♦  TSTEP 


Call  CALC 
updates  G.Qy.Qr, 
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Figure  9.  Flow  Chart  for  the  Computer  Program  SODAR? 


TABLE  2.  INPUT  DATA  FOR  PROGRAM  SODART  FOR  2  m  1 1  QUID  METAL  DWAHP 


1 . 

Vapor  Core  radius  (rv) 

=  0.635xl02  m 

2. 

Evaporator  length  (L£) 

-  0.375  m 

3. 

Adiabatic  length  (lA) 

-  0.745  m 

4. 

Condenser  length  ( L^ ) 

=  0.910  m 

5. 

Total  length  (Lp) 

-  2.030  m 

6. 

Evaporator  unit  thermal 
capacity  (Ct) 

=  911.717  J/m  K 

7. 

Adiabatic  unit  thermal 
capacity  ( CA) 

--  598.934  J/m  K 

8. 

Condenser  unit  thermal 
capacity  (Cc) 

-  886.340  J/m  K 

9. 

Molecular  weight  of 
sodium  (Mv) 

=  23 

10. 

Emissivity  of  condenser 
surface  (ea) 

*  0.66 

11  . 

Stefan  Boltzmann  constant  (0) 

=  5.669x10-8  w/m? 

12. 

Emissivity  of  vapor  core 

T  3 . 

14. 

15. 

16. 

17. 

18. 


surface  (eH) 

Shape  factor  between  hot 
arid  cold  zone  ( Ff1c ) 

Room  temperature  (lt) 

Chamber  sink  temperature  (Ts) 

NCG  charge  pressure  (Pj): 

fvaporator  input  ( Qi ) : 

Heat  input  increments 
( 1THLA1) 


19.  Solution  stop  time  ( I  FINAL) 


;  0.1 

-  0.01159 

-  300  K 

-  300  K 

0.05,  0.10,  0.20,  0.50, 
1.00,  2.00  torr 

100,  250,  500,  1000  W 


3000  sec 
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(31) 


Pv<THi)  =  PiHc)  £ 


Up  - 1() 


From  the  sodium  vapor  pressure  relation  CM-  (30)],  f  corre 

H  l 

sponding  to  P  (f  )  Is  computed.  Applying  the  energy 

V  H 1 

conservation  principle  to  the  hot  zone  and  neglecting  the  losses, 
one  obtains  the  equations  for  t ^  and  t  . 


ti 


(32) 


*s 


_ *1 _ 

(Cf4  f  CALA><THs 


(33) 


T  is  computed  from  the  SOOAR1  program  by  solving  the  governing 
Hs 

differential  equations  (Section  2.5.3) 

During  the  time  of  thawing,  vapor  mass  flow  rate  is  governed 
by  Qv.  If  m  is  the  rate  of  mass  depletion  in  the  evaporator  during 
startup,  Qv  can  be  written  as, 


-  irr  M  N  hf 
v  v  o  fg 


-  m  h 


fg 


or 


m  - 


(33) 


By  integrating  m  with  respect  to  time,  the  total  mass  of  sodium 
vaporized  f rom  the  evaporator  and  transferred  to  the  adiabatic 
section  is  obtained. 
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m  dt 


(34) 


t 

M-IS 

0 

This  integration  is  automat ica Ily  available  from  the  computa 

lion  in  SODARI. 

2.6.3  Res ul ts  of  Compu tatlons 

A  matrix  of  computations  was  done  varying  the  P.  and  Q. 
values.  For  every  value  of  (i.e.,  100,  250,  500,  and  1000  W) , 
six  P^  values  (0.05,  0.1,  0.2,  0.5,  1.0,  2.0  torr)  were  varied, 
and  a  total  of  24  computations  was  done.  In  each  set  of 

computation,  the  transient  results  t  versus  I.,  I,  Qv>  Qq,  Qr, 
and  M  were  computed  for  3000  sec  and  graphed.  The  time  step  used 
was  10  sec  and  the  CPU  time  consumed  in  a  COC  mainframe  computer 
was  220  sec. 

(able  3  lists  the  results  of  the  startup  parameters.  It  is 
observed  that  the  parameters  time,  temperature  and  mass  increase 

with  increase  in  NCG  charge  pressure,  P^ ,  while  the  location  of 
the  hot  front  decreases.  Figure  10  shows  the  variation  of  the  hot 
front  length  with  power  input  for  various  NCG  charge  pressures  while 
Figure  11  shows  the  same  result  in  different  form.  The  sodium  mass 
depletion  from  evaporator  during  startup  as  a  function  of  NCG  charge 
pressure  is  plotted  in  Figure  12.  The  important  design  input 
variable  for  heat  pipe  operating  condition  in  terms  of  selecting  the 
NCG  charge  pressure  is  obtained  from  these  results,  i.e.,  the 
integrated  sodium  mass  depletion  values.  Only  the  mass  depletion 
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TABLE  3.  RESULTS  OF  THE  STARTUP  ANALYSIS 


tion  of  hot  front  100  1.317  1.291  1.267  1.238  1.218  1.200 

000  sec  after  250  1.620  1.564  1.511  1.446  1.401  1.360 

time,  L  (meter)  500  1.901  1.846  1.784  1.695  1.626  1.560 

1000  2.011  1.999  1.980  1.935  1.886  1.882 


tiu 


200.0 


400.0  600.0 

Q  (WATTS) 


800  0 


1000.0 


Figure  10.  Heat  Front  Length  Versus  Heal  Input 


55 


IOOW 


OS  10  IS 
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Figure  11.  Heat  Front  length  Versus  NCG  Pressure 


NA  MASS  DEPLETION  (G) 


value  of  38.0  g  for  0.?  lorr  satisfied  the  design  conditions 

(evaporator  jrid  reservoir  wick  hold  38. 23  g);  hence,  this  value  of 
the  charge  pressure  was  recommended.  Other  charge  pressures 
(0.05  2.0  Torr)  also  could  be  tried  during  the  experimental  phase. 
The  startup  temperature  necessary  for  a  successful  start  with 
reference  to  the  characteristic s  of  liquid  metal  heat  pipe  startup 
(Section  2.1)  was  obtained. 

The  major  aspects  of  the  transient  results  of  the  computation 
are  presented  in  the  graphs  of  Figures  13  through  28.  Fach  graph 
shows  six  curves  corresponding  to  the  six  charge  pressures  and 

represents  one  of  th'  'our  parameters  (L,  1,  Q^,  or  Q^). 

Figures  13  through  16  correspond  to  -  100  W  and  show  L  vs  t, 
T  vs  t,  Qy  vs  t,  and  0Q  vs  t,  respectively.  Ihree  other  sets  of 
figures  (Figures  17  through  28)  show  the  results  for  250,  500  and 
1000  W  in  order. 

lhe  general  trend  of  the  Qv  vs  t  curves  is  importaril.  At 

the  starting  time,  0^  is  as  high  as  Q.  and  as  the  hot  front 

moves,  Qv  drops  exponentially  to  very  low  value  at  the  end.  lhe 

kink  in  the  curve  where  the  slope  changes  signifies  the  abrupt 

change  in  the  thermal  capacity  of  lhe  adiabatic  to  condenser 

sections.  On  the  other  hand,  Q  remains  zero  until  the  hot  front 

o 

reaches  the  condenser  and  increases  fast  to  a  maximum  toward  the 
end.  lhe  variation  of  hot  zone  temperature  during  startup  is  only 
marginal.  The  rise  temperature  is  as  high  as  the  startup 
temperature  and  this  signifies  the  frontal  nature  of  the  startup. 

These  results  are  compared  and  corroborated  with  the 

experimental  data  in  the  results  and  discussion  section. 
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Figure  13.  Heat  Front  Length  Versus  Time  for  100  W 
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Figure  16. 


Heat  Radiated  Out  Versus  lime  tor  100  W 
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Figure  17.  L  Versus  t  for  250  W 
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f- iyure  19.  Qy  Versus  t  for  260  W 
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figure  20.  Q0  Versus  t  for  250  W 
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Figure  21.  Qv  Versus  t  for  1000  W 
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figure  28. 
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OF  S  l  GN  ANALYSIS  OL  I  HE  LXPLRlMtNlAI  IttAl  PII'L 
3.1.  Design  Cons traints/Requi remen ts 

The  starting  point  of  any  heat  pipe  design  is  to  consider  the 
operating  temperature  regime  of  interest,  geometries  and  the  heat 
transport  capacity  requirements.  The  objective  of  the  present  work 
is  to  investigate  the  startup  characteristics  of  the  high 
temperature  DWAHP  in  the  1000  K  regime  with  a  selected  nominal 
laboratory  testable  load  of  2000  W.  1  he  most  appropriate  working 
fluid  in  the  temperature  range  400  1000"C  is  sodium  and  the 
compatible  container  easily  available  is  stainless  steel.  An 
overall  length  of  ?  m  and  a  tube  diameter  of  ?.??  cm  was  chosen  for 
the  physical  dimensions  for  the  heat  pipe.  I  he  geometry  of 

evaporator,  adiabatic  (transport)  and  condenser  sections  of  the  heat 
pipe  is  always  dictated  hy  the  application  requirements  of  t  tie 
device.  However,  in  a  research  device,  the  sizing  of  various 
sections  is  flexible.  the  evaporator  length  selected  by  the  input 
heat  t  lux  limit,  I ixed  the  type  of  heating  employed.  Radial  heat 
fluxes  up  to  100  W/cm  are  common  for  sodium  he.it  pipes.  Ihe 

condenser  length  is  dictated  by  the  type  of  heat  transfer  coupling 
available  between  the  heat  pipe  a'-J  (lie  sink.  Radiative, 
convective,  or  conductive  coupling  is  possible.  In  order  to 
simulate  space  radiator  conditions  and  for  simplicity  of  laboratory 


testing,  radiative  coupling  with  a  cold  shrouded  vacuum  chamber  was 
chosen . 

Since  the  present  study  is  an  offshoot  of  the  ?.  m  copper  water 
DWAHP,  the  mechanical  design  details  such  as  1  he  type  of  screen 
wick,  artery  groove  sine,  double  wall  feature,  etc.,  were  kept  the 
same  as  previous  conf  igur.rt  ions  for  comparison  purposes  [10,11  |. 
Some  of  the  special  features  included  in  the  experimental  heat  pipe 
design  are: 

(1)  Single  groove  artery  transport  section  to  study  the 
effects  of  orientation  of  the  groove  with  respect  lo  the 
gravity  vector  iri  I  g  testing; 

(?)  Nonroridensible  gas  ( NCfj)  loading  to  study  the  rapid 
startup  from  frozen  slate. 

from  a  fabrication  point  of  view,  the  geometrical  shape  and  size  of 
t  tic  groove,  vent  slots,  end  cap  design,  screen  wick  assembly, 
cleaning  and  filling  procedures  were  chosen  to  he-  simple,  safe  and 
within  the  limits  of  available  resources. 

The  basic  design  parameters  available  to  start  a  detailed 
design  are  summarized  as  follows: 

•  Working  fluid  Sodium 

•  Noncondensible  gas  Argon 

•  Temperature  range  ?5  -1000°C 

•  Heat  load  ?000  W 

•  Container  material  Stainless  steel 

•  Wick  material  SI  aidless  steel  screen 

•  lota  I  length  ?  m 
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•  Outer  diameter  ? .??  x  10  m 

•  Meat  input  mode  Direct  contact  electrical 

res  i  s  I  uru  e  coil  heat  ing 

•  Heat  removal  mode  Radiation  coupling  to  a  cold 

stiroud 

3.2.  Material  Selection  and  Compatibility 

Sodium :  Sodium  is  a  well  charact eri zed  substance  and  is 

available  in  very  pure  analytical  grades.  It  has  high  thermal 

conductivity  in  solid  and  liquid  phases.  the  low  vapor  pressure 

5  ?. 

even  at  high  temperature  (0.2x10  N/nr  at  1000  K)  makes  sodium 
pressure  containment  easy  since  a  relat  ively  thin  walled  container 
is  adequate.  Ihe  liquid  transport  factor  or  figure  of  merit 
(cy»  X/p  )  for  sodium  is  3X10”1  kW/crr/  at  1000  K.  It  is  compatible 
with  stainless  steel,  inconel,  nickel  and  niobium  which  are 
candidate  wall  materials  for  high  temperature  heat  pipes.  Ihe 

intrinsic  drawbacks  with  sodium  are  its  handling  hazards  and  the 
frozen  slate  at  room  temperature.  txtreme  care  is  required  in 
handling  sodium  in  order  to  ensure  personnel  safety  and  heat  pipe 
working  fluid  purity.  A  list  of  selected  properties  is  given  in 

Appendix  C. 

Stainless  Steel:  Ihe  steel  alloy  chosen  here  is  SS  304  which 
has  t tie  property  of  good  corrosion  resistance.  While  inconel  may 

have  been  preferred  over  stainless  steel  from  high  temperature 
strength  and  corrosion  resistance  points  of  view,  availability  in 
required  size  tubings,  bar  slock  and  woven  screen  mesh  dictated  the 
selection.  Ihe  ultimate  tensile  strength  of  SS  304  drops  from 

(.20x10  N/m  (MO  kl’Sl)  at  300  K.  to  207x10  ‘  N/m'  (30  kl'Sl)  at  1100  K . 


II 


The  end  cap  s  i  ze  and  wail  thickness  of  the  heat  pipe  envelope  are 
governed  by  this  high  temperature  strength  limitations. 


High  Temperature  (ill  Valve:  A  valve  at  the  condenser  end  ol 
the  research  heat  pipe  is  essential  for  processing  such  as 
evacuating,  filling  the  working  fluid  or  the  noncondensible  gas,  and 
burping  the  excess  I  luid/gus  if  necessary.  A  high  temperature, 
bellows  type,  all-metal  (stainless  steel),  all  welded  construction 
valve  manufactured  by  NUPRO  was  selected.  This  valve  has  a  stellite 
tip  for  the  valve  stem.  It  is  rated  for  service  at  6/19°C  and 
17,200  kPa  or  at  extended  ratings  of  8I6°C  and  340  kPa  if  self 
lubricating  fluids  such  as  liquid  metals  are  used.  A  secondary 
packing  of  grafoil,  lubricated  by  mo’y  disulfide  base  paste,  seals 
the  valve  for  protection  of  personnel  in  case  of  hazardous  fluid 
leaks . 

Argon:  fhe  noncondensible  gas  (NCG)  chosen  is  the  99.99% 

pure,  dry,  oil  and  moisture  free,  laboratory  grade  5.8  AIHCO 

supplied  specialty  gas,  argon.  The  requirements  of  the  NCG  in  the 

present  application  are  low  solubility  and  non  reactivity  with 

sodium  even  at  high  temperatures.  The  solubility  constant  (Ostwald 

-5 

coefficient)  for  argon  and  sodium  at  500 °C  is  ?.?/xlO  and  it  is 
considered  that  argon  is  insoluble  in  sodium  for  practical  purposes. 
Also,  argon  is  inert  and  does  not  form  any  compound  with  sodium.  In 
contrast,  nitrogen  could  be  used  as  the  NCG  since  the  solubility 

constant  is  very  low  (3.11x10  at  500°C);  however,  it  is  not 

considered  because  sodium  is  known  to  react  with  nitrogen  at  high 
temperatures  (>300°C)  forming  sodium  nitride. 
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3.3  Mechanical  Design  Details 

A  detailed  design  is  carried  out  based  on  the  available  design 
parameters  and  constraints  in  order  to  determine  the  following: 

(I)  Vapor  core  diameter  to  provide  a  vapor  Marti  number  of  0.? 
which  would  eliminate  vapor  compressibility  effects  and 
sonic  flow  limitations. 

(?)  The  container  details  such  as  evaporator,  adiabatic  and 
condenser  section  lengths,  wall  thickness,  and  end  cap 
d imens ions . 

(3)  Ihe  capillary  wick  details  and  flow  pressure  drops. 

(4)  Ihe  artery  groove  si/e  (depth  and  width)  and  screen  wrap 
detai  Is . 

(b)  The  heat  transport  limits  such  as  capillary,  sonic, 

boiling,  and  condenser  radiation  limitations. 

(6)  The  fluid  inventory  and  NCG  initial  charge  pressure. 

3.3.1  F  n velope  Design 

The  chosen  total  length  of  2  m  is  to  be  divided  into 

evaporator,  transport  and  condenser  lengths.  from  the  radiation 
limit  calculations  (refer  to  Section  3.S.1),  it  is  found  necessary 
to  have  a  minimum  of  0.91  m  long  condenser  section.  The  heat  pipe 
will  be  radiation  limited  if  the  condenser  is  shorter  than  0.91  m. 
The  lumped  thermal  capacity  analysis  (Section  3.4)  is  used  to 

determine  the  length  of  evaporator  needed  mainly  from  tfie  standpoint 
of  working  fluid  requirement  during  startup.  Accordingly,  the 

evaporator  length  is  0.37b  m  which  is  divided  into  a  reservoir  wick 
length  of  O.I?b  m  and  the  heated  evaporator  length  of  0.?b  m.  The 


remaining  heat  pipe  length  of  0.  /4b  m  is  available  lor  the  Iran-. purl 
se<tion  which  is  sufficiently  long  for  the  experimental  laboratory 


study.  lo  summarize  this  design,  the  total  length  (I  )  and 
effective  length  (I  ^)  are  written  as  follows: 


L 


P 


(3*>) 


or  2.03  m  -  0.3/5  m  i  0.  /45  m  i  0.91  m 


L  rr  0.6  I,  ♦  L  e  0.5  I 
eff  t  A  C 


(3b) 


=  1 .388  m 


The  minimum  vapor  core  diameter  required  for  an  axial  heat  flux  of 
2000  W  at  1000  K  nominal  operating  temperature  and  vapor  Mach  number 
limit  of  0.2  is  obtained  as  d^  8.3?  mm  from  f.q.  (37). 


dv  = 


20  q 


max 


i rp  V  J  y  R  T 
v  Vvvv 


(3  7) 


The  vapor  pressure  of  sodium  at  a  possible  maximum  operating 

temperature  of  1400  K  is  0.6061x10  N/m  and  ultimate 

6 

tensile  strength  of  SS  304  at  this  temperature  120./  x  10 

2 

N/m  as  given  in  Appendix  l). 

lhe  diameter  ratio  of  the  container  tube  is  obtained  from 

Appendix  t  as  D  /D.  1.0b  for  P  0.606  MPa  and  f  120./  MPa. 

oi  v  u 


8C 


An  available  standard  SS  '104  tube  stock  of  ?.??  cm  00,  1.89  cm  ID, 
and  IDS  mm  wall  thickness  is  selected  for  the  envelope.  Ihe  extra 
wall  thickness  provides  added  safely  for  this  design.  However,  the 
associated  heavy  thermal  mass  has  to  be  tolerated.  Ihe  flat  end  cap 
thickness  is  also  obtained  from  the  design  chart  of  Appendix  F  as 
t/D0  ;  0.0/S  for  the  above  vapor  pressure  and  strength  values. 
Therefore,  the  minimum  required  end  cap  thickness  is  t  0.07S  x 

7.72  cm  1.6/  nim.  Actual  thickness  used  is  6  mm  for  the 
conveniences  of  machining  and  welding. 

3.3.?  Wick  Design 

Ihe  wick  rorif  igurat  ion  chosen  is  a  version  ol  the  double 

walled  artery  concept  researched  by  the  author  in  the  recent  past 
110,11 j.  A  cross  sectional  view  ol  the  wick  geometry  is  illustrated 
in  Figure  ?9 .  In  general,  the  double  walled  wick  arrangement 
consists  of : 

(1)  A  cylindrical,  selectively  perforated  and  externally 
grooved  (axial)  inner  tube  wall; 

(?)  An  inner  annuli  screen  mesh  wick;  and 

(3)  An  outer  tube  envelope  with  internally  grooved 

evaporator  all  packaged  in  a  single  concentric 
assembly.  Ihe  inner  tube  contains  the  artery  transport 
groove  wick  and  firmly  presses  the  inner  annuli 
evaporator  supply  screen  wick  against  the  outer  tube. 
This  feature  ensures  that  the  screen  wick  remains  in  good 
mechanical  contact  will)  the  outer  tube  throughout  its 
operating  life,  thereby  enhancing  l  tie  radial  heat 
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transfer  in  the  evaporator  arid  condenser  sections.  It  is 
of  particular  importance  in  the  liquid  metal  heat  pipe 
evaporator  since  hot  spot  development  due  to  the  detach 
ment  of  screen  wick  from  the  wall  is  common  in  conven 
tional  screen  wick  designs.  Suitable  vapor  vent  cut  outs 
along  the  axial  grooves  of  the  inner  tube  provide  for  the 
vapor  injection  from  the  evaporator  and  the  vapor  suction 
into  the  condenser  screen  wick,  lhe  screen  wick  is  formed 
of  spirally  wrapped  multiple  layers  of  wire  inesh  with 
suitable  pore  size  to  give  the  required  capillary  head. 

lhe  special  feature  of  this  design  is  the  entrainment  -f ree 
vapor  transport  due  to  the  physically  separated  vapor  and  liquid 
flows  in  the  adiabatic  section.  Other  charact eri st ic  features 
uniquely  incorporated  in  the  present  wick  design  are: 

(1)  Single  groove  transport  artery  as  opposed  to  multiple 
grooves ; 

(?)  Screenless,  ventless,  and  groove  only  transport  section; 

(3)  Non  constant  groove  widths  for  evapora 1  or/condenser ;  and 
transport  sections. 

Inner  and  Outer  lubes:  lhe  inner  tube  diameter  and  the  wall 
thickness  are  determined  from  the  requirements  of  the  vapor  core 
diameter  and  the  artery  channel  size,  respectively.  Also,  the 
availability  of  standard  tube  stock  dictates  the  choice  of  the 
dimensions.  Hence,  two  standard  seamless  extruded  tube  sizes  were 
picked.  they  are,  ? .??  cm  00  and  1.65  mm  wall  outer  tube  and 
1.905  cm  0D  and  3  mm  wall  inner  lube,  lhe  )?.f  nun  10  of  the  inner 
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tube  more  than  satisfied  the  8.3?  mm  diameter  required  for  t tic  vapor 
core  and  the  3  mm  wall  was  adequate  to  machine  an  artery  groove  ol 
up  to  2.4  mm  (3/3?  inch)  square  (see  Section  3.3.3  for  artery 
dimension) . 

Screen  Wick:  The  screen  or  the  capillary  wick  radial 

thickness  depends  on  the  radial  temperature  differences  (AT) 

permissible  at  the  evaporator  and  condenser  sections.  larger  wick 
thickness  means  higher  thermal  resistance  and  consequently  larger 
Al  across  the  wick.  On  the  other  hand,  smaller  wick  thickness  may 
not  be  sufficient  to  hold  enough  working  fluid.  In  addition,  the 
main  parameter  controlling  the  heat  transfer  across  the  wick  is  the 
effective  thermal  conductivity  of  the  fluid  saturated  screen  wick. 
Several  analytical  and  empirical  models  are  available  in  the  heat 

pipe  literature  for  effective  thermal  conductivity.  Section  3.3.5 
deals  with  this  in  detail,  lhe  screen  wick  tube  diameter  ratio  can 
be  obtained  from  the  rad  i<il  heat  conduction  equation  for  the 

cylindrical  geometry  given  in  tq.  (38). 

2  it  1  k  AT 

Ur  <-  f  eff  r,s  (38) 

’  In  (d  /d.) 

'  o  is 

Qr  •-  ?000  W,  Lf  •  0.25  m,  k&^  39  W/m  K  (lor  sodium  saturated 

stainless  steel  wick  at  1000  K  using  series  parallel  model)  and 

AT  is  assumed  not  to  exceed  5  K.  It  may  be  noted  on  similar 
r,s 

calculations,  that  the  Al  across  the  outer  tube  wall  is  about 

15  K.  For  Al  =  5  K,  the  above  equation  gives  (d  /d  )  1.166 

r.s  or 
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and  knowing  d  1.89  cm,  d.  is  evaluated  as  1.6?  cm.  Hence, 

o  1 

the  total  thickness  of  the  screen  layers  is  0.13b  cm.  lhe  wire 

diameter  ol  the  40  x  40  cm  1  screen  mesh,  d  0.1143  mm. 

thickness  of  one  layer  of  screen  -  ?d  --  0.??8b  mm.  Number  of  layers 

of  screen  -  1.3b  mm/0.??86  mm  :  S.9.  In  order  to  be  consistent  with 

the  previous  copper  water  heat  pipe  design,  a  4  laye1'  screen  tube 

was  selected.  Ihis  incidentally  reduced  the  Al  to  3.3  K. 

r ,  s 

Capillary  Pumping  Head:  In  a  composite  wick,  the  effective 
capillary  radius  is  that  of  the  pumping  wick  (screen)  and  hence,  the 
capi I lary  radius, 


I 

?N 


(.39) 


1  0 . 000b  in.  =  0 . 1 ?  1  mm 


?x  1 00 


lhe  maximum  capillary  pumping  head  provided  by  r^  is  given  as: 


Po¬ 


em 


(40) 


At  400°C , 


cm 


?x_Q . 159  ,  ?b03 . 9  N/m 

0.1?/  x  10'3 


I  rom  the  thermophysical  property  data  given  in  lahlc  4  for  sodium, 
the  capillary  pumping  head  at  various  temperatures  can  be 
calculated.  lhe  normal  and  axial  hydrostatic  pressure  drops  given 


TABU  4.  f HEHMOPHYS 1  CAL  PROPtRlY  DATA  Of  SODIUM 
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by  Iqs.  (41)  <ind  (4?)  are  *>ubi r*ic t **tl  from  I*  to  obtain  maximum 
pumpin']  head  available. 


AP 

I 


p  g  (1  COS  <1; 
w  y 


(41) 


At  4 00 ° C 


AP  86S.9  x  9.81  x  I  .89  x  10  2  x  1  1  '>8 .  /  N/m? 

i 

AP^  p^g  sin  41  0  (for  horizontal  heat  (4?) 

pipe  operat ion) 

*  pm  *  cm  (43) 

Table  b  lists  these  values  for  various  t  empera  l  ures . 


lABIf  ‘j.  CAPIILARY  PUMPING  HEAD  AVAIt  ABLt 
WIIH  40x40  cm'1  MESH  AND  SODIUM 


Temperature 

P 

cm 

AP 

L 

P 

pm 

(°C) 

(N/m2) 

(N/m2) 

(N/m2) 

400 

?504 

T  58 . 7 

?34S . 3 

bOO 

?3  78 

153.5 

???4 . b 

GOO 

??3b 

149.3 

?08b .  / 

/  00 

?094 

141  .6 

19b?. 4 

BOO 

193/ 

140.4 

1796.6 

900 

17/9 

138.? 

1 640 . 8 

1000 

1  b38 

134.5 

1 b03 . b 
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3  .  .  3  Artery  Design 

The  artery  channel  design  procedure  followed  here  is  the  saun¬ 
as  that  used  in  the  2  m  copper  water  heal  pipe  |1,11,55].  Ihe 
groove  sues  used  for  water  heat  pipe  should  work  for  the  sodium 
heat  pipe  also  since  the  liquid  transport  factor  at  t he  rev, per  t  ive 
peak  operating  temperature  for  sodium  is  A.?  times  greater  than  that 
for  water.  In  principle,  even  larger  grooves  than  those  used  in 
water  heat  pipe  may  be  employed  for  sodium  heat  pipe.  However,  for 
the  sake  of  comparison  and  other  considerations  such  as  wettability, 
wetting  angle,  etc.,  the  dimensions  of  grooves  are  retained  the 
same.  Accordingly,  thf  number  of  grooves  on  the  evaporator  and 

condenser  sections  are  24  and  their  si/es  are  0.79  x  1.4  mm.  On  the 
transport  section,  only  one  2.4  x  2.4  mm  groove  is  necessary  as  a 
liquid  transport  artery.  Incidentally,  this  single  groove 

arrangement  helps  in  testing  f fie  heat  pipe  for  artery  orientation 
effects  (with  respect  to  gravity  vector)  and  acts  as  a  scaled  down 
model  for  easy  low  power  testing.  Ihe  approximate  scale  down 

criterion  used  is  as  follows: 

theoretical  transport  capacity  of  the  copper  water  heat  pipe  with  12 

artery  grooves  -  5142  W 

Capacity  for  1  groove  5142/12  428.5  W 

for  sodium  heat  pipe  with  I  artery  groove  of  the  same  si/e,  the 

capacity  --  428.5  x  4.2  (transport  factor  ratio) 

=  1800  W 

Ihis  is  close  to  the  nominal  2000  W  capacity  set  as  the  design  goal. 
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Ihe  Inner  lube  <jroove  and  vapor  vent  details  are  shown  in 
Figure  30.  Alternate  grooves  in  1  lie  evaporator  are  vented  whereas 
all  t he  grooves  in  the  condenser  are  vented.  No  vent  slots  are 
provided  in  the  adiabatic.  A  short  capillary  insert  is  installed 
within  the  adiabatic  groove  at  the  evaporator  exit. 

SeJX  Prirning  Requirement:  Arterial  heat  pipes  should  be 

capable  of  priming  with  the  working  fluid  within  the  capillary 
spaces  at  all  times.  11  an  artery  becomes  depleted  of  working  fluid 
due  to  external  influences  such  as  evaporator  dryout  or  unfavorable 
tilt  in  1  g  environment,  it  should  be  able  to  reprime  automatically 
once  the  external  influence  is  removed.  Such  an  artery  is  called  a 
self  priming  artery.  Self  priming  requirement  can  be  satisfied  by 
limiting  the  diameter  of  the  heat  pipe  to  within  the  maximum 
capillary  height  potential  (It)  given  by  the  relation, 


P  g«  =  .L 

l 


cos  0 
r 

c 


(44) 


lor  comparison  purposes,  H  has  been  calculated  for  several  wick 
structures  as  given  in  Table  6.  In  the  present  design,  the 
circumferential  flu’d  distribution  is  accomplished  by  the  screen 
wick  at  the  evaporator  and  condenser.  The  capillary  grooves  prime 
from  the  screen  and  the  larger  (?.4  x  ?.4  mm)  artery  groove  primes 
from  the  0.79  x  0.79  mm  grooves  in  the  condenser  section. 

Artery  Dimension:  The  artery  in  a  heat  pipe  could  be  a  screen 
tube,  groove,  a’ie-.ar  ^pace,  or  a  circular  tube.  In  the  double  wall 
design,  it  is  -./.enient  to  machine'  square  groove  on  the  inner  tube. 
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TABIE  6.  COMPARISON  OE  WICK ING  Ht I GH I S  EUR  SODIUM  IN  VARIOUS  WICKS 
(SCREEN.  GROOVE.  ANO  CIRCULAR  TUBE) 


Theoretical  Wicking  Height,  H  (cm) 


10x40  cm  1 
wire  mesh 


29.27 

28.30 

27.96 

26.0/ 

24.34 


1 . 59x1 . 59  mm 
groove 

r  1 . 59  mm 

L 

2.38x2.38  mm 
groove 

r  -  2.38  mm 

c 

2.34 

1.56 

2.26 

1.51 

2.24 

1.49 

2.09 

1.39 

1  .95 

1  29 

1  .84 

1.23 

In  order  to  arrive  at  the  width  of  the  groove,  .1  conservative  appro. u  h 
is  used  wherein  l  tie  groove  should  he  able  to  prime  to  the  height  at 
least  equal  to  t  fie  diameter  of  t  fie  he.it  pipe.  table  7  gives  I  tie 
artery  widths  at  various  temperatures.  At  /00°C,  the  calculated  value 

of  w  is  2.69  mm  and  the  chosen  width  of  2.4  mm  is  well  within  this 

estimate.  An  optimistic  estimate  based  on  the  priming  height  equal  to 

the  width  or  depth  of  the  groove  itself  is  also  made  and  the  numbers 

are  given  in  table  /. 


7A8IF.  /.  SQUARf  OROOVf  ARIfRY  SI /h  I  OK  SODIUM 


t  emperat  ure 

Artery  Width  w  -  A 

(°C) 

Conservative 

Optimist ic 

Design 

Desiyn 

(mm) 

(mm) 

O 

O 

2.85 

6.09 

600 

2.73 

5.99 

700 

2.69 

6.96 

800 

2.4/ 

6.  /S 

900 

2.2/ 

6.56 

1000 

2.12 

5.41 

Equations  used: 


1.  Conservative  Approach: 


80^  cos  0 
g(pt-pv) 


(45) 


2.  Optimistic  Approach: 


w 


2a ^  cos  0 
9(Pt~Pv) 


(46) 


9? 


3.4.  lumped  Ihermal  Capacity  Analyses 

3.4. !  Conf  i (jurat  ion  Summary 

•  Design  details  as  in  Section  3.3.1 

•  I  ,  /I  ./I  33//9/91  cm  I  „  ?03  cm 

r  A  l  P 

D  7.77  cm 
o 

•  Single  groove  artery  in  adiabatic;  multiple  groove 
screen  covered  evaporator  and  condenser 

•  Reservoir  wick  in  evaporator 

•  fxcess  metal  in  inner  tube  cut  out  in  adiabatic 

•  Sodium  working  fluid;  argon  as  the  noncondensible  gas 

•  l ro/en  start  from  room  temperature  (1q) 

3.4?  Objective 

lo  find  ttie  maximum  length  of  adiabatic  section  that  can  be 
tolerated  in  the  7  m  long  heat  pipe  without  drying  out  t tic 
evaporator.  Also  determine  the  si/e  of  the  liquid  reservoir  wick  at 
the  evaporator  end.  The  sketch  of  the  heat  pipe  is  shown  in 
I  igure  3 1  . 

3.4.3  Assumptions 

(1)  All  the  sodium  in  the  evaporator  is  transferred  in  the  form  of 

saturated  vapor  (no  superheat)  to  thaw  the  adiabatic  section 

and  (o  move  the  hot  z one  to  the  beginning  of  the  condenser. 

(?)  No  liquid  return  to  the  evaporator  during  startup. 

(3)  lor  a  given  noncondensible  gas  (NCG)  charge  pressure  (P..), 

there  is  a  fixed  evaporator  temperature  (1  )  to  which  the 

H 

evaporator  has  lo  be  heated  up  before  the  hot  /one  moved  to  the 
adiabatic  section. 
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RESERVOIR 

WICK 


LIQUID  RETURN 


SINGLE  CHANNEL 
ARTESL  NON  CONDENSIBLE 


Figure  31.  Single  Channel  Liquid  Metal  DWAHP  with  Reservoir  Wick 
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(4)  No  heal  losses  in  the  adiabatic  section. 

(b)  Only  vapor  diffusion  is  considered  for  energy  transport.  No 


axial  conduction  through  wall. 

(6)  Ihe  initial  power  input  is  less  than  the  sonic  limit 

corresponding  to  the  desired  startup  temperature  (T(<). 

( /)  Ihe  reservoir  wick  is  in  good  contact  with  the  evaporator  wick. 

and  there  is  fluid  communication  between  the  two  wicks. 

(8)  Sodium  working  fluid  is  uniformly  distributed  within  the  wick 
structures  (screen  and  grooves). 

3.4.4  Energy  Balance 

lor  the  hot  front  to  move  from  evaporator  exit  to  the 
beginning  of  the  condenser,  the  energy  supplied  externally  to  the 
evaporator  has  to  vaporize  the  sodium  in  the  evaporator.  Ihe  latent 
heal  of  vaporization  stored  in  the  vapor  is  transferred  to  the 
adiabatic  section  by  condensation  of  the  vapor.  H  lumped  energy 
balance  under  steady  state  can  be  written  as  follows: 


|  Energy  input  » 

I  to  Evaporator  (  * 


Energy  carried  by 

the  saturated  sodium 

vapor  from  evaporator 

2 

to  adiabatic  at 

temperature 

Energy  required  to 
raise  the  adiabatic 
section  from  T0  to 

th 


(47) 


Qi  4  mNa  W  4  CA(TH  -  V 


(48) 


9‘> 


For  determ  i  n  i  ng  the  exa<  t  c|  u.i  n  t  i  t  ies  of  '.odium  fill  arid  the  possible 
maximum  1^,  F  rj .  (4H)  <  an  he  written  a-.: 


V,  V, 


V1,, 


I  ) 
o 


fxpandinq  the  terms  we  can  write. 


m 


Na 


+•  m 


Na 


Evap 

screen 


+  m 


Na 


Evap 

groove 


x  (T  ) 
Na  H 


Reservoir 
wick 


=  (msscp,ssLA  f  mNacp,NaLA)(TH  '  [o) 


(49) 


where 


m 


Na 


m 


x 

C 

C 

I. 

1 

1 


ss 

Na" 
P .  Na 
p,ss 
A 
H 


o 


mass  of  sodium  in  the  wick  per  unit  length 
mass  of  stainless  steel  wall  per  unit  length 
heat  of  vaporization  for  sodium  at  temperature  I 
specific  heat  capacity  of  sodium  liquid 
specific  heat  capacity  of  stainless  steel 
adiabatic  length 
hot  zone  temperature 

room  temperature 


The  mass  terms  are  written  and  simplified  by  substituting  t  fie 
cross  sectional  areas  and  wick  porosities  as  follows: 
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m 


N  a 


f  vap 
screen 


A  I  r  £  p_ 
w  F.  vl 


0.:^s:is  i(,.# 


m 


NA 


I  vap 
groove 


V  1  I "ft 


0.?66!3  Ip 


m 


Na 


V'kVi 


Reservo  i  r 
wick 


LRpa 


the  IMS  of  fq.  (49)  becomes, 

( ( 0 . 3?S3‘j  i  0.7661 3)1..  •  7.11 674  l._)p  (IJ 

f.  R  8  Nj  H 

(0.S9M8  l.t  .  2.1I6?9  I  R)  VNa<  V 
(O.M829  I  y  i  I  .9661B  I^X^Cy 

3 

I  fie  density  of  sodium  p^  ~  0.9?/  g/cm  correspond ing  Lo  the 
phase  change  temperature  of  97.83°C  has  been  used  since  only  that 
much  of  sodium  can  he  stored  in  the  wick  pores.  Ihe  latent  heat  of 

sodium  x  depends  on  the  startup  temperature  I  which  can  he 

Nd  II 

varying  from  4()0°l  to  H00°C. 


9/ 


Similarly,  subs  t  i  tut  i  mj  the  spec  i  I  i  c  float  and  mass  values  <1  rid 
simplifying,  we  obtain: 

1  he  RHS  of  fq.  (49)  (5.91343  i  0. 0/5903)  I. (I,,  1  ) 

A  H  o 

5.98934  L  ( I  1  ) 

AH  o 

Now,  fq.  (49)  can  be  written  as: 

( 0 . 54879  I  (  f  1.9618  LR)  X^y  5.98934  lA(lH  y  (50) 

(or  sodium  vapor,  the  transitional  flow  regime  (from  free  molecular 
flow)  to  continuum  (low  regime  occurs  at  340°C.  In  a  gas  filled 
heat  pipe  such  as  our  present  design,  the  vapor  flow  attains 
(presumably)  continuum  flow  regime  before  the  hot  front  a(  the 
gas  vapor  boundary  starts  moving  away  from  the  evaporator.  Above 
5?0°C,  the  sodium  heat  pipe  switches  to  isothermal  operations. 

Substituting  lH  340°C,  I  >*5°C  and  y(340)  ;  ^393  J/g 
in  F.q.  (50)  we  obtain, 

0 . 548?9  I  f  1.9618  I „  -  0.4^946  I.  (51) 

F.  K  A 

From  the  design  details  we  have, 

I.,  e  I  »  I.  1 1 2  cm  (5?) 

FRA 
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for  practical  considerations  we  can  assume  that, 


I  heref  ore , 


or 

and 

Solving  l qs 


lo  give  a 
assume  1 


'l  ?,R  <53) 

tqs.  ('jl)  and  (S3)  become 

1  .09668  L  t  1.9618  l.„  0.43946  L. 

R  R  A 

3.06838  I .  -  0.43946  1 ,  (64) 

R  A 

3  r  i  I.  113  (66) 

R  A 

.  (64)  arid  (66) 

L_  =  11 .06  cm  > 

K  I 

L.  =  78.80  cm  >  for  T  =  340°C 
A  I  n 

Lf  =  22.12  cm  ' 

factor  of  safety  on  Ihe  startup  temperature  1  let  us 

H 

f)  -  400°C 

\  ( 400°C)  4380  .!/g  ;  I,,  I  =  376 

Na  Ho 


99 


fq.  (50)  becomes, 


0 . 54829  L  »  1.9618  I  -  0.5f?/flh  I  (51)’ 

r.  K  A 

Using  Eqs.  (52)  and  (53)  as  before,  we  obtain 

L  =  1 2 . 50  cm  , 

K  I 

La  =  74.50  cm  / 

>  for  Th  =  400°C 
lE  =  25.00  cm  I 

Lc  =  91 .00  cm 

We  select  this  result  for  siting  the  heat  pipe.  It  may  be  not ed 
that  even  though  this  analysis  does  not  predict  the  movement  of  the 
hot  front  as  a  function  of  time,  it  is  accurate  enough  to  tel)  that 
the  front  can  reach  the  condenser  before  the  evaporator  dried  out. 

3 . 5  Pe r f o rmance  Limit s 

The  performance  of  the  heat  pipe  is  judged  by  its  ability  to 
transport  a  desired  amount  of  thermal  energy  supplied  at  the 
evaporator  end  to  the  condenser  end  with  a  minimum  temperature  drop 
and  without  causing  a  dryout.  The  temperature  dependent  heat 
transport  capacity  of  t tie  heat  pipe  under  steady  state  conditions  is 
limited  by  several  other  factors  such  as  sonic  limitation,  capillary 
limitation,  boiling  limitation,  entrainment  limitation,  and 
condenser  radiation  limitation.  A  quantitative  estimate  of  each  of 
the  limitations  is  as  follows. 
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3.5.1  Condenser  limitation 

This  is  an  externa]  limitation  due  to  the  radiative  coupling 
of  the  condenser  to  Ihe  cold  wall  of  the  chamber,  lhe  net  radiant 
energy  exchange  from  the  heat  pipe  condenser  to  the  cold  chamber 
wall  is  given  by  Stefan's  formula. 

Qr  •-  Accco(TA  T*)  (56) 

1?  ?  4 

Here,  Stefan  Boltzmann's  constant,  o  •-  5.6?  x  10  W/cm  K  ; 

emissivity  of  the  oxidized  stainless  steel,  cc  0.6;  area  of 

2 

the  condenser,  Aj.  ir  0^.  1.^  -  4?5.6  cm  for  0^  2.22  cm 

and  l  -  60.96  cm;  the  chamber  temperature,  1  is  assumed  to  be 
C  s 

the  average  of  the  25°C  cooling  water  inlet  and  80°C  outlet 
temperature  or  I  -  52.5°C;  the  heat  pipe  operating  temperature,  1 
ran  vary  from  400  1000°C.  0p  has  been  calculated  for  l_c  -  30, 
60,  and  90  cm  and  tabulated  as  shown  In  Table  8.  This  limitation  is 
sketched  in  figure  32  along  with  other  limitations. 

3.5.?  Sonic  Li  ml tat i on 

Sonic  limitation  generally  occurs  when  the  heat  pipe  is 
operating  at  low  vapor  densities  (that  is  at  low  temperatures)  and 
high  vapor  velocities.  Critical  and  choked  flow  conditions  exist 
when  the  vapor  velocity  at  the  evaporator  exit  becomes  sonic 
(H^  1.0).  In  liquid  metal  heat  pipes,  the  sonic  limitation  is 

very  significant  deling  the  startup  phase.  At  this  limit  there  is  a 
maximum  axial  heat  transport  rate  due  to  the  choked  flow  and  a  fixed 
axial  temperature  drop  along  the  evaporator,  which  is  associated 
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HEAT  TRANSPORT  RATE,  WATTS 


400  500  600  700  800  900  1000 


TEMPERATURE,  °C 


Figure  32.  Heat  Transport  Limits  Q  Versus  T 


10? 


1  ABLE  8.  CONDENSER  RAD  l  A  I  ION  I.IMI  I  AT  ION 


temperature 

Radlat  1 

on  1 imit,  0 

r  (W) 

(°C) 

(K) 

Lc  =  30  cm 

L^  -  60  cm 

90  cm 

400 

6/3 

140 

280 

421 

500 

7/3 

250 

500 

751 

600 

873 

413 

825 

1237 

700 

9/3 

640 

1281 

1922 

800 

1073 

951 

1903 

2854 

900 

11/3 

1362 

_  _ 

1903 

2854 

1000 

1  273 

1893 

3786 

5680 

.  _  _ _ 

_ _  _ ..... 

with  any  given  evaporator  entrance  temperature.  Increasing  the  heat 
rejection  rate  heyond  the  sonic  limit  lowers  the  condenser 
temperature.  Induces  supersonic  vapor  flow,  and  creates  very  large 
axial  temperature  gradients  along  the  pipe  but  it  does  not  increase 
heat  transport  rate.  Consequently,  near  Isothermal  operation  is  not 
possible.  Levy's  equation  is  used  to  calculate  the  sonic  limit 
[  4  J .  I  hat  is. 


.0.5 


'sonic 


,  .  y  R  T 

A  o  X  |  v  v  0 

V  0 


2(Y/D 


(Levy's  Equation) 


(57) 
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Subscripts  v  and  o  refer  to  the  static  and  stagnation  states  of  t tie 


vapor,  respectively,  for  sodium  1.6/  and  Rv  R/M  361  J/k<)  K. 

?  4  ? 

The  vapor  core  cross  sectional  area  Av  ---  *d./4  1.267  x  10  m  . 

Usinq  the  property  table  for  sodium,  Q  ,  has  been  calculated  at 
*  sonic 

various  temperatures. 


Temperature  (°C) 

‘’sonic  <“> 

40? 

36.8 

452 

122.3 

502 

343.9 

552 

849.3 

602 

iD 

<T> 

CO 

r— 

Sonic  limit  curve  is  plotted  in  Figure  3?  along  with  the  other 
limitations.  It  can  be  observed  that  the  sonic  limit  determines  the 
present  sodium  heat  pipe  transport  capacity  In  the  400-560°C  range. 

Ihe  presence  of  NCG  in  the  heat  pipe  may  affect  the  sonic 
limit  calculations  based  on  Levy's  equation.  Ihe  effect  is 
favorable  to  a  frozen  startup  situation  since  the  NCG  slows  down  the 
vapor  molecules  and  possibly  prevents  them  from  attaining  sonic  or 
supersonic  velocity  until  the  vapor  pressure  becomes  equal  to  the 
NCG  charge  pressure.  Once  the  vapor  pressure  exceeds  the  NCG  charge 
pressure,  then  the  sonic  limit  by  Levy's  equation  will  be 
applicable.  Any  variation  in  temperature  varies  the  length  of 
active  condenser. 
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3.5.3  Cajpi  1  la ry  1.1ml  ta t  i on 

Capillary  wick  priming  considerations:  In  composite  wicks, 
capillary  pumping  Is  achieved  by  surface  tension  forces  created  by 
the  small  pores  of  the  screen  wick  at  the  liquid  vapor  interface  in 
the  evaporator  where  the  evaporation  of  the  liquid  takes  place, 
liquid  from  the  condenser  section  flows  through  the  artery  channel 
into  the  evaporator  grooves  and  then  radially  into  the  screen  layers 
to  make  up  for  the  depletion  of  fluid  by  evaporation,  this  process 
of  saturating  the  wick  Is  called  the  capillary  priming  and  this  has 
to  be  maintained  at  all  times  of  the  heat  pipe  working. 

Capillary  inserts:  As  evidenced  by  the  low  temperature  visual 
tests  on  glass  walled  copper  water  heat  pipe,  a  possible  vapor 
back  flow  into  the  liquid  channel  is  prevented  using  a  screen  insert 
at  the  evaporator  exit.  However,  this  Insert  offers  extra 
resistance  to  the  return  liquid  flow. 

L i ay_! d  flow  Resistance:  The  pressure  drop  associated  with  the 
liquid  flow  In  various  segments  of  the  wick  structure  Is  modeled  as 
a  set  of  series  and  parallel  resistance  paths.  The  total  pressure 
drop  is  calculated  from  these  resistances.  Eq.  (58)  gives  the 
liquid  flow  pressure  drop  In  the  wick  structure  for  elemental  length 
dx . 


f-  0  t  p  g  sin  <|r 

9.  I 


Ij  is  the  liquid  friction  coefficient  defined  as: 


or  J  .  (_t) 

K  A  p.\  K  A  ,  \ 

w  wl  w  w 


(58) 


(59) 
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The  capillary  transport  capacity  of  the  heat  pipe  is  inversely 
proportional  to  this  factor  F  .  lor  horizontal  operation  of  t tie 
heat  pipe,  ^  -  0  and  F.q.  (58)  simplifies  to: 


Upon  integrating  over  the  desired  tenth  Ax  of  flow  in  wicks, 


(60) 


-AP^  -  (FftAx)Q  (61) 

The  above  equation  is  analogous  to  an  electrical  resistance  circuit 
where  AP  corresponds  to  voltage  drop,  Q  corresponds  to 
current,  and  F^Ax  corresponds  to  resistance.  Hence,  the  liquid 
flow  resistance  due  to  the  various  segments  of  the  wick  structure 
can  be  added  in  series  to  get  the  total  resistance. 


(U) 


d  ef  f 


(F  Ax)  ,  , 
v  d  ' total 


( FaAr)evaporator 
screen 
radial 


c  ( F^ff/2) evaporator  screen 
ana  groove 
axial  effective 


f  (FdLl ) inserts 
axial 


f  i A^transport  groove 
axial 


f  ( ^dLc/^)condenser  screen  *  ( ^dAr)condenser 
and  groove  screen 

axial  effective  radial 

(6?) 


► 
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lhe  individual  resistances  which 

are 

f  unctions 

of 

the 

wick 

mechanlca 1 

properties  (permeability, 

a  rea 

,  etc.) 

and 

t  he 

fluid 

propert ies 

(kinematic  viscosity  and  latent 

heat)  are 

computed . 

lhe 

detai Is  are 

given  in  Appendix  G.  lhe 

result 

is, 

IP  v  p 

(FaL)eff  15. 2085  x  10  (VA)  [ ( N/m  )/W]  (63) 

Approximate  Capillary  limit  Calculation: 

An  estimate  of  the  capillary  limit  neglect  in<j  the  vapor  Mow 

friction  effects  (I-  )  is  made  in  order  to  obtain  a  first  cut 

v 

approxlmat Ion  using  fq.  (64) 


Pm 

Veff 


and  the  values  are  listed  in  lable  9. 

Vapor  FJgw^ Friction  Factor  (F 

lhe  applicable  equations  from  Reference  4  are. 


(64) 


I 

v 


f  Re  p 
v  y  v 

u  A  p  X 
h,  v  v'v 


?  r 


where  f  Re  16  if  Re  *  ?300  and  M  s;  0.2 

V  V  V  V 


Re  -  ?rhtv  Q 

v  A  v  X 
V  V 


M 

V 


A  p  X 
v  v 


T  R  T 

V  V  V 


(65) 


(66) 


(67) 
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The  vapor  core  dimensions  are. 


rh  v  =  rv  0 . 635  x  10  1  ... 

2  4  2 

A  ir r  I  .  266/  x  10  m 
v  v 

Re  100.255  0 

v  . 

pvX 

Re  is  calculated  in  terms  of  Q  for  various  temperatures  and 

listed  in  Table  10.  Ihe  limit  of  beat  transport  rate  allowed  to 

satisfy  laminar  vapor  I  low  condition  is  obtained  by  setting  the 

condition.  Re  *  2300.  Also,  M  is  calculated  in  terms  of  Q 

v  v 

by  substituting  1.6 1  and  R^  361  J/kg  K  for  sodium  in 

Fq.  (6/)  and  listed  for  different  temperatures  in  Table  10. 
Initially,  Q  is  unknown  to  determine  Re^  and  arid  hence  the 
approximate  estimates  of  Qcl  (Table  9)  is  used  to  start  the 
calculations.  The  conditions, 

s  0.2  (Incompressible  flow) 

*  0.2  (Compressible  flow) 

are  used  to  determine  the  compressibility  of  the  vapor.  Based  on 
the  sonic  limit  and  the  operating  power  limit  known,  the  vapor  flow 
is  lound  to  be  laminar  in  the  entire  range  of  operating  temperature 
and  compressible  up  to  600°C  and  incompressible  at  temperatures 
above 

600"C.  Ihe  applicable  equations  for  F^  are, 
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TABLE  10.  VAPOR  fLOW  FRICTION  COEFFICIENT  VERSUS  I LMPERA1 URL 
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8u 


r.  A  p  X 
h,v  v'v 


for 


Re  i  2300 

v 

M  S  0.2 

v 


(68) 

(Laminar  and  incompressible) 


8p. 


rT  A  p  x 
h,v  vrv 


Y  -1  2 

1  +  „v__  Hv 


-0.5 


for 


Re  s.  2300 
v 

M  >0.2 
v 


(69) 


(Laminar  and  compressible) 


Qr?  Ppm .  (70) 

<^Fv»eff 

After  calculating  I  .  the  capillary  limit  Q  „  is  calculated 

v  c2 

using  fq.  (10)  and  L  -  1.41  m.  1  he  results  are  given  in 

lable  9  and  figure  32.  It  should  be  noted  that  the  vapor  flow 
friction  is  very  high  during  the  startup  temperature  regime  (up  to 
540"C) . 

3.5.4  Bolling  Limitation 

Ihe  boiling  limitation  is  a  limitation  of  the  radial  heat  flux 
density  at  the  evaporator  due  to  the  format  ion  of  vapor  bubbles 
within  the  wick  structure.  Ihe  bubble  formation  is  undesirable 

because  it  can  cause  hot  spots  and  obstruct  fluid  circulation.  In 
liquid  metal  heat  pipes,  this  limitation  is  not  serious  since  liquid 
metals  have  high  thermal  conductivities.  Ihe  series  para  I  le 1  model 

given  by  Chi  [4J  for  the  effective  thermal  conductivity  (keff)  °f 
sodium  saturated  stainless  steel  wick  is  given  in  tq.  (71). 


k 


ef  f 


y(ki '  y  ■ (1  e)(ki_  ■  v] 
t('kl  rks)T (f ■'•"ej'('k~"ksjr 


(71) 


1 1  1 


The  thermal  conductivity  of  steel,  k^  is  1  / .  3  W/m  K  in  t  ho  nominal 
operating  temperature  range  of  r)0(J  I000°C  while  that  of  sodium 
liquid,  varies  from  70  <19  W/m  K.  Suhst  i  tut  ing  t  0.6/89 

for  the  <10x40  cm  *  wire  screen  with  0.1  M3  mm  wire  diameter  in 
Fq .  ( 71 ) .  one  obtains 


-  44.42 

W/m 

K 

at 

600°C 

or  34.27 

W/m 

K 

at 

I000°C 

Boiling  limit  is  given  by 


2»  I  k  t 
_  f  e  f  f  V  ( 2  a 

\  p  in  ( r  /r . )  r 
v  v  o  1  n 


»*C> 


(72) 


The  nucleation  radius,  r  ,  is  determined  using  the  nucleation 

n 

theory  of  Rohsenow  and  Bergles  [M]  for  smooth  surfaces  represented 
in  tq.  (73) . 


?a  T  .  k. ( V 
( _ sat  d'  v 


VJ  o  s 
) 


U  3) 


For  sodium  heat  pipes,  the  radial  heat  flux  density  values  reported 

in  the  literature  range  from  200-400  W/cn/.  Assuming  the  minimum 

value  of  200  W/cm2  for  qr  ,  rn  for  600  and  1000°C  are 

-4  4 

calculated  as  7.899x10  m  and  0.5046x10  m  respectively. 

These  values  are  much  higher  than  or  close  to  t fie  wick  pore  sizes  of 

4  1  4 

1.27x10  m  for  40x40  cm  screen  mesh  and  1.466x10  m  for 

the  triangular  grooves  on  the  evaporator  inner  wall.  lor  lower 
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radial  heat  fluxes  of  about  10  20  W/cm  applicable  to  the  present 

heat  pipe,  r^  predicted  by  Eq.  ( /  3 )  will  be  very  large.  It  is 

recommended  by  Marcus  [14]  that  this  equation  not  be  used  in  such 

cases.  Silverstein  [56|  has  suggested  that  the  effective  nucleation 

site  radius  in  properly  wetted  liquid  metal  heat  pipes  can  be  on  the 

order  of  10  ^  m  or  less.  An  upper  limit  of  r  3  to  7  x  10  ^ 

m  corresponding  to  a  hurnout  heat  flux  of  287  W/cn/  is  recommended 

by  him  [56].  Hence,  the  conservative  values  recommended  by  Chi  [4], 

r  2.54  x  10  m  for  gas  loaded  heat  pipe  and  r  2.54  x 

n  3  n 

10  m  for  conventional  heat  pipe  have  been  considered.  Using 
these  values  and  setting  I  0.25  m  and  assuming  that  the  pipe  is 

operating  at  capillary  limit  (P^  -  P  ),  is  calculated  from 

Eq.  (72)  at  various  temperatures  and  listed  in  Table  11.  for 

convertional  mode  operation,  the  boiling  limit  is  very  very  high 
whereas  for  gas  loaded  mode  the  limit  drops  sharply  at  and  above 
900°C . 

fxper [mental  Radial  .Heat  Flux  Oa  t  a :  The  experimental  data 
reported  in  the  literature  on  the  sodium  heat  pipe  radial  heat  flux 
density  for  various  wick  structures  are  summarized  below  [6,8]. 

Temperature  Radial  Heat  Flux 


Wick 

(°C) 

(W/cm^) 

1.  Stainless  Steel  Mesh 

760 

230 

2.  Various  Wicks 

850  950 

200  400 

3.  Stainless  Steel  Mesh 
(3  layer;  390  »im) 

925 

214 

4.  Stainless  Steel  Mesh 

lib 

1250 
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1  ABl  F  II.  BOI I  ING  LIMII  AN0CKM1CAI  ItMi'IIMIUlU  1)11  I  I  Kl  NCI 


temperature,  1^ 

(°C) 

Wick 

tf  f ec  t  i  ve 
Therma 1 
Conductivity, 
keff.  (W/m  K) 

Bolling  1  imi 

Convent iona 1 
Mode 

t.  Qb  (W) 

Gas  1 oaded 
Mode 

Al  .. 
cr  it 

(K) 

500 

44.42 

12. 24x10** 

98,131 

560 

600 

41.81 

2.908x10** 

23,316 

142 

700 

39.98 

0. /84/xlO6 

6,2/8 

39.9 

800 

38.53 

0.2959x10** 

2,368 

15.6 

900 

36.36 

0 . 1 248x10** 

999 

6.9/ 

1000 

34.2/ 

0 . 5508x10"* 

441 

3.27 

.  . 

....  .  .  _ 

...  _ _ 

.  .  .  . 

Boiling  of  sodium  in  sodium  hoot  pipes  is  impeded  as  a  result  of 
high  thermal  conductivity  of  1  tie  liquid  metal.  Fvaporation  ol  the 
working  f  luid  takes  place  at  the  free  surface  of  the  wick  pores. 
Assuming  a  sate  limit  of  100  W/cm^  radial  flux,  toe  radial  limit 
for  the  present  heat  pipe  is  obtained  as. 


0  .  ,  ir  0  L  q 

radial  F  F.  r 


(/6) 


2  4 

w  x  2.2?  x  10  x  0.25  x  100  x  10 


17,435  W 


However,  this  limit  is  too  high  to  be  encountered  since  the 

capillary  limitation  will  restrict  the  I  low  of  liquid  return  to  the 

evaporator.  Radial  heat  flux  density  expected  here  is  a  maximum  of 
2  2 

13  W/cm  based  on  Af  -  13/  cm  and  IJ  •  I  797  W  at  H00°C. 

3 . . ')  temperature  Drop 

Steady  state  temperture  differences  across  various  segments 
(pipe  wall,  wick,  and  vapor  p.issage)  of  the  heat  pipe  during  normal 
mode  of  operation  can  be  computed  using  Fourier's  law  for  conduction 
through  the  pipe  wall  or  wick  and  Clausius  Clapeyron  relation  for 
vapor  flow  [1 ] . 


1.  Pipe  wall  at  evaporator: 


ln(r  /r . ) 

A I  p  t  -  Q.  1  Q 

2w  I  £kp 


2.  Wick  at  evaporator: 


Hn(r./r  ) 

_ 1  .  S  Q 

2v  L.  k 
1.  ef  f 


3.  Vapor  passage: 


1  (P  P  ) 

Af  v  V  V,l  V.C  g 

P  X  .1  0 
V 


iy'v  (  f  '  '  '  C) 

(or  in  terms  of  lv)  h  A  6  q 

J 


4.  Wick  at  condenser: 


ln(r./r  ) 

Aiw  c  1  s  0 


5.  Pipe  wall  at  condenser: 


in{ r  /r .) 

ATp  c  .  o  i  .  q 

2.  Ickp 


Maximum  Al  is  expected  at  the  highest  power  transporting 
conditions  which  is  1/97  W  at  800°C  (lable  9).  The  A 1 1 s  are 
computed  using  the  physical  dimensions  and  the  property  values  at 
800°C  (table  4)  and  the  results  are  as  follows. 


10.83  C 


A 1  ,  3.9/  C 

W.t 


At  --  0.33  C 

v 


0.82  C 


2.98  C 


the  total  temperature  drop  across  the  heat  pipe  is  1/.82  C.  In 
gas  loaded  mode,  the  condenser  At  ’  s  will  vary  depending  upon  the 
active  condenser  length. 

3.b.b  f ntrainment  limitation 

In  conventional  screen  or  groove  wick  heat  pipes,  the 
counter  flowing  liquid  and  vapor  are  in  direct  shear  against  each 


1  lb 


other  and  this  causes  the  vapor  to  entrain  some  liquid  along  with  it 
causing  a  depletion  of  return  liquid  supply  to  the  evaporator.  Ihe 
heal  transport  limit  caused  by  this  process  is  called  the 
entrainment  limitation.  Here,  in  the  double  wa II  wick  structure 
there  is  no  direct  shear  between  the  liquid  and  vapor  flows 
par  icularyly  in  the  adiabatic  section.  Hence,  there  is  no 
entrainment  limitation  for  the  DWAHP. 

3.5.7  F luid  l n ye n to  ry 

Ihe  quantity  of  working  fluid  filled  inside  the  heat  pipe  is 
determined  on  the  basis  of  saturating  the  void  volume  available  in 
the  wick  pore  structure.  Due  to  the  variation  in  the  density  of 
liquid  sodium  with  temperature,  the  wick  may  operate  oversalurated 
or  undersaturated  depending  upon  the  operating  temperature.  Ihe 
calculation  of  Ihe  initial  fill  quantity  is  normally  based  on  the 
optimum  operating  temperature  at  which  the  liquid  transport  number 
of  sodium  is  the  highest.  Ihis  temperature  is  700°C  arid  the  density 

3 

is  t)./fi35  g/cm  .  However,  overfilling  is  tolerated  better  than 
underfilling  and  Ihe  recommended  overfilling  is  about  ^0%  of  the 
<  orrec  t  filling. 


Muid  inventory, 
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total  liquid  vapor 
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Substituting  the  physical  dimensions,  one  obtains 


m.  ,  ,  (0.59148  I  i  2.11629  I  „  (  0.056/04  I. 

total  t  k  A 


i  0.59148  L  )/>  t  (  I  .  2668  L,  )p 
r.  i  p  v 


lor  I  25  cm,  lR  12.5  cm,  I  74.5  cm  and  l.  91  cm. 


m.  ,  .  (14.78/  t-  26.454  <  4.224  i  53.825),-.  ♦  25/. 40/  p 

total  '  "4  'v 


(84) 


or  m 


tota  1 


99 .  29  pA  i  257.40/  P 
4  v 


(85) 


the  total  mass  of  fluid  required  will  vary  depending  upon  the 
temperature  at  which  the  density  value  is  picked.  In  order  to  see 
t  ho  variation,  a  table  of  values  are  calculated  at  various 
tempera  I ures  of  interest  as  shown  in  lahle  12. 

IABII  12.  KUIl)  INVtNlURY'  VARIAI  ION  Wl  lit  lf.MI'1  RAllJKt 


lemperature 

*4 

mp 

pv 

m 

V 

m.  ^  i 
total 

(°C) 

(g  /cm3) 

(g) 

(g  /cm3) 

(g) 

(9) 

97.83 

0.92/0 

92.04 

0 

0 

92.04 

500 

0.8281 

82.22 

0.003x10  3 

0.000// 

82.22 

600 

0.8054 

79.97 

0.013x10  3 

0.0033 

79.97 

/  00 

0.7635 

75.81 

0.050x10  3 

0.0128 

/  5.82 

800 

0.75/3 

/  5 . 1 9 

0. 134x10  3 

0.0344 

/  5.23 

900 

0. / 4  54 

74.01 

0.306x10  3 

0.078/ 

74.09 

1000 

0. /254 

/2 . 02 

0.667x10  3 

0.1/1  / 

72.19 
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Since  the  wick  is  to  he  fully  saturated  with  liquid  sodium  just 
before  freezing  (9/.Hl°C),  the  mass  of  liuqid  to  he  filled  should  ho 
determined  at  this  temperature.  If  it  is  done  so,  the  overfilling 
will  be  21.4%  more  than  the  correct  filling  at  /00°C. 

3.6  Noncondensible  Gas  loading 

One  of  the  means  of  facilitating  the  startup  of  a  liquid  metal 
heat  pipe  i  rom  frozen  state  is  the  introduction  of  noncundens  ib  le 
gases  into  the  heat  pipe.  I  he  amount  of  gas  in  the  pipe  in  the 
initial  (cold)  state  must  be  the  minimum  possible,  but  sufficient  to 
achieve  frontal  heating  on  startup.  Excessive  amounts  of  gas  would 
reduce  the  active  condenser  length  by  blocking  the  vapor  I  low  into 
the  condenser  /one.  A  predetermination  of  the  quantity  of  NCG 
charge  is  necessary  in  order  to  judge  the  a<  t ive/inacl ive  condenser 
lengths  at  a  desired  heat  input  rate  or  operating  temperature. 

A  simple  thermodynamic  calculation  is  done  to  obtain  the 

inactive  condenser  length  (or  NCG  slug  length)  as  a  function  of 

initial  NCG  charge  pressure  and  the  heat  pipe  hot  /one  temperature, 

1.  Assume  that  initially  the  working  fluid  is  frozen  within  the 

pores  of  the  wick  and  the  vapor  core  of  volume  is  filled  with 

NCG  al  a  pressure  and  temperature  1^..  After  the  heat  pipe  is 

started  and  allowed  to  attain  a  steady  state  condition  at 

temperature  1.  the  NCG  is  compressed  towards  the  condenser  end  by 

the  hot  vapor,  let  the  NCG  volume  which  is  still  at  the  temperature 

7  occupies  an  unknown  volume  V  .  The  pressure  of  the  vapor, 
c  2 

Pv(l)  and  that  of  the  NCG  are  the  same  for  equilibrium  condition. 


IPO 


for  an  ideal  gas  the  equation  of  state  can  be  written  for 


states  I  and  2  as, 


V, 


1>  V  --  R1 
2  2  2 


I ^  |  I  since  the  gas  temperature  remains  the  same. 


therefore,  P  V 


P, 

IK 


v2  !.  V] 


P .  ( l  ) 

.  1 . c  V! 

P  (f) 


tor  cylindrical  vapor  core  of  uniform  cross  section,  the  volumes  are 
proporl ional  to  the  lengths.  Hence, 


.  .  I  w 

•ncg  lp  J  fa 

v' 


(86) 


where 


NCG 


P. 


P  (I) 

v 


NCG  slug  length  at  any  temperature  I 
lotal  length  of  heat  pipe  (2  m) 

Initial  NCG  charge  pressure  at  temperature  l( 
Vapor  pressure  of  sodium  at  temperature  I 


Figure  33  illustrates  this  model  pictoriully. 
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Figure  33.  Nonrondensible  Gas  Compression  Model 

The  results  oF  Fq.  (86)  are  graphed  as  lengths  of  NCG  slug  versus 

temperature  for  a  number  of  assumed  NCG  charge  pressures  as  shown  in 

Figure  34 .  For  a  typical  operating  temperature  of  1000  K,  the 

non  condensible  gas  slug  is  compressed  to  negligible  lengths  of  less 

than  1  cm  if  P.  is  below  1  torr.  If  P.  is  between  b  and  P0  torr, 
l  l 

the  gas  slug  length  varies  I  ruin  b  to  ?b  nn.  Above  bO  lorr,  the  pipe 
cannot  operate  below  1000  K.  thus,  this  graph  facilitates  the  choice 
of  selecting  the  proper  NCG  charge  pressure. 


0  006 


a  *3aniv«3dW3i 


of  Temperature 


From  the  point  of  view  of  sodium  nuiss  depletion  (.it  I  he 

evaporator  during  startup)  and  the  length  ot  inactive  condenser 

tolerable,  it  is  recommended  to  use  f>.  ().?  Ion  as  the  initial  NC(i 

i 

rhdrge  pressure. 

3.7  Design  Summary 

A  summary  of  the  pertinent  design  and  performance  numbers  are 
listed  in  I  able  13. 


1ARI  I  13.  r.AS  I.OAOtO  LIQUID  ML  I  At  DWAHH  (SINGLL  ARILRY) 
1)1  SIGN  DM  All  SUMMARY 


MA  |J  RIAI 

Outer  tube,  inner  tube,  end  caps  Stainless  steel  304 


Screen  wick 

Fill  valve 

Working  I luid 
Noncondens th le  gas 


Stainless  steel  304 
wire  mesh 
Stainless  steel 
bel lows  valve  by 
NUPRO  with  stellite 
tip  stem 
Sod i urn 
Argon 


OVIRAII  DM  AIIS 

lemperature  range  (operating) 
Initial  temperature 
Heat  transport  capacity 
(circulation  limited) 

Sonic  1 imi tat  ion 


t  f  feet  i  ve  length 
total  length 

I  vaporator/Ad iabat  ic/Conderiser 
Overall  diameter 
Vapor  core  diameter 


400  H'jO°C 

Room  temperature  (?S°C) 
1800  W  at  800° C 

Not  appl icable  in  gas 
loaded  mode  (controls 
transport  upto  S60°C 
in  vacuum  mode,  refer 
figure  3?) 

1 . 388  m 
2 . 03  m 

0.37b  m/0. /4b  m/0.91  m 
2.22  x  10V  m 
1 .21  x  10 


m 


I  ABIE  13.  GAS  IOADHJ  IIQUID  Ml  I  At  OWAHP  (MNGII  AKIIHY) 
DESIGN  OH  All  SUMMARY  ( c  onl  i  nuo.l ) 


WICK 

lype 


Double  wd I  led  construct 
ion  with  screen  and 
groove  wick  for 
evaporator  and 
condenser  while 
adiabatic  is  screen 
less 


Screen  si?e  (evaporator  and  40  x  40  cm  mesh 

and  condenser) 

Number  of  spiral  wraps  4 

Evaporator  reservoir  wick  40  x  40  cm  SS  304 

screen  plug  1.H9  cm 
dia.  and  12.5  cm. 
long  with  poros i ty 


Capi llary  insert 


Artery  grooves:  Evaporator 

Number  24 

Width,  mm  0.79 

Depth,  mm  1.40 


75.25%. 

24x24  cm  mesh; 
0.24x0.24x1 .2/  cm 


1  ran sport 
1 

2.30 

2.38 


Condenser 

24 

0.79 
I  .40 


Vapor  Vents: 

Number  per  groove  10  (cut  on  none 

alternate  grooves) 


Si ze ,  mm  6 . 3  x  0 . 79 

Axial  spac ing,  mm  1.5 


36  (cut  on 
all  grooves; 
staggered) 
12.7  x  0.79 
12.7 


- 
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1ABI f  13.  GAS  lOAUfO  I 1QU1U  MKIAI  OWAHP  (S1NGI.K  ARHRY) 
Of  SIGN  U  1. 1 A 1 1  SUMMARY  (concluded) 

lit  A  i  IRANSPOR!  I  IM1 1  (Horizontal  inode  of  operation) 

2 

Available  ( tip  1 1  la r y  pumping  head  at  800°(.'  1/96.6  N/m  or  Pa 

I  iquid  pressure  drop  in  wicks  at  800°t;  0.9088  Pa  per  W 

Vapor  pressure  drop  in  vapor  core  of  800°C  0.0911  Pa  per  W 

Maximum  capi I lary  transport  limit  at  800°C  1797  W 

r I  U 1 1)  l NVtNfORY  ( Ana  1  y t  i  c a  1  g rade  sod  i  urn) 

Optimum  fill  for  operation  at  1000  K  92. 04  g 

(67.52  g  in  screen 
&  grooves ;  24 . 52  g 
in  reservoir  wick) 

NONCONDf.NSIBI  f  GAS  Argon  pressure  0.05  2.0  lorr 

at  25°C 


12/ 


CHAPltR  IV 


I  XPf  RIMINI  Al  WORK 

4.1  F  ab  r i c  a  t ion  of  Pa r t s 

The  experimental  heat  pipe  assembly  consists  of  an  outer  lube, 
an  inner  tube,  two  screen  tubes,  a  reservoir  screen  plug,  two  end 
caps,  a  capillary  insert  plug,  and  a  process  valve  with  a  fill  tube 
as  shown  in  Figure  35.  The  tube  and  bar  stocks  used  in  fabricating 
the  parts  are  certified  stainless  steel  material  (1304  Alloy) 
supplied  by  Kilsby  Roberts  Company.  Ihe  outer  and  inner  tubes  are 
cut  in  straight  and  continuous  lengths  and  machined  to  form  the 
assembly.  Ihe  inner  tube  contains  24  equally  spaced  longitudinal 
small  rectangular  grooves  each  on  the  evaporator  and  condenser 
sections  with  a  single  large  artery  groove  on  the  center  section 
(adiabatic)  connecting  them.  Ihese  grooves  are  machined  using 
precision  saw  cutters  on  a  horizontal  milling  machine.  the 
rectangular  vent  slots  (as  shown  in  Figure  30)  are  cut  by 
electro  discharge  machining.  Ihe  screen  tubes  and  reservoir  screen 
wick  are  formed  of  stainless  steel  304  wire  mesh  by  carefully 
trimming  along  the  edges  and  rolling  on  the  mandrel.  A  bellows  type 
high  temperature  valve  with  649°C  at  1 100  kPa  ruling  supplied  by 
NUPRO  Company  is  used  to  prucess  the  working  fluid  and  inert  gas. 
Ihe  end  caps,  fill  tube  and  valve  are  tungsten  inert  gas  (MO) 
welded  in  a  dry  box,  circulated  with  argon.  Ihe  evaporator  section 
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mm 


gure  35.  Photograph  of  the  Heat  Pipe  Parts  Laid  Out  Before  Assembly 


of  the  outer  tube  contains  fine  circumferential  threads  on  the  inner 
wall  to  a  length  of  37.5  cm.  The  inner  tube  wall  material  is 

machined  out  selectively  to  reduce  the  thermal  mass.  The  tolerances 
on  the  various  parts  are  such  that  there  will  not  be  any  gap  between 
mating  parts  after  assembly.  A  closeup  view  of  some  of  the  parts 
are  shown  in  Figure  36. 

4.2  Cleaning  of  Parts,  Assembly  and  Welding 

It  is  very  essential  to  clean  all  the  parts  of  the  heat  pipe 
before  assembling.  Machined  parts  and  tubes  have  contaminations 
such  as  oil,  grease  and  scales.  Those  are  rough-cleaned  initially 
with  vapor  degreasing  by  trichloroethylene  and  followed  by  metal 

wire  brush  scrubbing  and  thorough  cleaning  using  the  procedure  given 
below. 

(a)  Immerse  in  stainless  steel  deoxidizer  solution  (50% 
dilute)  for  20  minutes; 

(b)  Rinse  in  deionized  water; 

(c)  Dry  with  clean  air; 

(d)  Rinse  with  isopropyl  alcohol  and  dry  with  air. 

For  screen  mesh  cleaning,  ultrasonic  agitation  is  used. 

The  assembly  procedure  consists  of  rolling  the  screen  tubes 
over  the  inner  tube  and  inserting  the  composite  wick  structure  thus 
formed  into  the  outer  tube.  The  inner  edges  of  the  screen  mesh  are 

spot  welded  on  the  inner  tube  at  several  locations  in  order  to  get  a 

good  grip  of  the  screen  while  rolling  to  form  a  tube.  A  number  of 
short  sleeve  tubes  are  slid  over  the  screen  tubes  to  hold  them  in 
place.  Sterile  gloves  are  worn  while  handling  the  cleaned  parts  and 
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FILL  VALVE  ASSEMBLY 


Closeup  Views  of  the  Parts 


the  assembly  is  done  in  a  clean  room.  As  the  inner  tube  composite 
is  inserted  into  the  outer  tube,  the  sleeve  tubes  are  progressively 
slid  along  and  removed  at  the  other  end.  This  way,  the  assembly  is 
accomplished  without  damaging  or  displacing  the  screen  tube.  A 
short  pin  on  the  condenser  end  cap  locks  into  one  of  the  24  grooves 
cut  on  the  inner  tube  and  prevents  rotation  of  the  inner  tube 
composite  with  respect  to  the  outer  tube.  Also,  a  reference  mark 
corresponding  to  the  position  of  the  single  adiabatic  groove  is  made 
on  the  outer  tube  for  locating  the  groove  orientation  after  the  ends 
are  sealed. 

The  end  caps,  reservoir  screen  plug  and  fill  valve  are 
assembled  in  place  and  the  heat  pipe  is  made  ready  for  final  seal 
welding.  These  TIG  weldings  are  very  crucial  and  are  done  within  a 
dry  box  circulated  with  argon  of  purity  <5  ppm  oxygen  and  <0.5% 
water  vapor  in  order  to  avoid  oxidation  or  contamination.  A  number 
of  practice  weldings  were  done  on  samples  similar  to  the  heat  pipe 
end-seal  welding  and  examined  for  weld  penetration  and  grain 
structure  prior  to  the  heat  pipe  welding.  A  weld  penetration  of 
about  1.5  mm  was  obtained  which  is  considered  satisfactory  for  the 
1.65  mm  walled  outer  tube. 

4.3  Vacuum  Baking  and  Sodium  Loading 

The  heat  pipe  thus  assembled  and  welded  is  connected  to  a 
sodium  filling  rig  shown  in  Figure  37  which  also  serves  as  a  baking 
and  oxidizing  setup.  This  rig  consists  of  a  graduated  glass 
cylinder,  a  vacuum  pump,  an  argon  gas  cylinder  and  heating 
tapes/heaters  with  a  suitable  valving  arrangement.  Initially,  the 
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I  igure  3/.  Sodium  lining  Rig  Schematic 
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heat  pipe  is  pumped  down  to  about  10  torr  for  three  day.  and 
then  heat  is  applied  gradually  over  the  entire  heat  pipe,  fill  valve 
and  the  plumbing  connecting  to  the  glass  jar  for  thoroughly  baking 
and  degassing  of  the  internal  surfaces.  Ihe  pipe  is  treated  to  9()()°C 
for  4  6  hours  maintaining  the  vacuum  at  10  ^  torr  or  heller.  Ihis 
process  helps  to  solutioni/e  the  welds  and  also  oxidi/es  Ihe 
external  surface  to  improve  the  emissivity  needed  for  the  condenser 
heat  rejection.  After  cooling  to  room  temperature  and  making  sure 
that  the  glass  jar  is  leak  tight,  the  pipe  is  volved  off  and  Ihe 
glass  jar  is  removed  for  sodium  loading. 

the  sodium  loading  to  the  jar  is  done'  within  Ihe  dry  box  hy 
melting  and  draining  pure  analytical  grade  sodium  from  glass  vials. 
Four  vials  of  50  grams  each  are  transferred  to  the  glass  jar  arid 
sealed.  Ihe  glass  jar  is  taken  out  of  the  dry  box  and  mounted  on 
the  sodium  filling  rig.  Ihe  entire  lines  in  the  rig  including  the 
heat  pipe  are  evacuated  to  10  1  torr  or  better  when  the  sodium  is 
in  frozen  condition.  Now,  the  heat  pipe  is  ready  to  be  tilled  with 
sodium.  Ihe  sodium  jar,  heat  pipe  and  the  lines  including  valves 
are  heated  to  150°C  which  is  ~50°C  higher  than  the  melting  point 
of  sodium.  Ihe  heat  pipe  fill  valve  and  the  vacuum  isolation  valve 
are  shut  off  and  the  line,  up  to  the  valve  is  primed  with  sodium  by 
opening  the  valve  at  the  bottom  of  the  sodium  jar.  Argon  gas 
pressure  is  necessary  at  the  top  of  the  sodium  level  to  push  the 
sodium  into  the  line.  Ihe  meniscus  position  of  sodium  in  the  jar  is 
noted  and  required  volume  of  the  Sodium  liquid  is  made  to  flow  into 
the  heat  pipe  by  slowly  opening  the  fill  valve  and  observing  the 
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drop  In  meniscus  level,  lhe  vacuum  pressure  in  the  pipe  and  the 
gravity  head  available  to  the  sodium  are  sufficient  to  complete  the 
filling  process.  Once  the  filling  is  successfully  completed,  the 
fill  valve  is  shut  off  tight  and  the  valve  handle  is  removed  to 
prevent  accidental  opening  of  the  valve. 

The  filled  heat  pipe  is  removed  from  the  sodium  filling  rig 
and  weighed  to  countercheck  the  exact  quantity  of  sodium  filled  in. 
The  sodium  sticking  to  the  flange  joints  and  lines  is  cleaned  with 
suitable  solvent.  Now,  the  heat  pipe  is  ready  to  be  instrumented 
and  tested. 

4 . 4  lest  Setup  Description 

The  test  setup  consists  of  several  subsystems  which  include 
(1)  the  mechanical  support  and  mounting,  (?)  the  vacuum  chamber  and 
the  vacion  pump,  (3)  the  cooling  fluid  circulation  and  the  constant 
temperature  cool  bath,  (4)  the  electrical  resistance  heaters  and  the 
power  regulating,  controlling  and  measuring  circuit,  (b)  the 
thermocouple  temperature  measuring  circuit  and  data  logging 
instrumentation,  (6)  the  radiation  shields  and  the  condenser  load 
regulating  shutters,  and  (7)  the  NCG  filling  rig.  figure  38  shows 
the  schematic  of  the  test  setup  representing  all  these  subsystems, 
lhe  heat  pipe  is  completely  contained  inside  the  vacuum  chamber  in 
order  to  thwart  the  exposure  of  sodium  to  the  ambient  in  case  of  a 
leak  at  high  temperature  testing.  Photographic  views  of  the  test 
setup  are  shown  in  ( igure  39. 


COOLING  FLUIO 
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Figure  39.  Photographic  Views  of  the  Experimental  Setup 


137 


4.4.]  Mechanical  and  Vacuum  System 

The  heat  pipe  is  supported  at  three  positions  spared  at 
lengths  of  the  evaporator,  adiabatic  and  condenser  end  points  on 
brackets  having  ceramic  supports.  these  brackets  are  supported  from 
two  horizontal,  stepped  down  arbors  welded  to  the  vacuum  chamber 
flange.  Two  20  cm  diameter  stainless  steel  tubular  chambers  each 
122  cm  long  and  bolted  together  encase  the  heat  pipe  assembly  and 

provide  the  vacuum  environment.  A  turbo  molecular  pump  and  vac  ion 
pumping  unit  pump  the  chamber  to  10  ^  torr  or  better  vacuum.  One 
end  of  the  chamber  provides  the  necessary  vacuum  port  and  electrical 
and  thermocouple  feedthrough  ports.  Ihe  other  end  houses  the  rotary 
motion  feedthrough  for  regulating  the  shutter  opening.  A  few  view 
ports  fitted  with  quartz  windows  on  the  side  of  the  chamber  enable 
viewing  the  heat  pipe  assembly.  Ihe  whole  chamber  assembly  is 

mounted  on  a  metal  frame  with  a  fulcrum  at  the  center  point  and 
either  end  has  an  adjusting  screw  rod  for  adjusting  the  level  of  the 

frame  and  the  chamber.  Ihe  level  of  the  heat  pipe  assembly  is 

initially  checked  with  that  of  the  frame  before  assembling  tfie 
chamber.  Ihis  way,  the  level  maintained  on  the  frame  or  chamber 
represents  that  of  the  heat  pipe.  Ihe  inside  wall  of  the  chamber  is 
painted  with  a  special  diffuse  black  paint  having  an  absorptance 
(solar)  value  ^9**.  This  treatment  helps  to  obtain  good 

calorimetric  data  on  the  chamber  wall.  The  vacuum  environment  is 
necessary  to  eliminate  convection  heat  losses  and  to  prevent 
oxidation  of  the  heat  pipe  at  high  tempera  lures . 
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The  vacuum  chamber  wall  is  maintained  at  room  temperature  by 
circulating  chilled  water  from  a  constant  temperature  bath  through 
the  copper  coil  (6.3b  mm  diameter  tubing)  wrapped  around  the  chamber 
at  25 .4  mm  intervals.  Ihe  coil  is  tightly  wrapped  and  soldered  to 
the  chamber  using  soft  solder.  A  continuous  duty  centrifugal  pump 
circulates  water  through  the  coil  up  to  a  maximum  flow  rate  of  5.75 
liters/minute  and  a  rotameter  type  flowmeter  monitors  the  flow 
rate.  Three  thermocouples  mounted  on  the  cooling  coil  at  the 
evaporator  end  of  the  chamber,  the  condenser  beginning  and  the  end 
of  Ihe  chamber  provide  temperature  data  for  the  calorimetry.  The 
chamber  is  insulated  externally  with  fiberfrax  insulation. 

4.4.3  t lec.tr  1 ca I  Heat ing  System 

there  are  three  electrical  resistance  coil  heaters  mounted  on 
the  heat  pipe  as  shown  iri  Figure  38.  Ihe  main  heater  is  a  20  cm 
long,  closely  coiled  sheathed  coil  heater  and  is  mounted  directly 
over  the  26  cm  evaporator  screen  wick.  Ihe  auxiliary  heater  (10  cm 
long  and  similar  to  the  main  heater)  is  mounted  on  the  reservoir 
wick  portion  at  the  starting  end  of  the  evaporator.  Ihe  third 
heater  is  mounted  on  a  linear  trace  from  the  beginning  of  the 
adiabatic  section  to  the  condenser  and  back.  F.ach  heater  is  powered 
from  separate  variacs  and  each  of  them  has  a  power  measuring 
circuit.  Ihe  main  heater  has  a  current  shunt  (standard  0.01  ohm 
resistance)  and  a  signal  conditioner  circuit  to  measure  the  current 
through  the  heater.  Ibis  circuit  also  has  a  relay  to  cut  oft  the 


power  in  case  of  an  over  temperature  alarm  at  the  evaporator. 
Figure  40  shows  the  circuit  diagram. 

4.4.4  Thermocouple  Circuit  and  Oata  Logging 

There  are  28  chromel  alumel  (type  K)  high  temperature 
thermocouples  mounted  internal  to  the  vacuum  chamber  (? 2  on  the  heat 
pipe,  3  on  the  shields,  and  3  on  the  chamber  wall)  as  shown  in 
Figure  38.  In  addition,  there  are  two  other  control  thermocouples 
(channels  3c  and  4c)  attached  on  the  heat  pipe  and  positioned  on 
either  side  of  the  evaporator  main  heater.  These  thermocouple 
junctions  are  formed  by  TIG  welding  and  the  lead  wires  are  insulated 
with  porcelain  tubing.  The  junctions  are  positioned  on  the  small 
dents  created  on  the  heat  pipe  wall  and  spot  welded  by  resistance 
welding.  After  this,  the  junctions  are  covered  with  a  thin  bead  of 
high  temperature  cement  (Sauereisen  Cement  No.  78)  which  is  further 
covered  with  a  radiation  shield.  All  the  thermocouple  lead  wires 
are  run  along  the  heat  pipe  to  the  nearest  strip  connector  block 
from  which  the  extension  leads  are  taken  to  the  data  logger  through 
the  vacuum  glass  to  metal  seal  feedthroughs.  Three  other  type  K 
thermocouples  of  the  sheathed  type  are  used  to  measure  the  cooling 
fluid  temperatures  at  the  inlet,  midpoint  and  the  outlet  locations. 
All  the  thermocouples  are  connected  to  a  multichannel  programmable 
data  logger  (Fluke  Model  2280A). 

4.4.5  Radia t i on_ Shield $  a nd  Shut ter 

Multiple  layers  of  polished  stainless  steel  thin  sheet  metal 
are  used  to  shield  the  radiation  losses  from  the  heat  pipe  at  the 
evaporator  and  adiabatic  sections.  Four  concentric  layers  of 
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main  heater  hookup 
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Figure  40.  Circuit  Diagram  of  the  Heater  Power  Input 


0.26  mm  thick  SS304  foil  separated  by  f» . 3b  mm  (89  mm  diameter  to 
127  mm  diameter)  and  73.7  cm  long  shields  in  two  hemicyl indr ica  1 
halves  cover  the  adiabatic  section.  Hve  layers  of  similar 
construction  in  full  cylindrical  form  cover  the  37. 5  cm  evaporator. 
A  radiation  loss  calculation  based  on  the  shield  emissivity  value  of 
0.08  and  pipe  operating  temperature  of  1000  K  shows  only  5%  loss 
through  the  shield,  compared  to  the  direct  black  body  emission  from 
the  pipe.  The  91  cm  long  condenser  section  is  covered  with  a  set  of 
six  twin  layer  shutter  blades  to  form  a  hexagonal  cage  around  the 
condenser.  The  shutters  are  individually  pivoted  on  either  end  by 
their  central  shafts  supported  from  the  heat  pipe  support  brackets. 
All  these  shafts  are  connected  to  a  circular  disc  with  arc  like 
slots  through  L  shaped  cranks.  Ihe  disc  is  in  turn  connecled  to  a 
rotary  motion  feedthrough  (MFT)  mounted  on  the  chamber  end  flange. 
The  rotation  of  the  MFT  actuates  the  shutter  blades  in  unison.  The 
blades  can  be  positioned  in  fully  open  or  fully  closed  or  any 
intermediate  positions  and  thus  the  condenser  radiation  load  can  be 
varied  as  desired. 

4.4.6  NCG.Fi 1 1 inq  Rig 

The  noncondensible  gas  is  filled  into  the  heat  pipe  through 
the  fill  valve  in  -situ  using  the  filling  rig  while  the  free  end 
section  of  the  chamber  Is  dismantled.  A  schematic  diagram  of  this 
system  is  shown  in  Figure  41. 
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Figure  41.  Noncondensible  Gas  Filling  Rig 


4 .5  I est  Procedure 

4.5.1  System  Checkout  and  Preliminary  Test 

A  thorough  system  check  is  done  on  each  and  every  subsystem 

before  the  vacuum  jacket  is  assembled.  All  the  Internal  parts  are 
checked  for  cleanliness  and  greasy  spots  are  cleaned  with  alcohol. 
Any  outgassing  material  such  as  adhesive  tape  and  plastic  ties  are 
removed.  The  alignment  of  the  heat  pipe  is  checked  with  a  leveling 
instrument  and  adjusted  to  remain  horizontal  with  respect  to  gravity 
vector.  It  is  Important  to  check  the  continuity  of  all  the 

thermocouples  and  the  heaters  (in  order  to  detect  possible  breakage 
during  assembly  or  harnessing  the  wires)  and  to  repair  if 
necessary.  After  the  chamber  jacket  is  assembled,  the  turbo 
molecular  pumping  unit  is  used  to  rough  the  system  until  the  vacion 
pumping  unit  could  take  over  pumping.  If  any  leakage  at  the  vacuum 
joints  is  found,  these  joints  are  tightened  further.  A  vacuum  level 

of  10  k  torr  or  better  is  acceptable  for  this  test.  The  cooling 

fluid  circulation  is  checked  out  by  switching  on  the  centrifugal 
pump  and  circulating  the  coolant  through  the  coll  around  the  vacuum 
chamber.  Any  leak  in  the  line  or  joints  is  easily  detectable  and 
repairable.  The  rotameter  type  flowmeter  is  calibrated  using  stop 
watch  and  measuring  jar.  The  room  temperature  readings  of  all  the 
thermocouples  are  checked  against  a  known  reference  thermometer. 
The  heater  circuit  current  measurement  using  the  shunt  resistance  of 
0.01  ohm  and  the  data  logger  Is  verified  using  analog  meters  in  the 
circuit. 
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Ihe  initial  configuration  of  the  heat  pipe  is  vacuum  mode. 
Initially,  it  is  necessary  to  bake  out  the  heat  pipe  and  the 
internal  parts  of  the  chamber  by  slowly  warming  up  the  pipe.  About 
100  W  on  each  of  the  three  heaters  will  be  sufficient  to  raise  the 
temperature  to  ~500°C  for  the  baking  purpose.  The  power  to 
heaters  may  be  increased  or  decreased  as  necessary  to  maintain  the 
temperature  at  a  constant  level  for  about  8  hours.  In  this  process, 
the  outgassing  rate  of  the  vacuum  parts  Is  accelerated  and  more 
importantly,  the  heat  pipe  "wets  in"  with  the  working  fluid.  The 
main  purpose  of  including  the  trace  heater  is  to  uniformly  warm  up 
the  pipe.  An  added  advantage  of  this  is  to  redistribute  sodium  from 
the  condenser  and  adiabatic  area  to  the  evaporator  section  after  an 
unfavorable  shutdown  or  cooldown  of  the  pipe  from  a  previous  test. 
After  these  tests,  the  vacuum  level  in  the  chamber  should  improve. 
At  this  point,  the  pipe  is  ready  for  performance  tests. 

The  objectives  of  the  performance  test,  transient  as  well  as 
steady  slate  are  limited  to  two  modes  of  the  heat  pipe  internal 

configuration,  namely,  (1)  vacuum  mode  and  (2)  noncondensible  gas 

(NCG)  filled  mode.  Further,  the  tests  are  limited  to  the  level 

condition  (axis  of  the  pipe  is  horizontal)  and  to  the  horizontal 
orientation  (3  o'clock  position  with  respect  to  gravity  vector)  of 
the  artery  groove  for  both  of  the  above  modes  of  internal 

conf igurat ion . 

4.5.2  Vacuum  Mode  Testing 

The  vacuum  mode  represents  that  the  heat  pipe  does  not  contain 
any  NCG  witnin  its  vapor  core  and  the  wick  structure  contains  pure 
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sodium  in  its  pores.  The  objectives  of  this  test  are  twofold. 
Firstly,  the  initial  burn  in  test  for  welting  ami  distribution  of 
the  working  fluid  throughout  the  wick  structure  is  done  soon  after 
the  heat  pipe  is  instrumented.  This  is  followed  with  the  steady 
state  performance  tests  to  verify  the  transport  capacity,  operating 
temperature  and  overall  temperature  difference  (Al)  for  a  given 
power  input  for  a  number  of  input  powers  (250  W,  500  W,  1000  W, 
etc.).  The  steady  state  axial  temperature  plots  and  calorimetric 
input/output  measurements  to  determine  the  losses  are  done.  During 
these  tests,  the  pipe's  priming  ability  and  restart  capabilities 
after  a  cooldown  (frozen  fluid)  period  are  evaluated.  1  he  above 
tests  form  the  basis  for  design  verification. 

Secondly,  the  transient  tests  at  several  power  inputs  are 
performed  to  characterize  the  frontal  nature  of  startup  of  a  sodium 
heat  pipe.  The  frontal  startup  behavior  Is  mapped  by  temperature 
versus  time  plots  using  the  thermocouples  attached  to  the  heat 
pipe.  An  important  aspect  of  this  study/test  is  to  understand  the 
startup  behavior  of  a  heat  pipe  with  long  adiabatic  section  (the 
present  pipe  has  74.5  cm  long  adiabatic  section  without  vapor -liquid 
direct  contact  as  opposed  to  traditional  heat  pipes). 

4.5.3  Noncondens ible  Gas  FJ lied  Mode 

The  heat  pipe,  after  characterizat ion  in  vacuum  mode  as 
described  above,  is  filled  with  a  noncondensible  gas  (argon)  to  an 
initial  predetermined  pressure  (in  the  range  0.05  2.0  Torr  at  room 
temperature)  and  sealed. 
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The  main  objective  of  the  test  in  this  mode  is  to  conduct 
transient  startup  tests  from  frozen  state  and  to  determine: 

(a)  the  rise  time  and  temperature, 

(b)  the  startup  time  and  temperature, 

(c)  the  position  of  the  hot  front  as  a  function  of  time, 

(d)  the  hot  zone  temperature  as  a  function  of  time,  and 

(e)  active/inactive  condenser  lengths  (after  steady  state) 

for  any  given  power  input  to  the  evaporator.  These  startups  are 
attempted  from  liquid  as  well  as  frozen  state  of  the  working  fluid. 
A  trace  heater  attached  on  the  adiabatic  and  condenser  sections 
enables  the  thawing  of  sodium  when  powered.  The  measured  heat  pipe 
temperature  versus  time  data  are  used  to  determine  the  rate  of  hot 
front  propagation  from  the  evaporator  into  its  transport  and 
condenser  sections.  These  experimental  results  are  corroborated 
with  the  results  of  the  theoretical  analysis  done  using  the  binary 
vapor  gas  diffusion  model.  Adjustments  to  theoretically  calculated 
or  assumed  parameters  (such  as  evaporator  and  transport  heat  losses, 
unaccounted  losses  through  support  brackets,  emissivity  values, 
diffusivity  values,  etc.)  and  modeling  errors  are  done  as  necessary. 


147 


CHAPTER  V 


RESULTS  AND  DISCUSSION 

As  described  in  the  test  procedure  section  (Chapter  IV),  the 
2  m  sodium  heat  pipe  was  lirst  tested  In  vacuum  mode  and  then  in  gas 
filled  mode  with  1.35  Torr  argon.  Three  sets  of  test  results  were 
obtained.  Test  Sets  1  and  2  correspond  to  the  gas  filled  mode  tests 
with  the  condenser  shutter  closed  and  partially  open  condition, 
respectively.  Test  Set  3  corresponds  to  the  vacuum  mode  tests.  All 
test  sets  cover  steady  state  and  transient  tests  as  well.  The 
transient  tests  preceded  the  steady  state  tests  for  every  power 
input  level.  The  repeatability  of  results  was  good.  The 
orientation  of  the  pipe  axis  was  horizontal  to  the  gravity  vector, 
and  the  artery  groove  was  at  the  3  o'clock  position  as  shown  in 
Figure  42  in  solid  lines.  Both  frozen  state  startup  and  liquid 
state  startup  tests  were  carried  out  in  the  gas  filled  mode,  lhe 
results  follow. 

5 . 1  Vacuum  Mode  Test  Results 

5.1.1  Steady  State  Test  Results 

The  steady  state  test  data  were  used  lo  obtain  the  axial 
temperature  profiles  and  calorimetric  heat  loss  at  various 
evaporator  power  input  levels.  These  results  were  used  to  determine 
the  heat  transport  performance  of  the  pipe. 
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Figure  42.  Heat  Pipe  Test  Orientation 


Calorimetry:  Ihe  heat  input  (Q^)  to  the  evaporator  was 

split  three  ways  as  the  heat  pipe  transported  energy  along  its 
length.  Eq.  (87)  gives  the  heat  balance  for  steady  conditions. 


<v  Qt  m>a  *  Q0  (8'0 

The  quantities  Qr  ,  Q.,  and  Q  are  the  portions  of  the  heat 

F.  A  o 

lost  radially  from  the  evaporator,  adiabatic  and  condenser  zones, 

respectively.  Thus,  Q  is  the  actual  quantity  of  heat  transported 

o 

by  the  heat  pipe  and  this  quantity  has  to  be  compared  with  the  input 
0.  while  comparing  with  the  theoretical  results.  Q  and  Q 

1  t  n 

are  considered  as  losses  since  these  quantities  are  lost  along  the 
way  when  the  pipe  is  transporting.  Four  thermocouples  mounted  on 
the  coolant  line  at  the  entrance  and  exit  locations  of  each  of  the 
heat  pipe  zones,  monitored  the  temperatures  for  calorimetric 

calculations.  Figure  43  shows  the  percentages  of  heat  loss  in  each 
zone  as  a  function  of  the  heat  input.  At  low  heat  inputs  (<600  W) 
losses  in  the  evaporator  and  adiabatic  zones  were  higher  than  the 

transported  energy  which  was  less  than  50%.  At  high  inputs 
(>600  W) ,  0Q  was  stable  at  55-60%.  Ihe  heat  carried  by  the 

coolant  was  calculated  using  Eq.  (88). 


ycal 


p  C  Al  i  . 
y  p  coolant 


(88) 


Coolant  flow  rate,  V  -  2.3  liters/minute;  specific  gravity,  p  -  1.0; 
and  specific  heat,  Cp  -  4.2  J/g  K  for  water.  Substituting  in 
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HEAT  INPUT,  0;(W) 


Figure  43.  Heat  Loss  Distribution  from  Steady  State  Calorimetric  Test  Data 
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Eq.  (88),  Q  -  161.07  (Al  ,  )  watts.  There  was  good 

cal  coolant 

agreement  between  Q.  (power  measured  in  the  current  shunt  circuit) 
and  Q  ,  for  all  the  steady  state  test  data  reported  in 
Figure  43.  The  errors  were  within  the  instrument  uncertainty 
limits.  Appendix  H  shows  the  list  of  numerical  data  obtained  in  the 
steady  state  tests. 

Axial  Temperature  Profile:  Figure  44  shows  the  axial 

temperature  profiles  for  different  transported  power  levels 

(129-515  W).  The  pipe  operated  near-isothermal  from  end  to -end  at 

Q  -•  467  W  and  515  W.  The  average  temperatures  of  the  evaporator 
o 

and  condenser  for  0q  =  515  W  were  558. 1°C  and  534. 6°C, 

respectively,  and  the  total  AT  was  23.5  C.  As  seen  from 

Figure  44,  at  power  levels  below  Q  -  467  W,  the  hot  front  did  not 

o 

reach  the  condenser  end.  This  indicates  the  presence  of  a  trace 
amount  of  gas  in  the  pipe.  Hence,  the  vacuum  mode  tests  actually 
represent  near  vacuum  mode.  The  trace  amount  of  gas  could  have  been 
trapped  within  the  sodium  during  filling  and  was  not  pumped  out 
subsequent  to  the  filling  of  the  pipe. 

Transport  Capacity:  The  transport  capacity  data  attained 
experimentally  are  plotted  on  the  theoretically  predicted  capillary 
and  sonic  limit  curves  in  Figure  45.  The  experimental  points 

correspond  to  the  average  hot  zone  temperature,  Tu  and  the 

transported  power,  Q  .  The  pipe  operated  along  the  border  of  t fie 

o 

sonic  limit  curve.  High  capacity  and  high  temperature  test  data 
were  not  available  since  the  heater  could  not  be  operated  beyond 
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Figure  44.  Steady  State  Test  Axial  Temperature  Profiles  (Vacuum  Mode) 


1000°C  for  fear  of  melt  down.  Induction  heater  or  tantalum  heater 
could  solve  this  problem. 

5.1.2  Transient  Test  Results 

Frozen  State  Startup:  In  the  vacuum  mode,  the  pipe  was  nearly 

free  of  noncondensible  gas.  The  pipe  was  started  from  frozen  state 

of  the  working  fluid  which  was  at  room  temperature.  Figures  46  and 

47  show  the  axial  temperature  profiles  at  20,  30,  60,  and  150 

minutes  after  the  application  of  the  power  input  for  Q  =-  289.6  W 

o 

and  354.4  W,  respectively.  Figures  48  and  49  show  the  transient 
temperature  profiles  at  a  few  specified  axial  locations  for  the  same 
power  levels  mentioned  above.  It  may  be  observed  that  the  hot  zone 
was  not  isothermal  during  startup. 

5 • 2  Noncondensibl e  Gas  Filled  Mode  Test  Results 
5.2.1  Steady  State  Tests 

Two  sets  (Test  Sets  1  and  2)  of  tests  were  conducted 
corresponding  to  the  NCG  charge  pressure  of  1.35  lorr  at  various 
input  power  levels  up  to  1337  W.  The  calorimetric  heat  loss 
distribution  and  the  heat  transport  capacity  with  the  repeatability 
of  some  test  runs  have  been  shown  in  Figures  43  and  45  (Section 
5.1.1)  along  with  the  vacuum  mode  test  results.  The  numerical  data 
are  listed  in  Appendix  H. 

Axial  Temperature  Profile:  Figure  50  shows  the  axial 

temperature  profiles  for  various  transported  power  levels  (Q  from 

o 

145  to  718  W) .  Ihe  adiabatic  zone  temperatures  are  isothermal  at 
each  power  level  and  the  position  of  the  hot  front  lies  at  different 
axial  locations  in  the  condenser  zone.  The  frontal  nature  of  the 
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Temperature  Profiles  at  Specified  Time  from  0-150  Minutes 
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Figure  47.  Axial  Temperature  Profiles  at  Specified  Time  from  0-150  Minutes 
after  Power  Input  (Q0  *  354.4  W) 


oooi 
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Transient  temperature 


CONDENSER 


Transient  Temperature 


temperature  front  is  evident  from  this  figure.  The  slant  of  the 

front  is  off  vertical  and  the  temperature  of  the  cold  zone  is 

progressively  high  with  U  •  This  signifies  the  dominance  of  axial 

o 

conduction  over  the  diffusion  mechanism  at  the  steady  state 

condition  where  Q  --  0.  It  should  be  noted  that  the  marginal  rise 
v 

in  cold  zone  temperature  will  cause  corresponding  pressure  rise  in 

the  gas  slug  which,  in  turn,  will  affect  the  position  of  the  hot 

front . 

Inactive  Condenser  Length:  lhe  steady  state  test  data  are 
plotted  on  the  operating  temperature  versus  the  inactive  (gas 
slugged)  length  of  the  condenser  graph  as  shown  in  Figure  51  to 
compare  with  the  theoretical  values.  The  experimental  points  lie  in 
the  vicinity  of  the  theoretical  curve  corresponding  to  P.  --  1.0 

lorr.  More  data  points  fell  above  the  1.0  Torr  curve  than  those 
below  it.  This  closely  agrees  with  the  1.35  Torr  initial  gas  charge 

pressure.  The  gas  slug  length  varied  from  121  cm  to  28.3  cm  for  the 

transport  power  (Qq)  variation  of  32.2  W  to  782.8  W  (see  tables  in 
Appendix  H).  it  should  be  noted  that  this  plot  is  not  directly 

dependent  on  the  power  level;  Instead,  it  is  directly  dependent  on 
the  temperature.  For  a  given  power  level  If  the  pipe  is  operated  at 
a  higher  temperature,  the  gas  slug  length  will  be  reduced.  The 
experiment  could  not  be  run  at  high  temperatures  (>850  K)  due  to 
heater  limitation.  The  experimental  data  uncertainty  is  +0.05  m 
in  the  length  determinat ion  and  14.5  K  for  temperature 
measurements . 
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INACTIVE  CONDENSER  LENGTH  (L  -L)  cm 


Hot  Front  Length:  The  steady  slate  test  data  are  also  plotted 
on  the  I.  vs  Qq  graph  as  shown  in  Figure  5 2.  The  experimental 
points  lie  within  the  band  of  -  0.5  and  2.00  lorr  curves  and 

the  mean  curve  appears  to  fall  along  P..  -  1.35  lorr  which  is  the 
actual  charge  pressure  of  the  gas. 

5.2.2  Transien t  Tes t s 

Temperature  Profiles:  Transient  tests  were  conducted  in  order 

to  map  the  startup  behavior  of  the  heat  pipe.  The  temperature 

variations  along  the  length  as  a  function  of  time  were  graphed  for 

various  power  levels.  Figures  53  and  54  show  the  axial  profiles  at 

20,  30,  60,  and  150  minutes  after  the  application  of  power  input  for 

0  -  96.6  W  and  258  W,  respectively.  The  temperature  front  at  20 

o 

minutes  is  in  the  adiabatic  zone  and  the  profile  of  the  front  is 

relatively  « latter  than  that  in  the  condenser  zone.  At  150  minutes, 

the  front  is  completely  in  the  condenser  zone  and  is  more  slanted. 

This  is  the  expected  behavior  and  it  agrees  with  the  hypothesis 

proposed  in  the  analysis  section  (Section  2.3)  very  well.  The  vapor 

diffusion  activity  is  strong  during  the  transient  stale  and  the 

axial  conduction  is  dominant  in  the  steady  state.  Figures  55  and  56 

show  the  transient  temperature  profiles  at  a  few  specified  axial 

locations  for  Q  96.6  W  and  258  W,  respectively.  The  number 
o 

shown  along  each  of  the  curves  represents  the  thermocouple 
location.  Ihe  rise  time  (t^)  for  the  given  power  level  is 
obtained  from  this  plot  by  finding  the  time  (after  power  on)  at 

which  l.C  #4  attained  a  temperature  close  to  the  rise  temperature 
(1^)  within  2  K.  Also,  the  startup  time  (t$)  is  obtained  in  a 
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Figure  53.  Axial  Temperature  Profiles  at  Specified  Time  from  0-150  Minutes  after  Power  Input 
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Figure  55.  Transient  Temperature  Profiles  at  Specified  Axial  Locations  T  vs  t  for  Q0  =  96.6 


similar  procedure  applied  to  l.C.  ft  11  and  counting  the  time  from  the 

rise  time.  The  startup  temperature  corresponds  to  the  temperature 

l.C.  #11  reached  at  the  startup  time.  During  startup,  each  l.C. 

location  registers  a  sharp  rise  in  temperature  in  the  adiabatic  zone 

and  somewhat  slow  rise  in  the  condenser  zone  due  to  the  heat  loss  in 

that  zone.  As  the  beginning  and  the  end  of  the  heating  at  a  l.C. 

location  take  place  over  a  period  of  time,  it  is  hard  to  judge  the 

location  of  the  hot  front  precisely.  This  gives  rise  to  a  small 

deviation  in  comparing  the  experimental  data  with  the  theory.  This 

deviation  Is  explained  later  in  connection  with  the  i.  vs  t  plots. 

The  pipe  started  up  from  the  frozen  state  for  all  power  inputs 

up  to  Q.  =  564  W  (Q  -  306  W)  when  applied  as  a  suddenly  applied 

i  o 

load.  Beyond  this  power  level,  only  gradual  Increments  of  25  W  were 

possible.  In  addition,  the  heater  operating  temperature  limitation 

(I000°C)  restricted  the  heat  input  level. 

Comparison  of  Liquid _ State  and _ Frozen  State  Startup:  The 

liquid  state  startup  Is  faster  and  smoother  than  the  frozen  state 

startup.  For  a  typical  power  input  of  Q.  -  564  W,  these  two 

startups  were  conducted.  The  results  are  shown  in  Figures  57 

through  60.  It  is  evident  from  Figure  59  that  l.C.  if 3  registered 

very  high  temperature  and  the  power  had  to  be  reduced  slightly  to 

continue  the  test.  However,  In  the  liquid  state  startup,  there  is 

no  overshooting  of  temperature  Indicating  smooth  startup  (Figure  60). 

Hot  Zone  Temperature  Variation:  The  experimental  temperature 

variations  of  the  hot  zone  are  plotted  for  Q  =  96. b  W,  257.7  W, 

o 

and  306  W  in  Figures  61,  62,  and  63,  respectively,  lhe  experimental 
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Figure  57.  Comparison  of  Temperature  Profiles  for  Liquid  and  Frozen  State  Startup 
T  vs  L  for  Qc  =  306  m  -  Frozen  State  Startup 
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Figure  58.  Comparison  of  Temperature  Profiles  for  Liquid  and  Frozen  State  Startup 
T  vs  L  for  Q0  =  306  W  -  Liquid  State  Startup 
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Comparison  of  Temperature  Profiles  for  Liquid  and  Frozen  State  Startup 
T  vs  t  for  Q0  =  306  W  -  Liquid  State  Startup 


Figure  61.  Hot  Zone  Temperature  Variation  During  Startup  T  vs  t  for  Q0  *  96.6  w 
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Temperature  Variation  During  Startup  T  vs  t  for  Q0  =  257. 
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data  very  well  match  the  theoret lea  1 ly  calculated  values.  All  the 
points  lie  along  the  P.  -  1.0  to  ?.0  lorr  band  of  the  theoretical 
curves.  In  the  time  axis,  t  =  0  corresponds  to  the  rise  time, 
t..  Ihe  3000  seconds  of  transient  test  are  counted  after  the 

i 

lapse  of  t^ . 

Hot  front  Propagation:  The  transient  experimental  data  for 
plotting  the  L  vs  t  plot  were  processed  from  the  two  minute  scan 
group  data  recorded  in  the  data  logger/computer  during  the  test. 
Ihe  raw  data  available  In  the  form  of  degrees  celcius  for  each 
thermocouple  location  at  two-minute  time  intervals  for  2  to  3  hours 
of  test  were  processed  to  obtain  the  length  versus  time  data.  The 
thermocouple  positions  were  translated  into  axial  lengths  from  the 
evaporator  end.  Table  14  lists  the  input  data  used  in  the  computer 
program.  Since  the  hot  front  travels  at  varying  velocities 
depending  upon  the  input  power  and  inert  gas  charge  pressure,  the 
heating  of  a  particular  spot  along  the  length  to  the  startup 
temperature  will  take  a  finite  amount  of  time.  Ihe  beginning  of  the 
heating  process  is  marked  by  a  sudden  rise  in  temperature 
(designated  as  point  "8"  in  Figure  55)  and  the  ending  of  the  heating 
is  marked  by  more  or  less  constant  temperature  (variation  within 
?  K)  designated  by  point  "E*  In  Figure  55.  In  the  present  work,  the 
ending  point  "E"  is  taken  as  the  reference  to  locate  the  hot  front 
of  the  diffusing  front.  The  "B"  and  "l*  point  designation  applies 
to  all  the  thermocouple  positions  whether  the  slope  of  the 
temperature  profile  is  steep  or  otherwise.  However,  the  point  "E" 
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TABU  14.  THERMOCOUPLE  POSITIONS  MEASURED  EROM  1  HE 
END  OE  I  HE  EVAPORA I  OR 


T.C. 

Number 

Position 

(cm) 

T.C. 

Number 

Position 

(cm) 

1 

0 

11 

109.85 

2 

0.79 

12 

120.24 

3 

12.22 

13 

130.56 

4 

37.62 

14 

140.96 

5 

48.02 

15 

151.35 

6 

58.26 

16 

161.59 

7 

68.58 

17 

171.75 

8 

78.81 

18 

181 .99 

9 

89.21 

19 

192.3 

10 

99.45 

20 

202.2 

will  not  be  applicable  if  a  particular  T.C.  has  not  reached  the 
startup  temperature. 

Figures  64  through  66  show  the  comparison  of  the  experimental 
results  with  the  analytical  predictions  for  Qq  96.6  W,  251.7  W, 
and  306  W,  respectively.  The  agreement  between  the  experimental  and 


analytical  results 

is 

good . 

It 

may  be  observed 

that  the 

experimental  curves 

lie 

along 

a  mean  theoretical  curve  between 

P.  =  1.0  and  2.0  lorr. 

The 

smal  1 

deviation  in  each 

graph  is 

attributed  to  the 

1 imitat ions 

of 

length  determination 

(+0.05  m 

uncertainty)  and  temperature  measurement  (  +  4.5  K  uncertainty). 

In  the  I  vs  t  plots,  t  0  corresponds  to  the  rise  time  (t^) 
and  the  3000  seconds  of  transient  time  is  not  (ountlng  the  rise  time 
duration. 
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Figure  64.  Hot  Zone  Length  Variation  During  Startup  L  vs 


1.6- 
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Figure  65.  Hot  Zone  Length  Variation  During  Startup  L  vs  t  for  Q0  =  257. 


Rjse  TJme, _ Rise_  Temperature^  Startup  time,  and  Startup 

Temperature :  The  experimental  values  for  these  parameters  were 
determined  from  the  output  data  for  every  transient  test.  For  a 
typical  transient  test,  the  rise  time  corresponded  to  the  lime  al 

which  l.C.  #4  reached  the  peak  temperature,  lhe  power  "switih  on" 
time  was  considered  as  t  -  0.  As  explained  earlier,  the  point  "K" 

(Figure  55)  corresponds  to  the  rise  time.  lhe  rise  time  data  for 

various  power  levels  are  listed  in  Table  15.  With  increase  in 

power,  the  rise  time  decreased  as  predicted.  Comparing  these 

results  with  the  predicted  values  in  Table  3,  it  is  found  that  the 

experimental  values  are  almost  twice  that  of  the  theoretical  ones. 

The  reason  for  this  discrepancy  is  the  unaccounted  thermal  mass  (in 
the  form  of  heater  coll,  thermocouples,  clamps,  support,  etc.) 
attached  on  the  evaporator  length. 

The  rise  temperature,  startup  time,  and  startup  temperature 
data  are  also  determined  as  shown  in  Figure  55.  lhe  heat  pipe  is 
supposed  to  have  started  successfully  if  the  hot  front  has  crossed 
T.C.  #11.  These  three  parameters  are  in  good  agreement  with  the 

predicted  values  as  can  be  seen  in  Table  15.  This  table  also  lists 
the  results  for  the  location  of  the  front  in  3000  seconds  after  the 
rise  time  which  match  the  predictions  given  in  lable  3. 

5 . 3  Experiment a  1  Uncerta inty  Data 

The  experimental  data  were  subject  to  the  following 
uncertainty  in  measurement  due  to  systematic  or  instrumental  errors. 
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TASLt  lb.  1RANSUNI  SI  AH  1  UP  USI  RfSIlllS 


SI 

No . 

| 

USI 

Stl  1 

TEST  SE1  2 

SI AR 1  UP  PARAML 1 ERS 

P1  -  ’• 

35  lorr 

L  v  > 

.35  Torr 

°o  (W) 

1  STARTUP 

PARAML  UR 

Q0  (W) 

RESULTS 

32.2 

3000 

1 

Rise  lime. 

96.6 

2880 

2640 

96.6 

ti  (sec) 

145 

1920 

257 

1200 

1200 

241 

306 

960 

32.2 

. 

2 

Startup  Time 

96.6 

3000 

3840 

96.6 

t.  (sec) 

145 

1200 

257 

840 

960 

241 

306 

840 

32.2 

69/ 

3 

Rise  temperature 

96.6 

726 

726 

96.6 

’hi  <*> 

145 

732 

25/ 

732 

731.5 

241 

306 

840 

32.2 

4 

Startup  Temperature 

96.6 

751 

757 

96.6 

Tlts  <K> 

145 

756 

257 

774 

761  .6 

241 

306 

753.3 

b 

1 ocation  of 

32.2 

59  69 

hot  front  in 

96.6 

99  -110 

99  no 

96.6 

3000  seconds  after 

145 

110-120 

rise  time 

257 

141 

130-14) 

241 

L  (cm) 

306 

141  151 

NQTT :  1  he  experimental  settings  for  both  lest  Set  1  and  2  were  the 

same  except  for  the  following: 

lest  Set  1:  Condenser  shutters  closed;  zircar  ceramic  shield 
over  evaporator. 

lest  Set  2:  Condenser  shutters  partially  open;  stainless 
steel  radiation  shield  on  the  evaporator. 
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1.  temperature  measurement  using  chromel  alumel  thermo< oup les : 
0  1250°C  ...  <2.2"C  or  <0.75%  of  the  reading  whichever 


is  the  largest. 

2.  flectrical  power  measurement  through  current  shunt  and 
signal  conditioner: 

0  2000  W  ...  <1%  maximum 

3.  Water  flowmeter  (rotometer  type)  calibration: 

0  5  liters/min  ...  <0.25%  maximum 

4.  Pressure  measurement  using  thermocouple  type  ronvertron 
gage  calibrated  in  lorr  of  argon: 

I  999  Torr  ...  <0.5% 

0  999  ml orr  ...  <0.5% 

All  instruments  were  calibrated  before  use  in  the  test  setup.  By 
periodic  checks,  it  was  found  that  the  instruments  maintained  their 
calibration  within  specification. 
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CHAP] ER  VI 


CON Cl  US  1 ONS  AND  RECOMMENDATIONS 

6.1  Conclusions 

Ihe  aspects  of  starting  up  a  liquid  metal  heat  pipe  from 
frozen  state  of  the  working  fluid  have  been  investigated  in  detail, 
both  analytically  and  experimenta  1  ly.  Ihe  idea  of  gas  loading  the 
liquid  metal  heat  pipe  for  easy  startup  was  recognized  and 
explored.  A  vapor-gas  diffusion  model  was  developed  to  predict  the 
transient  startup  behavior  rather  than  the  axial  conduction  model 
used  in  the  low  temperature  variable  conductance  heat  pipe  (VCHP) 
theory.  A  computer  code  to  predict  the  propagation  of  the  startup 
front  much  Taster  than  the  state  of  the  art  codes  in  the  industry 
has  been  developed.  A  typical  run  in  the  CDC  3600  mainframe 
computer  took  only  220  epu  seconds.  Experimental  results 
corroborate  very  well  with  the  predicted  results. 

Ihe  project  was  carried  out  in  three  phases.  In  the  first 
phase,  called  the  preliminary  design  and  analysis  phase,  a  rough 
design  of  the  2  m  stainless  steel  sodium  heat  pipe  was  worked  out. 
This  was  followed  by  a  lumped  thermal  capacity  calculation  based  on 
steady  state  energy  balance  to  size  the  length  and  thermal  capacity 
parameters.  The  need  for  a  12.5  cm  long  reservoir  wick  in  the 
evaporator  to  meet  the  fluid  supply  during  the  startup  was 
determined.  The  two  dimensional  quasi  steady  state  binary  vapor  gas 


diffusion  model  determined  the  predominant ly  diffusion  controlled 
energy  transport  mechanism  of  the  vapor  at  the  diffusion  front. 
Axial  wall  conduction  was  neglected.  Ifie  key  point  in  this  model 
was  that  the  hot  vapor  front  did  not  move  until  the  vapor  pressure 
became  equal  to  the  initial  charge  pressure  of  the  inert  gas  filled 
in  the  pipe.  The  one  dimensional  transient  thermal  model  coupled 
the  diffusion  and  the  thermal  problems  because  of  the  diffusive  heat 
transfer  coupling  between  the  hot  and  cold  zones.  A  computer  code 
called  SOOAK1  was  written  to  solve  the  first  order,  nonlinear, 
ordinary  differential  equations  with  a  set  of  known  initial  and 
boundary  conditions.  The  code  determined  the  time  rate  of  change  of 
hot  front  lengths  and  temperature  of  the  hot  front.  lhe  mass 
depletion  of  sodium  from  the  evaporator  by  evaporation  during 
startup  was  also  computed  in  the  code  in  order  to  limit  the 
evaporator  heat  input  from  causing  a  dryout. 

In  phase  two,  called  the  detailed  design  and  fabrication 
phase,  a  scaled  down  single  channel  artery  double  walled  heat  pipe 
with  2.22  cm  outer  diameter  (1.27  cm  vapor  core)  and  74.5  cm  long 
adiabatic  section  was  designed  and  fabricated.  The  pipe  was  capable 
of  transporting  1800  W  at  1000  K  and  capillary  limited  at  high 
temperature  and  sonic  limited  at  low  temperature.  lhe  pipe  was 


filled  with 

93.4 

grams 

of 

sodium  and 

1.35 

lorr  argon  at  room 

temperature . 

The 

pipe 

was 

instrumented 

and 

tested  in  a  vacuum 

chamber  for  easy  calorimetry  and  safety. 
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In  the  final  phase  of  the  study,  the  heat  pipe  was  tested  in 
vacuum  and  gas  filled  modes  and  the  results  were  compared  with  the 
analytical/computat ional  results . 

In  summary,  the  following  major  conclusions  were  reached  from 
this  study: 

(1)  Noncondensible  gas  loading  of  the  liquid  metal  heat  pipe 
helped  the  pipe  to  start  easily  even  from  frozen  state.  Large  power 
inputs  (>600  W),  if  suddenly  applied,  caused  evaporator  dryout. 
The  gas  charge  pressure  could  be  predetermined  to  minimize  the 
inactive  condenser  length. 

(2)  The  long  adiabatic  artery  without  a  fine  capillary  surface  did 
not  pose  any  priming  problem  in  the  presence  of  noncondensible  gas. 

(3)  In  gas  loaded  mode,  the  heat  front  propagation  during  startup 
was  found  to  be  solely  diffusion  controlled,  while  in  the  steady 
state,  axial  conduction  determined  the  temperature  profile  of  the 
front.  The  axial  conduction  rate  accounted  for  approximately  16  W 
which  was  only  2  I0X  of  the  condenser  radiated  power. 

(4)  Ihe  variable  conductance  feature  would  be  an  added  benefit  of 
the  startup  solution. 

(5)  Ihe  theoretical  predictions  of  the  rise  time,  rise 
temperature,  startup  time,  startup  temperature,  heat  front  versus 
time  plot,  and  hot  zone  temperature  versus  time  plot  and  the 
experimental  verification  were  in  good  agreement. 

(6)  The  computat iona  I  time  required  for  the  transient  predictions 
was  small. 
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(7)  (he  unheated  evaporator  lengths  and  the  reservoir  wick  kept 
the  evaporator  end  relatively  cooler  than  the  evaporator  exit. 

(8)  The  temperature  of  the  hot  /one  remained  more  isothermal  in 
the  gas  loaded  mode  than  in  the  vacuum  mode.  As  expected,  the 
sLartup  from  the  liquid  state  of  the  working  fluid  was  smoother 
compared  to  that  from  the  frozen  state. 

(9)  The  steady  state  transport  performance  data  followed  very 
closely  with  the  predicted  sonic  limit  curve. 

(10)  A  completely  gas  free  mode  (perfect  vacuum  mode)  of  the 
present  heat  pipe  could  not  be  achieved  due  to  experimental  setup 
limitations.  Hence,  the  difficulty  of  frozen  state  startup  in  the 
vacuum  mode  was  not  demonstrated  to  contrast  the  easy  gas  loaded 
mode  startup.  The  nicrome  resistance  heating  was  another  drawback 
which  restricted  the  heater  temperature  to  less  than  1000°C. 

6 . ?  Recommendations 

(1)  The  heater  temperature  limitation  could  be  overcome  by  using 
tantalum  or  induction  heating. 

(2)  The  radiation  shield  at  the  adiabatic  section  must  also  be  a 
ceramic  (zircar)  type  for  better  insulation. 

(3)  A  gas  venting  arrangement  similar  to  the  one  used  for  the  gas 
filling  rig  may  be  needed  in  order  to  vent  the  pipe  in  hot  condition 
without  losing  the  working  fluid. 

(A)  The  thermocouples  on  the  pipe  may  have  to  be  placed  at  closer 
intervals  to  determine  the  exact  lengths  of  the  diffusion  zone 
during  startup. 
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(5)  Ihe  thermal  capacities  of  the  healer  coil,  thermocouples,  and 
support /mount ing  brackets  should  be  taken  into  account  in  the 
analytical  calculations.  Also,  the  heat  losses  in  the  transport 
section  should  be  accounted  in  theory. 

(6)  In  order  to  avoid  lengthy  iterative  calculation  procedure  of 
the  design  section,  an  integrated  computational  approach  to  couple 
the  hydrodynamics  of  the  heat  pipe  into  the  transient  analysis  may 
be  adopted. 

(7)  Ihe  possible  follow  on  tests  and  study  should  include  aspects 
such  as  pulse  or  programmed  heat  loads,  artery  orientation,  pipe 
tilt  angle,  NCG  pressure  variation,  variable  radiance,  shut  down 
tests,  condenser  load  variation,  etc. 
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APPENDIX  A 


MASS  0IFFUSIV1 I Y  CALCULATION  FOR  SODIUM-ARGON  MIX  1 URt 
The  mass  diffusivlty  D  for  binary  mixtures  of  nonpolar 

+ 

gases  Is  predictable  within  about  5%  by  kinetic  theory  [51]. 
For  gases  A  and  B  at  low  density,  the  Chapman-Enskog  relation  Is 
given  as 


cDA0  =  2.2646  x  10 


T(1/Ma  +  1 /Mg) 
°AB  nD,AB 


(A.l) 


where 

3 

c  -=  molar  density  of  the  mixture,  g -moles/cm 

2 

Dab  =  mass  diffusivlty  of  the  system  A  B,  cm  /sec 
T  -  temperature  of  the  system,  K 

Ma.  M0  -  molecular  weight  of  A  and  B  respectively,  g/mole 
°AB  =  collision  diameter.  Angstrom  units 

Ag  =  collision  Integral,  dimensionless  function  of 
temperature 


^Nonpolar  gases  at  low  density  have  rigid,  noriattracting, 
spherical  molecules  with  no  charge.  Example:  Ar,  Na,  Li,  K.  Polar 
gases,  on  the  other  hand,  have  positively  and  negatively  charged 
ends  and  have  highly  elongated  molecules.  Example:  H?0,  NH3, 
CH3OII,  NOC1 . 


196 


If  ideal  gas  law  is  applied  to  this  binary  gas  system, 
c  -  P/R!  where  P  -  pressure  in  atm.  and  R  is  the  gas  (onstanl 
82.05  cm^  atm./g  mole  K,  then  fq.  (A.l)  becomes 


D._  ,  0.0018583  7  <1/MA  ‘  ,/HB> 

AB  —  - 

P  °AB  fiD,AB 


(A.?) 


In  order  to  determine  <D  from  tq.  (A.l)  or  (A. 2)  as  a 

AH 

function  of  temperature,  <?,„  and  un  ,u  are  to  be  found  out 

AH  D,AB 

first. 

Q  is  a  dimensionless  function  of  t  fie  temperature  and 

D,  AB 

of  the  intermo lecular  potential  field  for  one  molecule  of  A  and  one 

of  B.  Tabulated  values  of  Q  are  available  as  a  function  of 

D ,  AB 

kT /c.nt  where  *  is  the  Boltzmann's  constant  1.380  x 

AB 
23 

10  Joules/K  and  cAg  is  one  of  the  tennard  Jones  parameters 

(the  other  is  o.D). 

Ad 

The  Leonard  Jones  (L-J)  parameters  «  and  e  could 

AB  AB 

be  estimated  fairly  satisfactorily  for  nonpolar,  nonreacling 
molecular  pairs  by  combining  the  I  J  parameters  of  species  A  and  B 
empirically  using  the  following  relationships: 


°A8  *  <"A  *  V 


'AB 


£acb 


(A. 3) 

(A. 4) 
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Ihc  o  and  t , 

the  collision  diameter  and 

t  tie  characteristic 

energy  of  interaction  between  the  molecules  respectively  are 

in  turn 

estimated  from 

Ihe  properties  of  the  particular 

fluid  (gas) 

at  the 

critical  point. 

c/k  0.7/  T 

c 

(A. 5) 

« 

tor  sodium: 

c  ?.44  (T(/Pr),/3 

(A. 6) 

1 

c 

?S00  K;  Pc  3/0  x  10^  H/m7 

3b'j .  16  atm. 

cNa 

-  0.  II  1  K 
( 

?:i 

=  0.  II  x  ?500  x  1  ,3R  x  10 

?6b6.b  x  10  ?3  J 

and  oNa 

1/3 

?-44  (1C/PC)'/J 

1  /3 

3.44  ( /500/36S .  16) 

• 

4.6331  A 

19/ 

1 

I 

For  argon: 


[Note : 


and 

Now  for 


fable  B 
values  of 


1  .  181  K ;  P  48  atm. 

t  r 

O  3.4IH  A;  t/k  1  ?4  K 


Data  taken  from  table  B  1  of  Reference  SI.] 


r.  =  124  x  « 

Ar 

-  124  x  1  .38  x  10  ?3 
171.1?  x  10  ?3  J 

*  3.418  A 

Ar 

the  binary  mixture  o!  sodium  and  argon 


CAB  cNa.Ar  lNa  cAr 

-  (26S6.S  x  10  73  x  171.1?  x  U)'P3)l/y 

=  6/4.22S  x  10  33  J 


(1 .  o 

AB  Na.Ar 


1/2  (<,Na  *  V1 

1/2  (4.8331  c  3.4)8) 
4 . 0?bb  A 


?  of  the  Reference  [SI]  lists  a  for  various 

u ,  Na .  Ar 

Ar)  which  has  a  range  of  0.3  to  100. 
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Now,  In  the  present  work,  the  operating  temperature  range  is 


300  K  to  I  POO  K.  cD  for  the  sodium  argon  mixture  is 

Na  .Ar¬ 
ea  leu  lated  for  this  temperature  range  interpolal  ing  t  tie  values  given 

in  lable  8  2  of  the  Reference  [51  J.  table  A.l  gives  the  data  thus 

calculated  using  Fq.  (A.l)  and  used  in  the  numerical  computation  of 

the  startup  analysis. 


I  ABl  f  A.l.  MASS  Oil  I  US  1 VI 1  Y  VA1.UFS  FOR  SODIUM  ARGON  GAS  MIXTURt 


1 

(K) 

k1 

cNa . Ar 
(N.D) 

flD,Na.Ar 

(N.D) 

cU..  A 

Na .  Ar 

(g  moles/cm  sec ) 

-  ..... 

- - 

.  ^  - - 

- - - - - 

J00 

0.6140 

1.861 

3.4035  x  10'6 

400 

0.818/ 

1  . 592 

4.5941  x  10  6 

500 

1 .0234 

1.424 

5.7424  x  10'6 

600 

1.2281 

1  . 109 

6.8431  x  10  6 

700 

1  .4328 

1.221 

7.9241  x  10  6 

800 

1  .6374 

1.157 

8.9398  x  10  6 

900 

1.8421 

1  .  10/ 

9.9104  x  10  6 

1000 

2.0468 

1.06/ 

10.8381  x  10  6 

1  100 

2.2515 

1  .033 

1  ,  12  x  10  b 

1200 

2.4562 

1.005 

12.6049  x  10  6 

units:  c  -  3. 

cm3 

2 

D  cm 
sec 
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PROGRAM  S  OD  ART  (  TAPES  ,  ?  A  P  E  3  ,  T  A  P  E  7  ,  TAPE  1  ,  TAPti  ,  TAPE*  ,  INPOT  ,  OUTPUT) 

C 

C  THIS  PROGRAM  NUMERICALLY  ESTIMATES  THE  PATE  AT  WHICH 

C  ENERGY  TRANSPORT  OT  VAPOR  OCCURS  BETWEEN  THE  HOT  AND 

C  COLD  ZONES  OP  THE  PIPE.  THE  VAPOR  TRANSPORT  IS 

C  CALCULATED  USING  A  SIMPLE  MODEL  or  THE  STARTUP  "T 

C  A  RADIATION  -  COOLED  HEAT  PIPE. 

r 

c  THIS  VERSION  WAS  CREATED  ON  FEBRUARY  1 J 8 8 

C  IT  CALCULATES  PV  USING  A  NEW  CURVE  TIT  EQN . 

C  SODIUM  MASS  DEPLETION  DURING  STARTUP  IS 

C  COMPUTED  IN  THIS  VERSION  AS  DELM 


C  SODAR2  IS  THE  CORRECTED  VERSION  OT 

C  SOD AR I  WITH  NEWER  EQUATIONS 

C 

c 

c  •••16  FEB  1988****  THIS  VERSION  STORES 

C  OUTPUT  DATA  IN  TAPES  FOR  PLOTTING  ON  TEKTRONIC3 

C  MACHINE  WITH  PLOTTER  HARD  COPIER 

C  •••15MARCH1988***  SODART3  13  THE  VER3ION  or 

C  SODART  WITH  SHORT  TIME  STEP  AND  PRINT  INTERVAL 

C 

C  • *  *  *  1 8 MARC HI  9  8  8 • • *  *  SODARTC  REVISED  VERSION  OF 

C  SOD ARTS  TO  INCLUDE  CD  AT  AVERAGE  TEMP. 

C  or  HOT  AND  COLD  ZONES. 

C 

C 

C 

C  REFERENCE: 

C 

C 

C 

c 

REAL  L, LE, LP , MV, LI 

DIMENSION  PV(10)  , HFG (10) , CD (10)  , TMP  (10) 

COMMON  /  PON 1  /  RW,  D,  LI,  LP,  PAI,  QKORM 
COMMON  /  PON2  /  C,  MV,  TC,TS 
COMMON  /  PON3  /  GAMMA,  EPSILON,  SIGMA 
COMMON  /  PON4  /  BETAI ,  AR,  AO,  PI 

COMMON  /  PONS  /  N3ET,  TMP  ( 1 0 )  , HTG ( 1 0 )  , CD ( I  0 )  ,  P V ( 1 0 ) 
COMMON  /  PON6  /  QHEAT(IO),  QTIMEUO) 

COMMON  /  PON7  /  L I ,  LC , C E , C I , CC , EV SB , THC , THET AS  4 , S 1 , S 2 
C 
C 

REWIND  1 
REWIND  2 
REWIND  4 
REWIND  3 
REWIND  7 


?0 2 


C 

C 

c 


c 

c 

c 


CALL  I  HP  SUB  (  Q,  TFINAL  i 


C 


C 

C 

C 

c 


c 


BETA  I 

-  2.4048255573 

ITRACK 

-  1 

TSTEP 

-  0.001 

PAI 

-4.0  *  ATAN 

TIME 

1 

o 

o 

ICOUNT 

-  1 

TC 

-  300 

IRESET 

-  0 

NUMB 

-  1 

DELM  - 

0  .  0 

NORMALIZE  VARIABLES 

DO  20  .1  -  1,  NSET 

20  TMP  (J)  -  TMP  (J)  /  TC 

PVBAR  -  PI  *  LP  /  (  LP  -  LE  ) 

AAA  -  PAI  *  SIGMA  *  (  TC  *  *  )  ) 
AO  -  0*  EPSILON  *  AAA 

AR  -  2 . 0  *  RW  *  rHC  *  EPSH  •  AAA 

L  -  LE 


CE 


CALL  STARTMP  (  PVBAR,  THETA  ) 

C 

ITOTAL  -  ITHEAT  +  1 

TO-THETA*TC 
RLLE  -  LE  *  RW 
TIM3U-CE«RLLE* (TO-TC) /Q 
WRITE  <  3,  100  )  PVBAR,  TO 

100  rORMAT  (  //10X,  *  INITIAL  VALUE  Or  PV  -  »,ri5.5, 

+  10X,  *  RISE  TEMPERATURE-  »,F1S.5,*  K  * 

WRITE (  3  , 30)  TIM3U 

30  rORMAT  (I0X,  'RISE  TIME  BETORE  STARTUP  -*,ri5  5 
«•*  SECONDS  *,///) 

Q  -  Q  /  QNORM 

QTIME  (ITOTAL)  -  TF INAL 

WRITE (3,444) 

444  FORMAT ( / 3  X ,  “TIME*,  6  X ,  *  H  E  AT  FRONT*,  5X,  *TEMP 

+6X,*G*,10X, *MASS  * , 9X ,  *QV»,  12X,*QO*,I2x,*QR*,/) 

DO  600  IT  -  1,  ITOTAL 
TriNAL  -  QTIME  (IT) 

C  Q  -  QHEAT(IT)  /  QNORM 

ISO  CONTINUE 

IF  (  I  RESET  .EQ.  1000  )  T3TEP  -  T3TEP  *  10. 0 

CALL  RUMGE  (  THETA,  L,  TSTEP ,  TIME,  Q, NUMB) 

IRESET  -  I  RESET  ♦  1 

IF  (  TIME  GE  TriNAL  )  GO  TO  500 

CALL  CALC  (  L,  THETA,  G, QV, HFGBAR,  QR,  QO,  NUMB) 

C 

C  DE-NORMALIZE 


*  BAR*,/ 

) 


?03 


QV  -  JV  *  JNORM 

DELM  -  OELM  *  <  QV/HrGBAR  I  *  TSTEP 

QR  -  QR  *  '.NORM 

0O  -  UO  *  OHORM 
RCALL  -  L  *  RM 

AD  1  L  -  ( LE  +  LI )  •  RW 

TEMP  R  -  THETA  ♦  TC 


IT  (  ICOUNT  . LT .  500  )  GO  TO  165 

ir  <  REALL  .  JT .  AD  I L  >  PELM  -  0.0 

WRITE  l  3  ,  300  )  TIME,  REALL,  TEMPR,  G,DELM,QV,  jo,  OR 

■00  rORMAT  ( 9 ( r 1 0 . 4 , 3X) ) 

WRITE  (  1,888  )  TIME  ,  REALL 

WRITE  (2,888  I  TIME  ,  :IMPR 

WRITE  (  4,888  )  TIME  ,  _V 

WRITE  (7  ,388  )  TIME  ,  jO 

388  rORMAT  <5X,  2ri5.4) 

162  CONTINUE 

ICOUNT  -  0 

165  ICOUNT  -  ICOUNT  +  1 

GO  TO  150 
500  CONTINUE 

ir  ( IT . EQ . ITOTAL)  GO  TO  600 
Q  -  QHE AT (IT)  /  QNORM 
IRESET  -  0 

TSTEP  -  TSTEP  /  10.0 

600  CONTINUE 

WRITE (1, 889) 

WRITE (2,889) 

WRITE (4,889) 

WRITE (7, 889) 

9  89  FORMAT (5 X, •/», //) 

STOP 

END 

C 

C 

SUBROUTINE  INPSUB  (  Q,  TF I NAL  ) 

C 

REAL  L  E , LP , MV ,  LI,  LC 

C 

COMMON  /  PON 1  /  RW,  D,  LE,  LP,  PAI,  QNORM 
COMMON  /  PON2  /  C,MV,TC,TS 
COMMON  /  PON3  /  GAMMA,  EPSILON,  SIGMA 
COMMON  /  PON4  /  BETA 1 ,  AR,  AO,  PI 

COMMON  /  PONS  /  NSET,  TMP(10),  HrO(lO),  CD ( 1 0 ) , P V ( 1 0 ) 
COMMON  /  P  ON  6  /  QHEAT  (10),  QTIME(IO) 

COMMON  /  PON7  /  L I , LC , C E , C I  ,  CC  ,  E P S H ,  r KC , THETAS 4 , S 1 , S 2 

C 

DO  10  I  -  1,10 
QHEAT (I)  -  0.0 

QTIME(I)  -  0.0 
PV ( I )  -0.0 
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10  CONTINUE 


READ 

(5,  *  1 

(ISET 

READ 

(5,  *) 

(  TMP 

(  I  ) 

READ 

(5,  *) 

( HFG  (  X ) 

READ 

(  5  ,  *  ) 

(CD 

1  X  ) 

1-1 . NSET) 
I-i, NSET) 
1-1 , NSET) 


IH 

118 


READ 

<5, 

*  > 

RW 

READ 

(5, 

*  ) 

D 

READ 

(5, 

*  ) 

LE , LI , LC , CE, C I , OC 

READ 

(5, 

*  1 

LP 

READ 

(5, 

*  ) 

MV 

READ 

(5, 

*  ) 

EPS  X  LON , EPSH 

READ 

(5, 

*  ) 

SIGMA,  THC,  TC,  T  3 

READ 

(5, 

*  ) 

PI,Q 

READ 

(5, 

•  1 

ITHEAT 

ir  ( 

ITHEAT 

■  EQ.  ')  )  GO  TO  118 

DO  114  J 

- 

1 , ITHEAT 

READ 

(5, 

*  ) 

QHEAT  ( J) ,  QTIME  ( J ) 

CONTINUE 

CONTINUE 


READ 

<5, 

* 

)  TriNAL 

CALL 

WRITEUP  (  PI,  Q, 

LE  - 

LE 

/ 

RW 

LP  - 

LP 

/ 

RW 

LI  - 

LI 

/ 

RW 

THETAS  - 

TS  /  TC 

QNORM  - 

TC  *  CE  «  RW 

THETAS  4 

- 

THETAS  *•  4 

31  - 

Cl 

/ 

CE 

32  - 

CC 

/ 

CE 

RETURN 

END 


SUBROUTINE  RUNQE  (  TOLD,  LOLD ,  DT,  T, 


Q, NUMB) 


20b 


o  n 


SEAL  1. ,  :  OLD 

riMENSIOH  M2),  VDRV  14,2) 


Ml)  -  TOLD 

'.M2)  -  LOLD 

II  NO  -  .  5  *  DT 

50  100  I  '  1,  4 

GO  TO  110,  20,  20,  20  )  I 

10  THETA  -  TOLD 
L  -  LOLD 

CALL  DERV  I  1,  L.  THETA,  YDRV,  Q , NUMB ) 
GO  TO  100 

20  THETA  -  'Ml)  +  YINC  *  YDRV  11-1,1) 

L  -  M2)  *  YItIC  *  YDRV  11-1,2) 

CALL  DERV  (  I,  L,  THETA,  YDRV,  Q , NUMB) 

CO  TO  100 


30 

100 

110 

120 


THETA  -  Ml) 
L  -  M2) 
CALL  DERV  < 
CONTINUE 


*■  DT  *  YDRV 
*  DT  *  YDRV 
,  L,  THETA, 


11-1,1) 

(1-1,2) 

YDRV,  Q , NUMB) 


DO  120  I  -  1,  2 

SUM  -  Y ( I ) 

DO  110  J  -  1,  4 

T ACTOR  -  DT  /  6.0 

IT  (  J.EQ.2  .OR.  J.EQ.3  )  F ACTOR 
SUM  -  SUM  4- YDRV  <J,I)  *  FACTOR 

Y ( I )  -  SUM 

CONTINUE 
TOLD  -  Y(l) 

LOLD  -  Y ( 2 ) 

T  -  T  +  DT 

RETURN 

END 


DT  / 


3 . 0 


C 

SUBROUTINE  CALC  (  LENS,  TEMP,  <3,  QV,  HFGBAR,  QR,QO,  NUMB) 

C 

REAL  LENG, LP, LE, LI, MV, LC 

COMMON  /  PON  7  /  LI , LC , CE , C I , CC , E P 9 N , T HC , THETAS  4 , 3 1 , S 2 
COMMON  /  PON1  /  RW,  D,  LE,  LP,  PAI,  QNORM 
COMMON  /  PON*  /  BETA 1 ,  AR,  AO,  PI 
COMMON  /  PON2  /  C,MV,TC,T3 
C 

TCN  -  TC/TC 

TAVG  -  (  TEMP  +  TCN  )/2.0 

CALL  INTPOLT  (  TEMP  ,  TAVG,  HFGBAR,  CDBAR,  NUMB) 

CALL  PVCALC  (  TEMP,  PVBAR  ) 
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PVBAR  -  PVBAR  *  1.0E-05 


999 
1  1 1 


EXTRA1  -  LP  -  LENG  -  ( 

LP  *  P  I 

)  / 

PVBAR 

EXTRA2  -  1.339  +  BETA1 

*  EXTRA1 

G  -  0 

.  0 

ir  ( 

ABS ( EXTRA2 1  . LT . 

20  )  G  - 

EXP 

( EXTRA2 ) 

AV  - 

PAI  *  MV  *  CDBAR 

•  HFGBAR 

/ 

(TC  •  CE 

QV  - 

AV  *  G 

QR  - 

AR  *  (TEMP  *  *  4 

-1.0)  * 

LENG 

QO  - 

AO  *  (  LENG  -  LI 

-  LE  )  * 

(  ( 

T  EMP  *  *  4 )  - 

ir  ( 

LENG  .LE.  (  LE  +  LI )  )  QO  -  0. 

0 

RETURN 

WRITE 

<  3,111  ) 

M  VALUE 

or 

LOG  (  PV 

STOP 

END 

20 

30 


SUBROUTINE  DERV  (  I,  LING,  TEMP,  YD,  Q , NUMB ) 

REAL  LENG, LE, LI , LELI 

COMMON  /  PON1  /  RN, D, LE, LP , PAI ,  QNORM 
COMMON  /  PON  7  /  LI, LC,CE,CZ, CC, E  P  S  N , FHC, THETAS  4 , S 1 , S  2 

DIMENSION  YD  (4,2) 


CALL  CALC  (  LENG,  TEMP,  G , QV, HFQBAR  , QR, QO, NUMB  ) 
DENM1  -  LE  ♦  SI  * (LENG  -  LE) 

DENM2  -  LE  ♦  S  1  *  L  I  +  S  2  *  (LENG-LE-LI) 

LELI  -  LE  4  LI 
QOQ  -  Q  -  QV  -  QR 
ir  (  LENG. GT. LELI  )  GO  TO  20 
YD  (1,1)  -  QQQ  /  DENM1 

YD (1,2)  -  QV/ (31* (TEMP-1 . 0) ) 

GO  TO  30 

YD  (1,1)  -  (QQQ-QO) /DENM2 

YD (1,2)  -  QV/ (32* (TEMP-1 . 0) ) 

CONTINUE 

RETURN 

END 


SUBROUTINE  INTPOLT  (  A ,  B ,  C,  D,  NUMB) 

COMMON  /  PONS  /  NSET ,  TMP  (10),  HTG  (10),  CD  ( 1 0 )  , P V ( 1 0 ) 

IT  (  A  .LE.  TMP ( 1 )  )  GO  TO  100 
IT  (  A  .GT.  TMP (NSET)  )  GO  TO  300 

I  -  NUMB 
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5  CONTINUE 

TBEGIN  -  TMP  (I) 

ir  f  I  .EQ-  NSET  )  GO  TO  300 

TEND  -  TMP  <  T  4- 1  > 

IT  I  A  . LE.  TEND  I  GO  TO  10 
I  -  I  +  l 
NUMB  -  I 
GO  TO  5 

10  TDirr  -  TEND  -  TBEGIN 

SLOPE  -  (  A  -  TBEGIN  '  /  TDirr 

SLOPES  -  (  B  -  TBEGIN  )  /  TDIFT 

C  -  HFG(I)  +  <HrG(I  +  ll  -  H  EG  ( I )  )  *  SLOPE 
D  -  CD  (II  +  (CD  (  1  +  1  )  -  CD  ( I )  )  *  S  LOP  E  2 
RETURN 
100  CONTINUE 

A  -  TMP ( 1 ) 

C  -  H  E  G  (  1  ) 

D  -  CD ( 1 ) 

RETURN 
300  CONTINUE 

WRITE  (  3,111  )  A 

111  rORMAT  (  /V  3  X ,  *  * ********  ERROR  # * * * *******  *  *  , 

+  / 3X, *  »  TEMPERATURE  IS  OUT  Or  #  *,  / 

+  3  X ,  *  ♦  RANGE  Or  THE  TABLE  ♦  »,/ 

+  3  X ,  *  ****************************  *,/> 

STOP 

END 

C 

c 

c 

SUBROUTINE  WRITEUP  (  PI,  3,  TFINAL  ) 

COMMON  /  P  ON 1  /  RW,  D,  LE,  LP ,  PAI,  QNORM 

COMMON  /  PON3  /  GAMMA,  EPSILON,  SIGMA 
COMMON  / P  ON  2  / C , MV , TC  ,  TS 

COMMON  /  PONS  /  NSET,  TMP(IO),  HrG(10),  CD ( 1 0 ) , P V ( 1 0 ) 

COMMON  /  P  ON  6  /  QHEAT(IO),  QTIME(IO) 

COMMON  /  P  ON  7  /  L I , LC , C E . C I  ,  CC  ,  EP S K  ,  r HC , THET AS « , S 1 , S 2 

C 

WRITE  (3,111) 

111  rORMAT  (  / / / / 5  X ,  *  -  STARTUP  ANALYSIS  - *,/5X, 

+  * -  TOR  A  HIGH  -  TEMPERATURE - *,/5X, 

+  * -  GAS  -  LOADED  HEAT  PIPE  - *,////) 

C 

c 

WRITE  (3,112) 

112  rORMAT  ( / / / 5  X ,  *  INPUT  DATA  LIST  *,///) 

C 

C 

WRITE  (3  ,  777  )  RW , D ,  L E ,  L I ,  LP , CE , C I , CC , MV , S I GMA 


rORMAT 

(  /  3  X  ,  • 

INSIDE 

RADIUS  Or  WICK 

-  * , T15 . 5 , *  M 

* 

t 

+ 

/  3  X  ,  * 

PIPE  OUTSIDE  DIAMETER 

-  * , r is  .  5 , * 

M  *  , 

♦ 

/  3X,  * 

LENGTH 

Or  EVAPORATOR 

-  • ,  r  1 5  .  5 ,  * 

M  *  , 

+ 

/  3X,  * 

LENGTH 

or  INSULATED  SECTION  -  *,rl5 

.  5  ,  *  M  *  , 

/  3X,  * 

LENGTH 

or  HEAT  PIPE 

-  *,  ns  .5,  * 

M  *  , 

+ 

/3X,  * 

HEAT  CAPACITY  PER  UNIT 

EVAP . LENGTH-* , 

FIS . 5 , *  J/M 
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*  /IX,*  HEAT  CAPACITY  PER  KNIT  ADIA  LENGTH  -*,FI5.5,* 

►  3  X ,  *  HEAT  CAPACITY  PER  UNIT  GOND  LENGTH  -  *  ,  F  1  5 . 5  ,  * 

3  X  ,  *  MOLECULAR  WEIGHT  O  T  VAPOR  -  *  ,  F  1  5 . 5  , 

»  / 3X ,  *  STEFAN- BOLZMANN  CONSTANT  -  *,E15.5,*  W/SQ  CM 

DO  220  I-1,NSET 
TEMP  -  TMP  (  I )  /  TC 
CALL  PVCALC  (TEMP  ,  PVBAR  I 
PVBAR  -  PVBAR  *  1.0E-O5 

PV(I)  -  PVBAR 
220  CONTINUE 
r 
r 

WRITE  (3,222)  TMP , PV, HFG, CD 
222  FORMAT  (/3X,  *TEMP.  X  - * ,  1 0 F  1  2 .  4 / 3 X ,  * P V  BAR  -  *, 

*  l  0  E  1  2  .  4 / 3  X ,  *HFG  J/G  -  * ,  1 0 T 1 2 . 4 / 1 X ,  *  CD  G-M/CM-S-*, 

* 10EI2 . 4///) 

r 

WRITE  (3  ,  3  3  3  )  P I  ,  Q,  EPSILON , EPSH, THC, TC. TS ,  ITHEAT , TFINAL 
333  FORMAT  (3X,*  INITIAL  NCG  CHARGE  PRESSURE  -  *,ri5.5,* 

*  3X ,  *  HEAT  INPUT  TO  EVAPORATOR  -*,F15.5,*  W 

*  / 3  X ,  *  EMI 33  I V I TY  OF  CONDENSER  NALL  -  *,F15.5, 

♦  / 3X, ‘EMISSIVITY  OF  HOT  20NE  INNER  WALL  -*,ri5.5, 

+  / 3  X ,  *  SHAPE  r ACTOR  OF  HOT  TO  COLD  ZONE  »»,F 15.5, 

♦  / 3  X ,  *  AMBIENT  TEMPERATURE  -*,ri5.5,»  K  *, 

♦  /  3X ,  •  COLD  SHROUD  TEMP.  AT  COND .  ZONE  «*,ri5.5,*  K  *, 

+  / 3 X ,  *  NUMBER  OF  HEAT  INCREMENT  STEPS  »*,F15.S, 

♦  / 3X, *  SOLUTION  END  TIME  -  *,FX5.5,*  SECONDS*,///) 

C 

WRITE  <3,  443) 

443  FORMAT  (//10X,  •  OUTPUT  DATA  :  *,//) 

RETURN 

END 

C 

C 

SUBROUTINE  PVCALC  (TIN,PVOUT  ) 

C 

COMMON  /  PON2  /C,MV,TC,T3 

C 

TIN  -  TIN*TC 
T3Q  -TIN  *•  0.5 
TD  -  5367.0/TIN 
COEF  -  2.29E+11 
TD  -  COE F/  ( 1 0 . 0  *  *TD ) 

PVOUT  -  TD/T3Q 
TIN  -  TIN/TC 
RETURN 
END 
C 
C 

c 

SUBROUTINE  STARTMP  (  PVBAR,  THETA  ) 

COMMON  /  P  ON  2  /  C,  MV,  TC,T5 

C 

C 


J/M-K 

J/M-K 

-  K  ~  4  *  , 


BAR  *  , 


?09 


TINC  -  n  .2 
ERRLMT  -  0  . 0  0 1 

TOLD  -  b  .  / 3  . 

ITER  -  0 
10  CONTINUE 

CALL  PVCALC  (  TOLD,  PVOLD  > 

PVOLD  -  PVOLD  *  1.0E-05 

IF  (  PVOLD  . GT .  PVBAR  )  GO  TO  20 

ITER  -  ITER  +  1 

TOLD  -  TOLD  +  TIHC 

GO  TO  10 

20  ERR  -  (  PVBAR-  PVOLD  )  /  PVBAR 

ERR  “  ABS  (  ERR  1 

IT  (  ERR  .LE.  ERRLMT  )  GO  TO  30 
ITER  -ITER  +  1 
TOLD  -  TOLD  -TIHC 

TIHC  -  TINC  /  10.0 

IT  (ITER  -GE.  100)  GO  TO  222 
GO  TO  10 
30  THETA  -  TOLD 
RETURN 

222  WRITE  333 

333  rORMAT  ( / 2X , *  INITIAL  TEMPERATURE  DID  NOT  CONVERGE  •) 
STOP 
EHD 


imuuaat 

10 

300 . 400 . 500 . 600 . 700 . BOO . 900 , 1000 . 1 100 , 1200 
4600.4550.4474.4401,4322.4234.4139.4040.3949.3847 

3 . 4035E-06 , 4 . 594 IE- 6 . 5 . 7424B-6 . 6 . 84 3 IS- 6 . 7 . 9241E-6 , B . 9398E-6 . 9 . 9104B-6 . 
10 . B381E-6 . 11 . 7412E-6 . 12 . 6049E-6 . 

.00635. .022225. .375. . 745 . . 910 . 911 . 717 . 598 . 934 . 886 . 340. 2 . 03 . 23 . . 66 . . 1 . 
S.6696B-B. .01159.300. 300. 66 . 661E- 5 . 250. 0. 3000 
END  OP  PILE 


?10 


APPENDIX  C 


Sflf.CTf.O  PROPtRiY  DA  I A  f-OR  SODIUM 


Atomic  number 

11 

Atomic  weight 

23 

Specific  heal  ratio 

1.6/ 

Melting  point 

3/0.83  K 

Boi ling  point 

1  IS?.?  K 

Critical  temperature 

?S00  K 

Critical  pressure 

370  x  10 

Dens i  ty : 

Solid  ill  ?94  K 
Solid  at  370.83  K 
liquid  at  3/0.83  K 
l  iquid  at  I  ISO  X 


0.9684  g/cm3 
0.9S14  g/tm3 
0.0?/0  g/cm3 
0.740  g/cm3 


Ditrhburn  and  Gi Imour  relation  for  vapor  pressure  f  3 


11  ,  0?>w,„-5Sbb/l 

?.?9  x  10  (1  ) ( 10  ) 


where  py  is 


in  N/m 


arul  I  in  degrees  Kelvin. 


■mper  a  t  ure 

hfq 

(K) 

( N/m^ ) 

(.)/<!) 

300 

0.36/0 

X 

10  8 

4600 

400 

0. 1386 

X 

10  3 

4560 

soo 

0.  76?? 

X 

10  i 

4474 

600 

0 . 49?6 

X 

101 

4401 

700 

0.9648 

X 

10^ 

43?? 

800 

0.8903 

X 

103 

4?34 

900 

0.49/9 

X 

I04 

4139 

1000 

0.  1963 

X 

1 0^ 

4040 

1  100 

0.600! 

X 

I05 

3949 

1?00 

0.161/ 

X 

I06 

38-1  / 

Heat  capacity  (Cp)  at  9/.83°C:  Solid:  1.364?  kJ/kg  K 

f  1 iqu id :  1 .38S0  kJ/kg  K 

(Nearly  constant  in  the  range  ?6  11?6°C) 

Prandtl  Number  (liquid):  0.011  0.004 

( 100  /00°C) 

Maximum  liquid  transport:  ? . 3  x  I013  W/m^  at  bOO’C 

factor  ( ^ 
or  merit  number 


0 


500 


1000  1500 

Twnparatura.  R 


2000 


2500 


Dili  nu<  ( e 
(  I  fleg 


[ensile  Strength  ol  Several  '.olid  Materials. 
R  O.Sbbb  K,  1  kspi  b .  89b  x  106  N/m?) 
[Reproduced  from  deference  4 .  J 
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Maximum  vapor  prwaura. 


AIWNDIX  t 


Dt  SIGN  CHARI  f-OR  HtAI  PIPh  C0N1AINKR  IUBKS 


Design  Chart  for  Heat  Pipe  Container  lubes. 

(I  ps  i  6. 9b  x  10^  N/m? ,  I  kps i  6.89S  x  10^  N/m^) 
[Reproduced  Iron)  Reference  4.] 


Maximum  vapor  pc 


APPENDIX  f 


Df SIGN  CHARI  I  OR  Ml  A I  IM  PL  l  NO  CAPS 


HEAT  PIPE  DESIGN  PROCEDURES 


Design  Chart  for  Heat  Pipe  f.nd  Caps. 

6.895  x  I0-1  N/m'>,  1  kps  i  6.896  x  I0b  N/m2) 
| Reproduced  from  Reference  A.) 


?l  / 


(I  ps  i 


APPENDIX  G 


RESISTANCE  MOOEL  FOR  LIQUID  FLOW  IN  THE  DWAHP 

The  liquid  flow  resistance  network  is  sketched  in  Figure  G.l. 


Figure  G.l 


Definition 


R1  -  Radial  flow  resistance  through  4  layers  of  40  x  40  cm 
screen  wick  in  evaporator. 


Axial 

flow 

resistance 

in 

the  evaporator  screen  tube 

11 

CO 

QC 

Axial 

flow 

resistance 

in 

evaporator  grooves. 

SO 

11 

Axial 

flow 

resistance 

in 

capi 1 lary 

insert. 

R5 

Axial 

flow 

resistance 

in 

adiabatic 

groove. 

V 

Axial 

flow 

resistance 

in 

condenser 

screen  tube. 

R7  " 

Axial 

flow 

resistance 

in 

condenser 

grooves . 

R  =  Radial  flow  resistance  in  condenser  screen  tube. 
8 
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R)  =  FtAr 


Ar  1  n ( d  /d . )  \> 

- 9 — 1_  (_i) 

2ir  Lj-Ar  X. 


With  d 


1 .89  cm,  d . 


1.71  cm,  k 


(6.1) 


1 .9342x10  ^  m^ , 


L_  =  0.33  m  and  Ar  =  0.091  cm 

L 


0  v  2 

R-]  =  2.5343x10  (JL)  N/m  per  watt 


LJ2 
R2  -  -E _ 


ksAs 


ivM 

\ 


(6.2) 


The  cross-sectional  area  (A^)  used  here  is  1  /1 2th  of  the  total 

screen  tube  area  in  order  to  account  for  2  out  of  24  grooves  in  the 

-4  2 

evaporator.  =  0.04311x10  m  . 


12  v„  2 

R2  -  197.878x10  {_!)  N/m  per  watt 


Lc/2  va 

R3  =  -I -  (_!) 

k  A  x 

g  g 


Mean  radius  of  grooves  (Figure  6.2), 


d  4 

r_  =  _o  -  _ 
m  2  2 


=  LIL— 1  OJi  =  0.785  cm 
2 


(6.3) 


8  *1.397  mm 

...... , 

H  r~  < a  *0.794  mm 

Figure  G.2 
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Porosity,  e 


(2*  r  )/l? 
m 


(ass ume  flow  through 
single  groove  only) 


=  0.1931  (19.31%) 


Hydraulic  radius,  rh  a  -  JUL  0.5063  x  10  m 

*  w  e  A 


Permeability,  k  =  2c  rh,g. 

5  \  "'t 


lor  w/A  '  0.568  and  laminar  flow,  faRe^  --  15.5  from  Figure  G.3. 


kg  =  0.6325  x  10“8  m2 


A  =  2ffV  -  0.05742x10  4  m  2 
q  _ 


12 

R3  =  4.5432x10  (_!) 

\ 


N/m  per  watt 


R4  *  _Jl_  (!i) 

Kj  Aj  \ 


(0.4) 


(f  •  Re> 


-  0.0127m,  K j  -  6.7432  x  10  ^  for  24  x  24  cm  screen  and  Aj  = 

-6  2 

Ar  5.67x10  m  for  one  groove. 

G 


1 2  u  2 

R4  -  3.3215  X  10  (!i)  N/m  per  watt 

X 

R5  =  (!i)  (G.5) 

^G  ^ 

-7  2 

i-O  79  m,  K_  -  1.1115  x  10  m  for  2.38  mm  square  groove 
A  G  ’ 

usinq  f  Re  =  14.25  (Figure  G.3)  if  w/6  =  1  and  r  =  <5/2  in  the 
3  ft  ft  h  ,ft 

**6  2 

permeability  equation,  and  A^  =  5.67  x  10  m  for  one  groove 

( Figure  G.4) . 


cu  *  2. 38  mm 


12 

R5  =  1.2528  x  10 


v  2 

(_!)  N/m  per  watt 
X 


Re 


Vi 

K  A. 


V 

(M 

x 
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(G.6) 


The  condenser  screen  and  qroove  geometry  is  the  same  as  that  of  the 

evaporator  except  the  length.  Here,  1  ^  0.91  m,  K  -  1  .934^  x 

-10  ? 

10  m  and  A^  is  same  as  in  . 


545.6635  x  10 


1? 


(  -M 

X 


N/m  per  watt 


«7 


VL 

K  A 

g  g 


(G.7) 


Mean  radius  of  grooves  shown  in  Figure  G.5  (same  as  in  R^) 


8*1.397  mm 


T 

cu  *  0.794  mm 


Figure  G.5 


r  =  0.785  cm 
m 


Porosity,  c  =  - QM _ 

(2wrm)/12 

=■  0.3863 


n  =  2  (Effectively  two 
condenser  grooves 
feed  to  one 
adiabatic  groove) 
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or  38 .63% 


Hydraulic  radius,  rh  ^  --  0.5063  x  10  J  m  (same  as  in  R.^) 


Permeabi 1 i ty 


?  c  rw 


f.  Re» 
l  5. 


-8  2 

1  .  ?65  x  1 0  m  ( twice  as 

that  of  R  ) 


Axial  flow  area  is  same  as  in  R  .  Hence, 


12 

R7  *  b . 2641  x  10  (_!) 

X 


N/m  per  watt 


Ar  1 n(d  /d . ) 
_ o  1 

2 ir  L_  Ar 


(G.8) 


This  is  computed  similar  to  R 


8  «  ? 

Rg  -  0.9190  x  10  (  JL)  N/m  per  watt 

X 


Summari zing, 


8  \1 

Rt  =  2.5343  x  10  (JL) 

X 


197.878  x  101?  Cl) 
'x 


1 2  « 

4.5432  x  10  (_1) 

X 


12  v, 

R4  3.3215  x  10  (_!) 

X 


1 2 

R5  -•  1  .  2528  x  )0  (..A) 

X 


Rft  3.3215  x  10i?  (\) 

X 

1? 

R7  --  6.2641  X  10  (_A) 

X 

Rg  0.910  x  108  Cl) 

X 


(R-)  <ind  Ro  are 
orders  of  magnitude 
small  compared  to 
the  rest  and  hence 

may  he  neglected) 


Since  R  and  R  are  in  parallel,  their  equivalent  resistance  is 

c  0 


calculated  as  R?3  =  R2R3 


R„  +  R„ 
equi  2  3 


Thus 


r2-3 

equi 


12  vn 

4.4412  x  10  {A) 

X 


Simi larly. 


r6-7 

equi 


12 

6.193  x  10  (_*) 

X 


The  total  liquid  flow  resistance  is. 


FtLeff 


=  ♦  R2_3  +  r4  *  r5  +  r6_7  ♦  r8 


(G.9) 


equi 


equi 


1 2  v  2 

i.e.,  RiLeff  =  15.2085  x  10  (I?-.)  N/m  per  watt 
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